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APPENDIX A2. METHOD OF MOMENTS 


Tables A2~-1 ~- A2~-57; 0 < ka ¢ 14.0, h/A = 0.25, 


b/A = 107°, Nel 


Tables A2-58 - A2-66; Experimental monopole antennas, 
h/A * 0.25, b = 0.25 in., N=4 
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APPENDIX A3. OBLATE SPHEROIDAL WAVE FUNCTIONS 
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APPENDIX A4&. VARIATIONAL METHOD 


Page 


Table A4-1; 28.0 < ka <0, h/A= 0.25 128 


127 


TABLE A4-1 , 


2 wali 


28 .2743340 
29.0597324 
29.8451309 
30.6305294 
31 .4159279 
32 .2013245 
32.9867249 
33.7721214 
34 .5575218 
35.3429184 
36.1283150 
36 .9137154 
376991119 
38.4845123 
39.2699089 
40.0553093 
40.8407059 
41.6261024 
42.4115028 
43.1968994 
43.9822998 
44 .7676964 
45 .5530930 
46 .3384933 
47 .1238899 
47 ..9092903 
48.6946869 
49 .4800873 
50.2654839 
51 .0508804 
51.8362808 
52.6216774 
53.4070778 
54.1924744 
54.9778748 
55.7632713 
56.5486679 
57 .33406 83 
58.1194649 
58.9048653 
59 .6902618 
60.4756584 
61.2610588 
62.0464554 


RR. 
(OHMS ) 


-0.0428709 

1 .0697033 
-0 .0340973 
~0 9435477 
-0.0364609 
0.9636334 
-0.0294509 
-0.8572898 
-0 .0314970 
0.8766313 
-0.0257698 
-0.7855343 
-0.0275676 

0.8039876 
-0.0228035 
-0.7249014 
-0 .0243846 

0.7424243 
-0 .0203678 
~0.6729873 
-0 .0217679 

0.6895890 
-0 .0183338 
-0 .6280351 
-0.0195912 
0.6437514 
-0.0166221 
-0 .5887301 
-0.0177491 

0.6036100 
-0.0151622 
-0.5540695 
~0.0161776 

0.5681658 
~0.0139050 
-0.5232750 
-0.0148283 
0.5366415 
-0.0128094 
~0.4957 332 
-0.0136525 
0.5084220 
-~0.0118541 
-0.4709536 


X=- 
00 
(ORMS ) 


-1.0999534 
0.0411037 
0.9679080 
0.0328219 

-0.9881375 
0.0351067 
0.8773347 
0.0284486 

-0.8968816 
0.0304365 
0.8023182 
0.0249676 

-0.8210008 
0.0267142 
0.7391611 
0.0221447 

-0.7569175 
0.0236855 
0.6852531 
0.0198227 

~0.7020830 
0.0211916 
0.6386985 
0.0178772 

-0.6546321 
0.0191055 
0.5980861 
0.0162310 

-0.6131703 
0.0173351 
0.5623448 
0.0148274 

-0.5766323 
0.0158252 
0.5306469 
0.0136123 

-0.5441912 
0.014522] 
0.5023422 
0.0125548 

-0.5151959 
0.0133842 
0.4769124 
0.0116311 


R 
(OHMS ) 


36 .4966393 
37 .6092110 
36 .5054131 
35.5959625 
36.5030479 
37.5031433 
36 .5100594 
35.6822205 
36.5080109 
37.4161415 
36.5137405 
35.7539749 
36.5119400 
37 .3434982 
36.5167046 
35.8146057 
36 .5151253 
37.2819328 
36.5191422 
35.8665199 
36.5177422 
37 ..2290993 
36.5211754 
35.9114723 
36.5199165 
37.1832619 
36.5228882 
35.9507790 
36.5217590 
37.1431198 
36 .5243454 
35 .9854393 
36 .5233307 
37.1076736 
36.5256042 
36 .0162354 
36 .5246811 
37 .0761490 
36.5266991 
36 .0437775 
36 .5258560 
37.0479317 
36 .5276566 
36 .0685539 


128 


x 
(OHMS ) 


20.1576080 
21.2986641 
22.2254696 
21.2903824 
20.2694225 
21.2926674 
22.1348953 
21.2860088 
20.36067 96 
21.2879963 
22 .0598793 
21.2825279 
20 .4365597 
21.2842751 
21.9967213 
21.2797050 
20.5006428 
21.2812462 
21.9428139 
21.2773838 
20.5554771 
21.2787514 
21.8962593 
21.2754383 
20 .6029282 
21.2766666 
21.8556461 
21 .2737923 
20.6443901 
21.2748966 
21 .8199062 
21.2723885 
20.6809292 
21.2733860 
21.7882080 
21.2711735 
20.7133694 
21.2720833 
21.7599030 
21.2701149 
20.7423649 
21.2709446 
21.7344723 
21 2691917 


d( 7/2) 
(NUMERIC) 


0.8214043 
0.7971051 
0.8212068 
0.8421881 
0.8212600 
7993595 
.8211023 
8401521 
8211484 
8012182 
-8210195 
- 8384660 
.8210601 
8027768 
- 8209529 
8370466 
8209884 
8041025 
8208981 
-8358350 
-8209296 
8052436 
8208524 
8347888 
8208807 
8062363 
.8208139 
-8338761 
8208393 
.8071076 
8207812 
-8330729 
8208039 
8078786 
8207529 
- 8323606 
.82077 36 
-8085655 
8207283 
8317245 
-8207472 
8091813 
.8207068 
-8311532 


ooooocoococooococooeocecodcoecoocceoeocoeocoooocoocoooc*oo 


D( 7/2) 
(dBi) 


-0.8544304 
-0.9848440 
-0.8554745 
-0.7459092 
-0.8551931 
-0.97257 86 
-0.8560272 
-0.7564207 
-0.8557835 
-0.9624919 
-0.8564651 
-0.7651453 
-0.8562508 
-0.9540519 
-0.8568177 
-0.7725036 
-0.8566298 
-0.9468860 
-0.8571074 
-0.7787943 
-0.8569410 
-0.9407272 
-0 8573493 
-0.7842339 
-0.8571995 
-0.9353768 
-0.8575531 
-0.7889849 
-0.8574187 
-0.9306855 
-0.8577261 
-0.7931700 
-0.8576057 
-0.9265392 
-0.8578759 
-0.7968851 
-0.8577662 
-0.9228480 
-0.8580059 
-0.8002049 
-0.8579059 
-0.9195414 
-0.8581198 
-0.8031892 


TABLE A4-1 (cont.) 


27a/xr 


98 .1747742 
98 .9601746 
99 .7455673 
100 .5309677 
101.3163681 
102.1017609 
102.8871613 
103 .6725616 
104 .4579620 
105 .2433548 
106 .0287552 
106.8141556 
107 .5995483 
108. 3849487 
109.1703491 
109 .9557495 
110.7411423 
111.5265427 
112.3119431 
113.0973358 
113.8827362 
114.6681366 
115 .4535370 
116 .2389297 
117 ..0243301 
117 .8097305 
118.5951233 
119.3805237 
120.1659241 
120.9513168 
121.7367172 
122.5221176 
123.3075180 
124.0929108 
124.8783112 
125.6637115 
126 .4491043 
127.2345047 
128.0198975 
128.8052979 
129 .5906982 
130.3760986 
131.1614990 
131.9468994 
132.7322998 


R-R 
(oums )~ 


0.3115083 
-0 .0058392 
-0.2945590 
-0.0061717 

0.3017562 
-0.0055717 
-0.2856542 
-0.0058889 

0.2925949 
-0.0053280 
-0.2772731 
-0.0056270 

0.2839723 
-0.0051016 
-0.2693710 
-0.0053838 

0.2758422 
-0.0048909 
-0.2619078 
-0.0051617 

0.2681635 
-0.0046944 
-0.2548479 
-0 .0049507 

0.2608999 
-0.0045108 
-0.2481594 
-0.0047535 

0.25401 86 
-0.0043353 
-0.2418137 
-0.0045691 

0.2474902 
-0.0041742 
-0.2357852 
-0.0043961 

0.2412883 
-0 .0040230 
-0.2300505 
-0.0042375 

0.2353889 
-0 .0038 806 
-0.2245889 
-0.0040811 

0.2297704 


X - Xo 
(OHMS ) 


-0064025 
- 2968726 
0057676 
- 3041370 
.0061040 
. 2878293 
.0055081 
- 2948331 
-0058280 
2793218 
-0052675 
. 2860803 
.0055678 
~2713037 
-0050439 
.2778310 
-0053304 
0.2637343 
0.0048357 
-0.2700431 
0.0051093 
0.2565766 
0.0046414 
-0.2626789 
0.0049032 
0.2497980 
0.0044634 
-0.2557049 
0.0047106 
0.2433694 
0.0042932 
-0.2490909 
0.0045303 
0.2372638 
0.0041335 
-0.2428097 
0.0043575 
0.2314578 
0.0039869 
--0.2368367 
0.0041988 
0.2259296 
0.0038422 
-0.2311501 
0.0040531 


oooocoococooooooqoso 


R 
(OHMS ) 


36 .8510170 
36.5336685 
36 .2449493 
36 .5333366 
36 .8412666 
36.5339355 
36 .2538528 
36.5336189 
36. 8321037 
36.5341797 
36 .2622375 
36 .5338821 
36 .8234825 
36 .5344086 
36.2701378 
36.5341263 
36.8153496 
36 .5346184 
36.2775993 
36 .5343475 
36 .8076706 
36 .5348129 
36 .2846603 
36.5345573 
36 .8004074 
36 .5349998 
36 «2913513 
36 .5347557 
36 .7935257 
36 .5351753 
36 .2976952 
36 5349388 
36 .7869987 
36 .5353355 
36.3037224 
36.5351143 
36 .7807961 
36.5354843 
36 .3094597 
36.5352707 
36.7748985 
36.5356293 
36 .3149185 
36 .5354271 
367692795 


129 


4 
(OHMS ) 


21.2639637 
21.5544338 
21.2633286 
20.9534245 
21.2636642 
21.5453892 
21.2630692 
20.9627285 
21.2633896 
21.5368824 
21.2628288 
20.9714813 
21.2631285 
21.5288639 
21 .2626038 
20.9797306 
21.2628918 
21.5212955 
21.2623959 
20.9875183 
21.2626705 
21.5141373 
21.2622013 
20.9948826 
21.2624645 
21.5073586 
21.2620239 
21.0018559 
21..2622719 
21 .5009308 
21.2618542 
21.0084705 
21 .2620907 
21.4948254 
21.2616940 
21.0147514 
21.2619190 
21.4890194 
21.2615471 
21 .0207233 
21.2617588 
21.4834900 
21.2614021 
21.0264111 
21.2616138 


d( 7/2) 
(NUMERIC) 


0.8135052 
0.8205717 
0.8271082 
0.8205792 
0.8137205 
0.8205657 
0.8269051 
0.8205728 
0.8139229 
0.8205602 
0.8267139 
0.8205669 
0.8141135 
0.8205551 
0.8265338 
0.8205614 
0.8142934 
0.8205504 
0.8263638 
0.8205565 
0.8144633 
0.8205460 
0.8262030 
0.8205518 
0.8146240 
0.8205418 
0.8260506 
0.8205473 
0.8147764 
0.8205379 
0.8259063 
0 .8205432 
0.8149209 
0.8205343 
0.8257692 
0.8205392 
0.8150584 
0.8205310 
0.8256387 
0.8205357 
0.8151891 
0.8205277 
0.8255146 
0.8205322 
0,.8153136 


D( 7/2) 
(dBi) 


-0.8963966 
-0.8588346 
-0.8243767 
-0.8587951 
-0.8952473 
-0.8588665 
-0.8254435 
-0.8588289 
-0.8941671 
-0.8588954 
-0.8264478 
-0.8588602 
-0.8931505 
-0.8589226 
-0.8273939 
-0.8588892 
-0.8921910 
-0.8589475 
-0.8282872 
-0.8569154 
-0.8912850 
-0.8589709 
-0.8291324 
-0.8589402 
-0.8904280 
-0.8589929 
-0.8299333 
-0.8589639 
-0.8896159 
-0.8590138 
-0.8306924 
-0.8589857 
-0.8888455 
-0.8590326 
-0.8314133 
-0.8590065 
-0.8881130 
-0.8590504 
-0.8320996 
-0.8590251 
-0.8874166 
-0.8590677 
-0.8327526 
-0.8590438 
-0 8867530 


TABLE A4-1 (cont. ) 


27a/id 


62.8318558 
63 .6172523 
64 .4026489 
65.1880493 
65.9734497 
66 .7588425 
67 .5442429 
68.3296432 
69.1150436 
69.9004364 
70 .6858368 
71.4712372 
72.2566299 
73.0420303 
73.8274307 
74 .6128311 
75 .3982239 
76.1836243 
76.9690247 
77.7544174 
78.5398178 
793252182 
80.1106186 
80.8960114 
81.6814117 
82 .4668121 
83 .2522049 
84 .0376053 
84 .8230057 
85 .6083984 
86 .3937988 
87.1791992 
87 .9645996 
88.7499924 
89 .5353928 
90 .3207932 
91.1061859 
91 .8915863 
92.6769867 
93 .4623871 
94 .2477798 
95 .0331802 
95 .8185806 
96 .6039734 
97 .3893738 


R-R 
{OHMS ) 


-0.0126236 
0.4830137 
-0.0110089 
-0 .4485402 
-0.0117172 
0.4600173 
-0.0102629 
-0.4281690 
-0.0109139 
0.4391052 
-0.0095979 
-0.4095726 
-0.0102046 
0.4200068 
-0.0090020 
-0.3925286 
-0 .0095630 
0.4024961 
~0.00846 56 
-0.3768499 
-0.0089857 
0.3863834 
-0 .0079807 
-0.3623789 
-0 .0084639 
0.3715079 
-0.0075354 
-0.3489809 
-0.0079905 
0.3577328 
-0.007 1347 
-0.3365408 
-0.0075594 
0.3449399 
-0.0067684 
-0.3249591 
-0.0071708 
0 .3330283 
-0 .0064327 
-0.3141500 
-0.0068096 
0.3219100 
-0.0061239 
-0.3040383 
-0.0064776 


X- X 
(onms ) 


~0 .4891256 
0.0123905 
0 .4539403 
0.0108128 

-0.4655595 
0.0115103 
0.4330855 
0.0100862 

-G.4441538 
0.0107323 
0.4140679 
0.0094376 

-0 .4246245 
0.0100381 
0.3966545 
0.0088557 

-0.4067361 
0.00941 56 
0.3806503 

. 0083369 

3902901 

0088547 

- 3658908 

. 0078620 

~3751190 

0083475 

. 3522363 

0074307 

. 3610804 

.0078812 

3395666 

.0070373 

- 3480521 

0074616 

3277787 

0066775 

3359287 

-0070780 

- 3167838 

-0063474 

3246197 

.0067261 

- 3065044 

0060481 

3140456 


ooeoococooocjoocoocjcococooqocooqcooqcjoo9c;so 


R 
(OHMS ) 


36 .5268860 
37 .0225220 
36 .5284996 
36 .0909691 
36.5277901 
36.9995270 
36 .5292473 
36.1113396 
36 .5285950 
36.9786148 
36.5299110 
36.1299362 
36 .5293045 
36 .95951 46 
36 .5305061 
36 .1469803 
36 .5299454 
36 .9420052 
36 .5310440 
36 .1626587 
36.5305214 
36 .9258919 
36 .5315285 
36.1771317 
36 .5310440 
36.9110184 
36.5319748 
36 .1905289 
365315170 
36 .8972435 
36 .5323753 
36 .2029686 
36 .5319481 
36 .8844490 
36 .5327415 
36 2145500 
36 .5323372 
36.8725357 
36 .5330772 
36 .2253571 
36 .5326996 
36.8614197 
36 .5333862 
36.2354698 
36.5330315 
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X 
( OHMS ) 


20.7684345 
21.2699509 
21.7115002 
21.2683735 
20.7920017 
21.2690716 
21.6906471 
21.2676468 
20.8134079 
21.2682934 
21.6716290 
21. 2669983 
20.8329353 
21 .2675991 
21.6542149 
21.2664165 
20 .8508244 
21 .2669754 
21 .6382103 
21. 2658978 
20 8672714 
21.2664146 
21 .6234512 
21 .2654228 
20 .8824425 
21 .2659073 
21.6097965 
21.2649918 
20 .8964806 
21.2654419 
21.5971279 
21.2645988 
20.9095078 
21.2650223 
21 .5853386 
21.2642384 
20 .9216328 
21 .2646389 
21 5743446 
21 .2639084 
20.9329414 
21.2642860 
21.5640659 
21.2636089 
20 .9435158 


d( 7/2) 
(NUMERIC) 


0.8207241 
0.8097367 
0.8206878 
0.8306370 
0.8207037 
0.8102400 
0.8206710 
0.8301684 
0 .8206857 
0.8106982 
0.8206561 
0.829741) 
0 .8206697 
0.8111171 
0.8206428 
0.8293499 
0 .8206553 
0.8115016 
0.8206307 
0.8289903 
0.8206424 
0.8118557 
0.8206198 
0.8286587 
0.8206307 
0.8121828 
0.8206097 
0.8283519 
0.8206201 
0.8124861 
0.8206007 
0.8280673 
0.8206103 
0.8127679 
0.8205925 
0.8278025 
0.8206016 
0.8130305 
0.8205850 
0.827355) 
0.8205935 
0.8132756 
0.8205780 
0.8273246 
0.8205860 


D( 7/2) 
(dBi) 


-0.8580283 
~-0.9165620 
-0.8582201 
-0.8058872 
-0.8581358 
-0.9138634 
-0.8583090 
-0.8083379 
-0.8582314 
-0.9114082 
-0 .8583879 
-0.8105739 
-0.8583160 
-0.9091644 
~0.8584585 
-0.8126222 
-0 .8583920 
-0.9071064 
-0 .8585226 
-0.8145053 
-0.8584604 
-0.9052117 
-0.8585804 
-0.8162431 
-0.8585226 
-0.9034622 
-0.8586333 
-0.8178512 
-0.8585787 
-0.9018409 
-0.8586810 
-0.8193436 
-0.8586302 
-0.9003348 
-0.8587244 
-0.8207328 
-0.8586765 
-0.8989318 
-0.8587642 
-0.8220285 
-0.8587194 
-0.8976224 
-0.8588012 
-0.8232408 
-0.8587589 


TABLE A4-1 (cont.) 


2m7a/X 


133.5176849 
134.3030853 
135.0884857 
135.8738861 
136 .6592865 
137.4446 869 
138.2300873 
139.0154724 
139.8008728 
140 .5862732 
141.3716736 
142.1570740 
142.9424744 
143.7278748 
144.5132599 
145.2986603 


314.1591797 
314.9445801 


471.2387695 
472.0241699 


628.3183594 
629.1037598 


785 .3981934 
786 .1835938 


942 .4777832 
943.2631836 


10**99 
INFINITY 


R-R 
(oums )~ 


-0.0037433 
-0.2193809 
-0.0039409 
0.2244135 
-0 .0036202 
-0.2144097 
-0 .0038018 
0.2193003 
~0.0034976 
-0.2096590 
~0.0036771 
0.21441 42 
~0.0033878 
-0.2051150 
-0.0035593 
0.2097409 


-0.0011151 
0.0962590 


-0.00061 21 
0.0640954 


-0.0004024 
0.0480331 


-0 .0002733 
0.0384038 


~0.0002102 
0.0319889 


0.0000000 
0 .0000000 


X - X 
(oums )~ 


0.2206601 
0.0037124 
-0.2257294 
.0039089 
.2156310 
.0035835 
.2205569 
.0037729 
. 2108267 
-0034680 
2156154 
0036510 
. 2062321 
.0033526 
2108901 
.0035293 


oooooo°cococeococ.e 


-0.0965002 
0.0010780 


-0 .0642021 
0 .0005776 


-0.0480928 
0 .0003687 


-0 .0384425 
0.0002777 


-0.0320156 
0 .0002084 


0.0000000 
0.0000000 


K 
(OHMS ) 


36 .5357666 
36.3201294 
36 .5355682 
36 .7639236 
36.5358887 
36.3250999 
36.5357056 
36.7588081 
36 .5360107 
36 .3298492 
36.5358315 
36.7539215 
36.5361214 
36 «3343925 
36 .5359497 
36.7492485 


36 .5383949 
36.6357689 


36.5388985 
36 .6036034 


36 .5391083 
36 .5875435 


36 .5392342 
36.5779114 


36.5392990 
36.5714989 


36.5395090 
36 .5395090 
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x 
(OHMS ) 


-4782200 
2612724 
-0318317 
2614689 
~4731922 
2611446 
0370045 
. 2613335 
- 4683876 
- 2610283 
-0419445 
2612114 
4637928 
- 2609138 
.0466709 
- 2610893 


- 1610603 
. 2586384 


. 1933594 
. 2581387 


. 2094669 
. 2579288 


2191181 
. 2578392 


» 2255459 
- 25776 86 


-2575610 
- 2575610 


d( 7/2) 
(NUMERIC) 


0.8205246 
0.825396] 
0.8205290 
0.8154324 
0.8205218 
0.8252832 
0.8205259 
0.8155459 
0.8205191 
0.8251753 
0.8205231 
0.8156543 
0.82051 66 
0.82507 21 
0.8205205 
0.8157580 


0.8204656 
0.8182849 


0.8204542 
0.8190039 


0.8204495 
0.8193634 


0.8204467 
0.8195792 


0.8204452 
0.8197229 


0.8204405 
3.2818500 


D( 7/2) 
(dBi) 


-0.8590841 
~0.8333758 
~0.8590605 
~0.8861202 
~0.8590986 
-0. 8339702 
~0.8590769 
~0.8855162 
-0.8591131 
~0.8345379 
~0.8590918 
~0. 8849385 
-0.8591264 
-0.8350810 
-0.8591059 
-0.8843864 


-0.8593965 
-0.8709548 


-0.8594564 
-0.8671401 


-0.8594813 
-0.8652343 


-0.8594962 
~0. 8640909 


~0.8595040 
-0.8633294 


-0.85952902 
5.16118730 


APPENDIX A5. METHOD OF MOMENTS COMBINED WITH GEOMETRICAL 
THEORY OF DIFFRACTION 


Page 


Table AS-1; 27fa/A= 20.0, h/A = 0.25, b/A=10-° 134 


133 


“£T4G 
SEL01 
£809 
Ate t 
Ovbeed 
64862 
£3188 
cELLb 


ESIBC © 


ECOL 
6TS6vE 
TévO00 
80629 
OLZEZLE 
GOT460 
ELEES 
86709 
1Z90¥ 
eceve 
498TI 
cIgseo 


00000 © 


9ZLZ¢vTO 
60680 
sesctce 
Oc? Iv 
02629 
8461CO 
GEStY 


BvreEGs - 
9T6tS © 


L6Z2ES 
ZZ100 
E9908 
ptocs 


8SSE89 © 
96V09 © 


8466S 
OcZst 
vrBcl 
OZLEEO 


22890 — 


97618 
6S$Z20E 


(GP) 


y3MOd SAILV13Y 


CORRS 
‘790S0 
cellt 
1S¢S6Z8 
6C4TE 
pe4e0l 
GEéFt 
ZzETZS 
90e96 
9PSBEE 
67609 
SZEv0 
I969E 
66029 
E9ZLS6 
9é6rvITCc 
tloeove 
S49 
91S¢08 
OO00ES 
ZScTlTO 
698v0 
ESEE0 
5SE96 
osecs 
4¢79E9 
6~769E 
OZ9c0 
cEEO? 
T2S4h60 
ESS64hv 
22S 18 
TS&900 © 
90202 
$G86e — 
68ve9 O- 
ZeE9GS I- 
6C1ltv e- 
TESET E- 
GL6ZL9 E- 
TOSGS6 E- 
80020 t- 
BLO -E— 
O6ESGe E- 


2 
! 


CANANNETTNNNNNH DODNHNONNNH TTT INNAMWNAKOSO 


(TEP) 
ALIATLOAYIG 


9728S 
6¢0ZLB6 
66001 
GCEZS 
Loose 
St9er 
£90209 
vi 208 
T6226 
B69Gt 
96EvE 
CELEG 
2OsEed 
92vES6 
BSTET 
V2CCE 
LZZ0G6 
G0999 
tv908 
2vZtés 
19266 
S$6S$20 
6422TO 
vpElrd 
2Z0€8 


Tel9o9- 


60Ett 
SZLIgt 
v2988 
9G29¢ 
Srvece 
LESTS 
tlco9 
SSEte 
9TTLO 
66€98 
£8869 
S6EZG6 
2vssr 
L£S82%v 


Sv&6E © 


LE96E 


BIST t © 


CLeLb 


(DT YAWN ) 
ALIAILISAYIAG 


DODO DODCDCO aH RR ANAANMNANMNETTANANANMNANANANK Kener B HOO 


ar 
Le 
98 
ss 
vA 
eg 
28 
18 
0B 
6L 
BL 
Le 
9 
GL 
vl 
EL 
ZL 
TZ 
OZ 
69 
89 
49 
99 
$9 
v9 
E9 
29 
19 
09 
&S 
as 
L¢ 
9 
S¢ 
vs 
ES 
2s 
1S 
OS 
ob 
8b 
Lp 
9p 
Se 


(930) 
STONY 
NOI LYA313 


vOEéS b&- 
tEbvilc <- 
O8el8 G- 
p9IB96& G- 
6E€ce G- 
S9BC9 ¢v- 
CEST!S v- 
SITIO v- 
SG9vEO b- 
SvSOE t- 
26968 t- 
Ev60l G- 
Sttoé 9- 
Ecezy G- 
66LSt¢ OT- 
26808 2t- 
Ssecditv GI- 
coe2s Zi- 
246692 61- 
Tetes 6t- 
BZIGTt O2- 
G9TEE S6I- 
62186 GI- 
PL9ZO vI- 
2sesr ITt- 
EvO09r 6- 
84v00 8- 
L4v660 Z- 
OLZLEZL ‘9- 
Leb6vé 9- 
OZ'28 ‘Z- 
96¢GS 6- 
96629 2I- 
6119S BSI- 
82960 Z2- 
Ov2Z8r 9I- 
822249 It- 
Geert 6- 
2vocrvs L- 
vSS9v L- 
OSTES L- 
9SESGE 6- 
O2@ZEIt St- 
919¢Z2 LI- 
SALINISNI- 
(EP) 
¥3MOd SAILV SY 


sfvceso 2 
S1999 | 
TtTtTZZ O 
$0080 C 
o2G6cs8 0 
L£00ev I 
9€cees | 
cSZEO 2 
porto ¢ 
veOvd I 
SLIGI 1 
Fcé6ee O 
9veG8 O- 
¥IEECE E- 
O€60b t- 
82008 9- 
ZLIiLeY *6- 
SEULE T= 
Sctce ECt- 
cless e€t- 
60E01 tvi- 
96C8S El= 
O9cEes& OT- 
¢08coO 8- 
EtOvey G- 
VEIT €= 
6E966 I- 
8Z£0¢60 I- 
TOé6é89 O- 
8S006 0O- 
COELLZ I- 
28606 ‘€- 
Ss808s 9- 
OGEIS E€I- 
60800 Té- 
TZEEv OI- 
60629 G- 
99E60 E- 
ClI6l t= 
SS9tv [- 
cceSes I- 
S8v0E €- 
TS880 ‘9- 
Z¥ZL<2Z9 ‘Tt= 
SALINIANI- 
(TEP) 
ALIAILOSYIA 


O266S 0 tr 
LEla9 9 Cv 
fEZzes 9 cv 
o9sto | lv 
62602 | Ov 
ai98E 1 of 
s8vres ft se 
S986¢6 I Le 
cO0éS fT 9E 
Sszérv | Se 
GZSTe ft ve 
GZ1so ft 3 2 
82128 0 ce 
ZEzLZE O TE 
OEZ9IE O oe 
26802 0 62 
OTtevIt oO = Pad 
97990 O Le 
E9L~vO O 92 
06070 0 Fad 
S88E0 ‘0 v2 
96970 0 ee? 
89080 ‘0 é2 
Zezgt O Té 
€ZSo82 0 oe 
¢sccr 0 ot 
vELE9 0 et 
60682 0 Zi 
6Z2E¢68 0 91 
2Zz2ts 0 om | 
T8v99 0 vl 
S09tt O et 
vZ6le O at 
Z09S60 0 Tt 
E6200 ‘0 fen § 
06060 ‘0 6 
SSEL2S 'O 8 
6v06t 0 Z 
641199 0 9 
EVT2Z OO S 
St8v9 Oo v 
T2L9v O £ 
219¢Z 0 ro 
96290 0 T 
00000 0 0 
(DT YSWON ) (930) 
ALIAILOSYIC STONY 
NOILVA313 


(SWHO) 20+3v0I90ZLT02 0 f+ 2043820461 0E6E 0 =3ONVG3SdWI LONAGNI 
(HLONSIZAVM) OOTOOOOO 0 
(HLONSTISAYM) COO000G2 0 
(SNIGVYH SNVIdGONNOYS G3SZIWWYON) CO000000 02 =N1Sd3 


=SNIGVY A TdONOW 
=HLSN3S1 3 1dONOW 


LWS3y YTIVOVMY 


[-SV BTAVL 


134 


2£026 
92060 
c842r 
9028S 
£0020 
OOGcSYe 
v686¢6 
rI9b6r 
82262 


EtOGE - 


vZ860 
r9tco 
rL9e8 
TZESE 
£8628 
BEcrs 
Lccvs 


91906 © 
EL026 © 
25066 © 


T6ETT 


629G6E © 
LEGGE 
88289 - 
1ZZE8 © 
£8620 © 
2vBét - 


BG2ZLt 


TOESt © 
8686 © 
90866 — 
2862 - 
16089 
OvO0OS - 


EZ26LS 


98Crl - 
66208 


LEGLI 


verse 


vIécEe © 


62266 © 
66029 © 
68trr © 
vSEtst © 
82266 ° 


rag £04) 


GLO066 © 


(GP) 


y3MOd 3SAILV13uY 


rVI26 [6- 
BSTOO II- 
ELBE ¢- 
Z£EBES E- 
BETeO I- 
TE9Or O- 
€é20sc 0- 
$4060 I- 
40rte E- 
PRILE 9- 
$0060 éI- 
942e0 LE- 
SOBZZ ET- 
eG0coe s- 
BITES &- 
OL962 t- 
Boréc t- 
£pigz G- 
rOeZB L- 
retre tt- 
Eecso 1é- 
o7z0e te- 
B990E ZI- 
étrey 6- 
ecorec F- 
pITEeO 9- 
EZétl 9- 
6SE¢T i- 
Gerert 4- 
66862i E€1- 
BE46re 0F- 
ETTt&e Ge- 
2ece9 vri- 
T£tSte OT- 
6coce B- 
£1260 i- 
OE66/ F- 
699et L- 
9iGEe B- 
spose oT- 
O768068 EI- 
O6TL6 O0€- 
STE6E EE- 
Steet Gt- 
O7e04 OI- 
Ertor 6- 
£0ér6 9- 
(Tap) 
ALIAILIAYNIG 


vOv000 
Ov4szo0 
99682 
OvzdGcs 
Ev06d 
69016 
0018s 
66612 
OBELY 
ELcES 
LEcC9O 
02000 


OStvO © 


viel 
viigve 
SETEE 


WS IEE © 
C9G9S © 
iat he 


vvc90 
Z8200 


L£S200 ° 
62860 © 
S69ET © 
OL60Z © 
6E60zZ © 
B9Zv2 
T6EST 
90221 © 
ZLZ2S60 € 


18600 


OOE00 - 
ZvveO 
21060 © 
OCOCt © 
TIG6T © 


O601Z 


6ZLE6T © 
ZtOGt © 
GLE6O © 
O9TvO © 
£2800 © 
86 700 © 


L90E0 
zttso 


OStet © 
76902 © 
(OT y3SWAN) 
ALIATLIAYIG 


COCCDCCCDOCDCDDDDDDDODDOCOCOCOCCOOOCOOOCCOCOOCCOCOOCCOOOOOCOO 


ost 
62! 
82 
LEI 
9ZLI 
cét 
vit 
elt 
cit 
TZt 
OLt 
691 
Sgt 
L991 
991 
CFI 
vot 
Est 
zat 
tot 
ost 
6¢t 
8ct 
Lot 
9¢61T 
ect 
vst 
est 
cct 
TST 
Ost 
évt 
Srt 
Lvl 
9¢vT 
Stet 
ryt 
cvt 
crt 
tet 
Ort 
6E1 
BET 
Let 
9E1 
Get 
vet 
(934d) 
TTONY 
NOTLVYA313 


close pi- 
$2068 cI- 
£6926 II- 
ccerg ITI- 
19209 T1- 
O2Ze90 eI- 
Ev6lZé 2t- 
L9LZ%o vI- 
60g 9I- 
ESteZs &t- 
66160 ve- 
6t0GE S2- 
ZE966 O2- 
OvItve Lt- 
66982 tvi- 
T9TGS St- 
ZIGO9 Tt- 
vEvtS OI- 
E2268 6- 
EZzOv 6- 
Z1t60 6=- 
6086 8B- 
E672z6 8- 
O86Z0 6- 
EcGre 6- 
OB9cS 6- 
OLTSé 6- 
61Gtv Ol- 
vG826 Ot- 
BL99¥ TI- 
06966 TI- 
66rlLy ZI- 
e6cecs CI- 
ev2s80 Et- 
LZOEt Et- 
12886 2I- 
CvLL9 EI- 
cCeves? Zt- 
997202 TI- 
CB6TT Tt- 
9IStS OT- 
86206 6- 
O660€E “6- 
6LL2L B- 
O9r9t 8- 
ecetez9 L- 
92260 L- 
(GP) 


y3M0d 3AILVTSY 


(SWHO) 


EvlOv B- 
L020B 9- 
B8226 G- 
vG6rtGc G- 
EF6ESGS G- 
26StO 9- 
vwL0E& 9- 
868ZE &- 
eeces ot- 
vIess Et- 
T2EvO BI- 
OSTtOE é6I- 
E9Lvh& vi- 
eLe6e II- 
O6ZLEZ 8- 
26206 9- 
Sr9cse G- 
S9S9IG ‘t- 
vOée8 E- 
wGEGE E- 
Svcev0 E- 
Ové868 ¢- 
v2ZvL8 2- 
2tteé 2- 
wS96lt E- 
tTIt8O0¢6 E- 
tTOEO06 E- 
OS99E “¥- 
666286 tv- 
608Tr ¢- 
teZré G- 
OE9SY “9- 
6tr08 9- 
ELEEO L- 
e¢tso £- 
200v6 9- 
ELB29 “9- 
v9cst 9- 
LLEG9 G- 
Ettzo G- 
Lyo9r v- 
616¢68 E- 
t2t9e% E- 
Ot6L9 “2- 
T6éGTTt SZ 
Ese2ze t- 
LE8v0 “I- 
(TEP) 
ALIATLOISYNIC 


OSer! O Gel 
26902 O vel 
6662 O eet 
L£g062 0 cet 
9EBLS O Tel 
4c06e 0 Oet 
ELz02 0 621 
S2Grvl 0 82t 
09880 0 Lat 
ZBz%r0 O 9Z1 
69610 0 Get 
vZtto oO vzt 
TOZEO 'O ezt 
920140 0 22t 
ezeet o Tet 
vovo2 oO o2t 
02822 ‘0 étt 
6¥v6rE 0 Stitt 
ét2tv oO Ltt 
0029" 0 9TT 
TE96e O ott 
Ttyt¢s O tt 
T6GTS O Ett 
LEeeos 0 Ztt 
TO06ZLe O ttt 
S8crr 0 Ott 
OTLOv O 601 
68S9E 0 B0l 
OTSZE O Lot 
ozzBe 0 90T 
Eevee 0 sOT 
OLLZ2 0 vOt 
ELB0Z 0 cot 
86261 0 cot . 
TESéT O Tot 
OEZOS “O Oot 
CELTS 'O 66 
B90t2 O 86 
E0ZzLz 0 L6 
60TTE O 96 
9SLGE 0 Gé 
ESTIv O vé 
C4tLy O £6 
ZI96ES O cé 
vErvtI9 O Té 
Z2969 0 06 
67GB 0 66 
(1 Y3SWON (930) 
ALIATLOIZBYIA J TONY 
NOI LVA3T3 


ZO+AVOTPOLIOZ 0 f+ 20+38Z0610E6E 0 =39NVA3dWI LNANI 
(HLONATSAVM) OOTOOOOO °O 
(HLONATSAVMN) OCOOCOCOOGS O 
(SNTAVe ANVIdONNOYS G3ZIWWYON) COO00000 02 #=N18d3 


=SNIQve 3 IdONOW 
2HLON3ST 3 1dONOW 


Lwss3y vTwavey 


(*quov) T-<SV eTQPL 


135 


APPENDIX B. 


Bl 
B2 
B3 
B4 
BS 


LISTINGS OF COMPUTER PROGRAMS 


BARDEEN (Integral Equation) 


RICHMD1 and RICHMD2 (Method o€ Moments) 

MONOPL (Oblate Spheroidal Wave Functions) 

MONOSTOR (Variational Method) 

AWADALLA (Method of Moments Combined with Geometric 


Theory of Diffraction) 
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143 
179 
Ley 


201 


APPENDIX Bl. LISTING OF PROGRAM “BARDEEN” ({ntegral Equation) 
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BARDEEN (1 of 3 pages) 


Cc 
Cc APPENDIX A .- PROGRAM "BARDEEN" 
CO eo eke eee ee eee ee eee ake hee keke keke oie ake eke oak keke ekki 
Cc 
THE PRDGRAM "BARDEEN", WRITTEN IN FORTRAN, COMPUTES 
THE RADIATION RESISTANCE AND THE DIRECTIVE 
GAIN FOR A MONOPOLE ELEMENT OF LENGTH, h , MOUNTED 
AT THE CENTER OF AN ELECTRICALLY-SMALL CIRCULAR 
GRDUNDPLANE DF RADIUS, @ , UTILIZING A BARDEEN-MITRE 
INTEGRAL EQUATION METHOD. 


2222202 ££ eee eee eee eee EERE SE SESS SE SERS RES SSR RSS ESSE SSS SSS SS SS 


THE INPUT VARIABLFS OF THIS PROGRAM ARE: 


L = 27h/, (ELEMENT LENGTH IN WAVENUMBERS) 


€ = EPSLN = 277a/A  (GROUNDPLANE RAOIUS IN WAVENUMBERS) 


A = EXCITATION WAVELENGTH 


eee ede oe ete ake kde ake ake eke ke ahead ake ake kde dete dekh teed idee kek ake ako 


THE DUTPUT VARIABLES OF THIS PROGRAM ARE: 

THETAD = @ , ELEVATION ANGLE (DEGREES) 
RADRES = R , RADIATION RESISTANCE (DHMS) 
orrec = d(@), o1RECTIVITY (NUMERIC) 


OIRECO = D(@), DIRECTIVITY (O81) 


* 
* 
we 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 

Wek kk kok kok kkk kaki kkk ko kkk kkk kk kkk kkk ke kkk ke kteake kak keke keke keke ake 


C 
C 
Cc 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C WHERE 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
Cc 
C 


REAL*8 Z,21,22,G,G1,G2,Y,Y1,Y¥2,F4(250),F2(250),T 
REAL*8 1,Q,EPSLN,OCAORE ,A,8.GSUM1,GSUM2,H,HEIGHT,L 
REAL*8 TSUM1,TSUM2, TYSUM1, TYSUM2,BRKT1,YSUM1, YSUM2 
REAL*8 BRKT2,BRKT(89),wW( 12), THTA( 12) ,BRKT3,8RKT4,B8RKT5 
REAL*8 THETA(89),PI,CNVRT,MBRKT(89),DENOM 
REAL*8 RAORES,OIREC(89), THETAD(89) 
REAL*8 DIRECR(89) ,.DIRECD(89),DIRMAX 
PI=3.141592653600 
L=PI/2.0D0 
HEIGHT =L/(2*PI ) 
C*eeeee REE KEKEEKCONVERT DEGREE TD RADI AN * * & ® & & % Hk kk kk 
oo 1 II=1,89 
K=II*2 
THETAO(II )=OFLOAT(K) 
THETACII )=PI* THETAD(II )/ 180.000 
1 CONTINUE 
CHEEK REEEEEEKEEKEKREKEREEKRERKERKEEKREEAEKKEKEEKKEKEKEKKKKRKEEKREK ERE EE RE KK ee 
EPSLN=2.2500 
DD 26 II=1,5 
Cee kk RR RO ROKR ECOMPUTE WITH RESPECT TO EVERY TWO DEGREE******* 
DD 2 INDEX=1,89 
K=4 
T=0.000001D0 
IF (EPSLN.EQ.0.000) I=0.0D0 
IF (EPSLN.EQ.0.000) Q=0.0D0 
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000000 10 
000000206 
00000030 
00000040 
00000050 
O0000060 
00000070 
00000080 
00000030 
00000 100 
00000 1 10 
00000 120 
00000 130 
00000 140 
00000 150 
00000 160 
00000 170 
00000 180 
00000 190 
00000 200 
000002 10 
00000220 
00000230 
00000240 
00000250 
00000260 
00000270 
00000280 
00000290 
00000300 
000003 10 
00000320 
00000330 
00000340 
00000 350 
O0000360 
00000370 
00000380 
000003390 
00000400 
000004 10 
000004 20 
00000430 
000004 40 
00000450 
00000460 
00000470 
000004 80 
00000490 
00000500 
000005 10 
00000520 
00000530 
00000540 
00000550 
00000560 
00000570 


BARDEEN (2 of 3 pages) 


33 


10 


C******CQOMPUTE RAOIATIDN RESISTANCE AND DIRECTIVE GAIN*******444444 04 
44 BRKT3=OSIN(L*DCDS(THETA( INDEX) ) )-DCDS(THETA( INDEX) ) *DSIN(L)- 
*2.000*1/PI 


23 


15 


Ceereereer eee eee ee eWRITE STATMENT S$*# Se Ft eeseeees eee ese eee eee eeee es 


80 


81 


82 


83 


IF (EPSLN.EQ.0.000) GD TD 44 

Z1=DSINH(T ) *DSINH( EPSLN*DCDS(THETA( INDEX) ) ) 
Z22=DSINH(EPSLN) *OSINH(T*OCOS(THETA( INOEX) ) ) 
2=(Z1-22)/OSINH(EPSLN) 
¥1=(H*DSINH(T)*OSIN(L)+T*DCDS(L)*DCDSH(T))/(L*L4+T*T ) 
¥2=0COS(L)*DCOSH(T)/T 

¥=V1-Y2 

G1=(DSINH(T)*OSIN(L)-T*DCDSH(T )*DSIN(L))/(T*T) 
G2=(T*OCDSH(T)*DSIN(L )-H*OSINH(T)*OCDS(L))/(L*L+T*T) 
G=G1+G2 

F1(K)=Z*Y 

F2(K)=Z*G 

T=T+0.0100 

K=K+ 14 

IF(T.LE.EPSLN)GO TO 33 

GSUM120.0D0 

YSUM1=0.0D0 

M=K-1 

DD 10 J=1,M 

Y¥SUM1=YSUM1+F 1( J) 

GSUM1=GSUM1+F2(J) 

CONTINUE 

1=¥YSUM1*0.01 

Q=GSUM1*0.01 


BRKT4=0.000 

BRKT5=0.0D0 

IF (EPSLN.EQ.0.000) GO TO 23 
BRKT4=0SIN(L) *OCOS(THETA( INOEX) ) 


BRKTS=DSIN(L) *ODSINH(EPSLN*DCDS(THETA( INDEX) ))/OSINH(EPSLN) 


BRKT 1*BRKT3+BRKT4-BRKT5 
BRKT2=OCDS(L*DCDS( THETA( INDEX )) )-2.000*Q/PI-OCDS(L) 
BRKT ( INOEX ) =BRKT 1*BRKT 1+BRKT2*BRKT2 

MBRKT ( INDEX ) =BRKT( INOEX ) /OSIN( THETA( INOEX) ) 
CONT INUE 

DENOM=0.0D0 

DD 3 K=1,89 

OENOM=OENDM+MBRKT(K) 

CONTINUE 

DENDM=DENDM* THETA( 1) 
RADRES=14.989625 14D0*DENDM/ (DSIN(L)*DSIN(L) ) 
OIRMAX=0.0D0 

OD 4 K=1,89 


DIREC(K)=2.0D0*BRKT(K)/((DSIN( THETA(K)))**2.000*DENDM) 


DIRECOD(K)=10.00*0LOG10(DIREC(K) ) 

IF (OIRMAX.LT.OIRECO(K)) OIRMAX=DIRECO(K) 
CONTINUE 

00 15 J=1,89 

DIRECR(J)=DIRECO(JU)-OIRMAX 

CONTINUE 


WRITE(6,80) 
FORMAT(1H1,’ INTEGRAL EQUATION METHOD’ ) 
WRITE(6,81) EPSLN 


FORMAT( 1X, “EPSILON= ’,F6.4, ’(NORMALIZEO GROUNDPLANE RAOIUS)’ ) 


WRITE(6,82) HEIGHT 

FORMAT( 1X, ‘MONDOPOLE LENGTH= ‘,F8.6, ‘(WAVELENGTH)’ ) 
WRITE(6,83) 

FORMAT( 1X, “INFINITELY THIN MDNDPOLE RAOIUS’ ) 
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BARDEEN (3 of 3 pages) 


WRITE(6,84) RADRES 
FDRMAT( 1X, ‘RADIATION RESISTANCE = 
WRITE(6,85) 


84 


85 
WRITE(6,86) 

86 
*, RELATIVE PDWER’ ,5X)) 
WRITE(6,87) 


‘,* 10.6,° 


FDRMAT(//,1X, “ELEVATION’ ,53X, ‘ELEVATIDN’ ) 


(DHMS)’ ) 


87 FDORMAT(14X,2(’(DEG)’,6X, ‘(NUMERIC)’ ,7X,’(DBI)’, 11x 


*,.’(DB)’ ,15X)) 
WRITE(6, 88) 

88 FDRMAT(//,2x,’ OO. 
DO 16 L=1,44 


8X, 0.00000’ , 5X, ‘- INFINITY’ ,6X, 


WRITE(6,89) THETAD(L),DIREC(L) ,DIRECD(L) ,DIRECR(L) 
*, THETAD(L+45) ,DIREC(L+45) ,DIRECD(L+45) ,DIRECR(L+45) 
89 FDRMAT(1X,2(1X,F5.1,6X,F9.5,5X,F9.5,6X,F9.5, ##X)) 


16 CONTINUE 


FORMAT (1X,2(’ANGLE’ ,6X, ‘DIRECTIVITY’ ,5X, ‘DIRECTIVITY’ ,5X 


‘~INFINITY’ ) 


WRITE(6,90) THETAD(45) ,DIREC(45S) ,DIRECO(4S) ,DIRECR(45) 
90 FDRMAT(1X,1X,F5.1,6X,F9.5,5X,F9.5,6X,F9.5, 12X, ‘180. 


*,8X,’0.00000’ ,5X, ‘- INFINITY’ ,6X, ’-INFINITY’ ) 


EPSLN=EPSLN+0. 25D0 
26 CONTINUE 

STDP 

END 
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0000 1190 
0000 1200 
0000 12 10 
0000 1220 
0000 1230 
0000 1240 
0000 1250 
0000 1260 
0000 1270 
0000 128D 
0000 1290 
0000 1300 
0000 13 10 
0000 1320 
0000 1330 
0000 1340 
0000 #350 
0000 1360 
0000 1370 
0000 1380 
0000 1390 
0000 1400 
0000 1410 
0000 1420 


APPENDIX B2. LISTING OF PROGRAMS “RICHMDL” and "“RICHMD2"™ 
(Method of Moments) 


Page 
RICHMD1L (Input impedance and current distributions) 144 
RICHMD2 (Far-field radiation pattern) 174 
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RICHMD1 (1 of 30 pages) 


PROGRAM "RICHMD1" 
St HE TE HE HE St TE TE TE TE TE HE TE He He 2G HE HE HE HE EE HE EE HE EEE EEE HEHEHE EE EE 
THIS COMPUTER PROGRAM, IN FORTRAN LANGUAGE, WAS WRITTEN BY 
PROF. RICHMOND OF OHIO STATE UNIVERSITY. IT USES A SINUSOIDAL- 
GALERKIN METHOD OF MOMENTS TO COMPUTE THE INPUT IMPEDANCE AND 
THE CURRENT DISTRIBUTIONS OF A MONOPOLE ELEMENT OF LENGTH 
h AND RADIUS b AT THE CENTER OF A CIRCULAR GROUNDPLANE OF RADIUS a 
A DETAILED DERIVATION IS PUBLISHED IN REFERENCE 1. IN 
ORDER TO ADAPT THE PROGRAM INTO OUR COMPUTER SYSTEM 
(VAX 11/780), THE PROGRAM HAS BEEN CHANGED FROM SINGLE PRECISION TO 
DOUBLE PRECISION. 


UTILIZING THE CURRENT DISTRIBUTIONS DETERMINED FROM THIS 
PROGRAM, A RADIATION PATTERN CAN BE OBTAINED FROM THE MITRE 
COMPUTER PROGRAM "RICHMD2". AN INPUT FILE, “INPUT1. DAT", TO THE 
COMPUTER PROGRAM "RICHMD2" IS AUTOMATICALLY GENERATED FROM THIS 
PROGRAM. 


REFERENCE 
1.) J.H. RICHMOND. "MONOPOLE ANTENNA ON CIRCULAR DISK" 


OHIO STATE UNIVERSITY, ELECTROSCIENCE LABORATORY, 
TECHNICAL REPORT 711639-1, JULY 1979. 


eeegRgkaegeEKRAeKeEKREeEKE KEKE HE HEH HK 


St Hh AG Sb tt 3b HE Se a St St St HE OE HG Se tt He aE OE EE HE HE HE HE He HE EE HE HE HE He HE HE HE HE HE HE HE HE HE HE HE HE TE EE HE HE HE HE HE HE HE HE HE HE He He HH HE Ht Ht tt 


NOTE: THE MAXIMUN NUMBER OF SEGMENTS (NSD+NSW) CAN BE USED IN 
THIS PROGRAM IS 31. IF NUMBER OF SEGMENTS USED IS OVER 
31, THEN ALL THE ARRAYS OF DIMENSION OF 30 AND 30x30 
SHOULD BE CHANGED TO ARRAYS OF DIMENSION OF (NSD+NSW) 
AND (NSD+NSW) x (NSD+NSW). 


THIS COMPUTER PROGRAM REQUIRES SEVEN INPUTS WHICH ARE ENTERED # 
FROM AN INPUT FILE NAMED "MON. DAT" WITH FORMAT 17. +t 
* 

b = AL = MONOPOLE RADIUS IN WAVELENGTHS # 

h = HL = MONOPOLE LENGTH IN WAVELENGTHS * 

€ = 277a/) = EPSLN = NORMALIZED GROUNDPLANE RADIUS IN WAVENUMBERS # 
b,/b = BAR = RATIO OF OUTER TO INNER CONDUCTOR RADII * 

OF THE COAXIAL LINE FEED + 

N = NSW = NUMBER OF MONOPOLE SEGMENTS # 

M = NSD = NUMBER OF GROUNDPLANE CONCENTRIC ANNULAR ZONES * 

X = IX = NUMBER OF SUBDIVIDED ANNULAR ZONES FOR EACH NSD # 

+ 

+ 

+t 

* 

4 

+ 


HAA PWOoONnQAanAaAaraGAHhNonAaAAniaadthrtaGaanriAQaninganrangaaqanoarioaeonannh 


tHE HE a UG HE tt aE a TE UE tHe HE HE a SE a aE a HE He HE HE HE aE HE HE Se HE a He HE HE EE HE HE HE a aE Ha HE HE HE HE HE aE TE HEH HE HE EE HE HH HE HE EE HEE Ht tt 
IMPLICIT REAL#8 (A-H), (P-Z) 

COMPLEX#16 CU(30), ZJU(30), VIU(30, 30), ZIU(30, 30), J11 

COMPLEX#16 Y11,EC,D11,D12, D21, D22. DZ1J, DV1, W11 

COMPLEX#16 P11,P12,P21, P22, ZDD, ZDM, ZMD. ZMM, Z11. ZD, Z22. Z12, Z21 
DIMENSION FB(S500) 

OPEN(UNIT=4, FILE=’INPUT1 ’, STATUS= ‘NEW ) 

OPEN(UNIT=9, FILE=’RICH’, STATUS® “NEW” ) 

OPEN(UNIT=6, FILE=’MON. DAT’, STATUS ‘OLD’, READONLY ) 

DATA EO, UO0/8. 85418533677E-12, 1. 25663706144E-6/ 

DATA ETA, PI. TP/376. 730366239, 3. 14159265359, 6. 28318530718/ 

DATA ICC. IFB/30, 5S0OO/ 

FORMAT (1X, 215, 2E15. 8) 

FORMAT (/16X, 2F15. 4) 

FORMAT (1X, 215, 4F 10. 4) 

FORMAT (1HO) 

C34 te HE Se aE HE Se Ee HE aE FEE HE SE HE SE EE HE HEE FE ETE HE HE HE HE HE HE EE HE HE HE He HE HE HE HE EE HE HE HH HH EH HH HH a 
C INPUTS # 


cn G) AQ om 
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C 3 EE EH HE HE OE EE EE He He HE HE a EE aE HE Ha a a a 2a a a a a 0 a a a a Oa HE 
READ(6, 17) AL. HL, EPSLN, BAR, NSD. NSW, IX 
17 FORMAT (3F 12. 9,F6. 3, 312) 
101 BL#EPSLN/ (2. #PI1) 
OPEL ELISE SEPT TTT PTT TTT PT TET TTT TTT TT TT TT TTT ET TT TTP TP TTT TPT TT TTT SS 
C INPUT WRITE STATMENTS * 
C #9 Ht 3 0 0 HE SE 9 Ht OE OE He Oe HE HE HE HE HE He HE He HE OH HE HE HE EH aE HE OE Oe aE HE OE EE ae a a a a a a a a a OE a a a a 0 a a 0 a OE a 
103 WRITE (9, 102) 
102 FORMAT (1H1,14X, ‘DOUBLE PRECISION’) 
WRITE (9,80) AL 
80 FORMAT (15X, ‘“MONOPOLE RADIUS(WAVELENGTH), AL= ‘,F15. 8) 
WRITE (9,82) HL 
82 FORMAT (15X, ‘MONOPOLE LENGTH(WAVELENGTH), HL= °.FS. 4) 
WRITE (9,91) EPSLN 
91 FORMAT (15X, ‘DISC RADIUS(WAVELENGTH)., BL® ’,G, ‘/(2"PI)’) 
WRITE (9, 84)NSW 
84 FORMAT (15X,. ‘NUMBER OF MONOPOLE SEGMENTS, NSW ’, I2) 
WRITE (9, 85)NSD 
85 FORMAT (15X, ‘NUMBER OF DISC SEGMENTS, NSD=’, I2) 
WRITE (9, 88)BAR 
88 FORMAT (15X, ’COAXIAL RADIUS RATIO, BAR#® ',F4. 2) 
NLI=NSD+NSW 
WRITE (4, 666)AL, HL, EPSLN, NSD. NSW, NLI, IX 
666 FORMAT(E12. 4, 2X,F8. & 2X, F4. 4, 2X, 4( 12) 2X)) 
81 CONT INVE 
DO 400 IBL#=10, 10 
TL=AL/100. 
IWC J=1 
IZIJ=0 
IWZ=O 
NPH=6 
NEG=NSD+NSW-1 
AK2=TP#AL 
BA=TP#*BL 
HA=TP #HL 
TASTP#TL 
DAD=( BK-AK) /NSD 
DAW=HK /NSW 
RH2=AK+DKD 
IF (RH2. LT. BAR#AK)GO TO 400 
TDKD=2. *DKD 
CDKAD=DCOS ( DKD) 
SDKAD=DSIN(DKD) 
CDK=DCOS (DKW) 
SDK=DSIN( DKW) 
MAX=NSW- 1 
NA#NSD+1 
CALL QMM(AK, DKD, DAW, CDKD, SDKD, CDK, SDK, TK, IWZ, NPH, Z11) 
ZIJ(1,1)22Z11 
IF (NSD. LE. 1)»GO TO 100 
S1=AK 
DO 60 J=2, NSD 
S2=S1+DKD 
S3=S1+TDKD 
T1=AK 
DO 50 I22, J 
T2=T1+DKD 
T3=T1+TDKD 
CALL @DD(CDKD,. SDKD, $1, S83, T1. Ta. TK, IWZ. NPH, Z22) 
ZIVCI,J)#Z22 
390 T1=T1+DKD 
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CALL GDM( AK, DAD, DAW, CDAD, SDKD, SDK. S1, S3, TK, IWZ, NPH, Z12) 
ZIJ(1, J)#Z12 
60 S1=51+DKD 
100 IF(NSW.LE.1)G0 TO 200 
CALL SPART(AK, DAD, DAW, MAX, IWZ, ZU; CU) 
L=0 
DO 160 I=NA, NEG 
DO 150 J=I, NEG 
K=J-I+1 
150 ZIJCT, Jy=ZJ(K) 
L=L+1 
ZIJC1, T)*=CUCL) 
160 CONTINUE 
1768 IFC(NSD. LE. 1)GO TO 200 
Z2=. ODO 
DO 190 J=NA, NEG 
Z2=Z2+DKAW 
S1=Z2-DKW 
S3=Z2+DKW 
RH2=AK 
DO 1680 I=2, NSD 
RH2=RH2+DKD 
T1=RH2-DKD 
T3=RH2+DKD 
CALL SKEWT(AK, S1,S3, T1, T3, CDA, SDK, CDKD. SDKD, IWZ, Z12) 
160 ZIJC(I, J)=Z12 
190 CONTINUE 
200 CALL GRILL(AK, BAR, DAD, DAW, NEG, NSD, NSW. TK, CU) 
DO 210 I=1,NEQ 
DO 210 J=I, NEG 
210 VIJCTI, JySZIJUCT, J) 
OTE CLEC Ce Tee eee eT eT TTT Eee TTT TTT TP Te Tee eS 
C MATRIX INVERSION FOR A NSW+NSD-1 BY NSW*+NSD-1 SQUARE MATRIX # 
OTE EE EEE ETS EEE eT PTET PTET Tree eT TTT Pee Tee Tere ST 
CALL CROUT(ZIJ, CJ, ICC, 0. IWCJ, 1, NEG) 
¥11=CU(1) 
Z11=1. /Y11 
320 CONTINUE 
J11=Z11 
CM FE HE aE OE IE aE bE HE HE OE EE He HE HE HE HE HE SE HEE HE EE SE IEE EE Ea aE HE aE a aE aE Hae Ee aE aE Ee de Ht at 
Cc OUTPUT WRITE STATEMENTS * 
Cat HE TE HE EE IEE SEE de EE HEE bE OE ba de dd HEH Ea dd de OE J OE EE a AH SEE a Eb Ha aE EE a a dt 
400 CONTINUE 
WRITECS, 92) 
WRITE(9, 93) 
92 FORMAT (//15X, ‘INPUT IMPEDANCE Z11 (OHMS) ‘) 
93 FORMAT (24X, ‘RE(Z11)°, 8X, ‘IM(Z11) %) 
WRITE(9, 2)J11 
CLOSE (UNI T=9, STATUS= ‘KEEP “ ) 
CLOSE (UNI T=6, STATUS# ‘KEEP ” ) 
SOO CALL EXIT 
END 
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CH Bt He HE HE ae EE HE HE a HE a SH a a a a a a a 2 a a 2 a a a HH a a a a a a a 0 a a a 2 a a Ha a 
C SUBROUTINE CISI + 
CC 96 0b He aE HE He SE He HE Se SE Ea HH a a a 0H 0 HS I 2 a a a a a a a a a a 0 Ha a a a a a 2 a a a a a a 
SUBROUTINE CISI(CI,CIN. SI, X) 
IMPLICIT REAL#8 (A-H), (P-Z) 
DATA GAM, P2/. 537721566, 1. 37079632/ 
A=DABS( X) 
IF CA. GT. 4. )GO TO 10 
IF(CA. GT..1)¢G0 TO 3 
IF (A. GT. 0. )GO TO 2 
C1=. oDO 
CIN=. ODO 
SI=. ODO 
RETURN 
2 X2zBAHA 
ST=X#((. OS#X2—-1. )#X2/18. +1. ) 
CIN=. 25#X2H#( (X2/45. —-1. )*#X2/24. +1. ) 
GO TO 8 
3 Y=(4. -A)#(4. +A) 
SI=X#(((C(1. 753141E-9#Y+1. S6898BE-7) #¥+1. 374168E-5) #Y+6. 9398B9E-4) 
C#Y+1. 964B8B2E-2)#Y+4. SISSO9E-1) 
CIN= A#tA#((00(1. 366985E-10#Y+1. 5864996E-8) #Y 
C+1. 725752E-6) #Y+1. 185999E-4) #Y+4. 99O092Q0E-3) #Y+1. 3I1S30GE-1) 
8 CI=GAM+DLOG(A)-CIN 
RETURN 
10 SI=DSIN(A) 
Y=DCOS(A) 
Z=4./A 
U=((((000(€4. 048069E-3#Z-2. 2791 43E-2) #Z+5. S15070E-2) #Z-7. 261 642E-2) 
C#Z+4. 987716E-2) #Z-3. 332519E-3)#Z-2. 314617E-2) #Z—-1. 13495BE-5) #Z 
C+6. 25001 1E-2)#Z+2. SB3989E-10 
V=(((00(000(0-5. 108699E-3#Z+2. 8191 79E-2) #Z-6. S3B72B3E-2) #Z 
C+7. 902034E-2) #Z-4. 40041 6E-2) #2-7. 945556E-3) #Z+2. G01 293E-2) #Z 
C-3. 764000E-4)#Z-3. 122418E-2)#2Z-6. 646441E-7) #Z+2. SOOOOOE-1 
CI=Z#(SIT#V—-V#U ) 
SI=-Z#(ST#U+yveV) +P2 
IF(X. LT. .ODO)ST=-SI 
CIN=GAM+DLOG(A)-CI 
RETURN 
END 
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TESS RSE SEE ETE ET ET TT TTT TT Te PT TP PT PP PP Ph he 
C SUBROUTINE CROUT % 
C2 HE HE TE aE HE TE TE HE He HE HE HE I ae He He HE HE HE HE Oe HE He HE TE HE HE HE aE SE a IE HE HE He HE He aE Te HE EE ae HE He TE TE EE 
SUBROUTINE CROUT(C,S, ICC, ISYM, IWR. 112, N) 
IMPLICIT REAL#8 (A-H), (P=-Z) 
COMPLEX#16 C(ICC, ICC), S(1) 
COMPLEX#16 F,P,SS.T 
FORMAT(15X, ‘LIST OF THE LOOP CURRENTS CUJ(I)’%) 
FORMAT(5X, ‘1’, 3X. ‘MAGNITUDE’, 3X. ‘PHASE’, SX, ‘REAL ’ 
2, 5X, “IMAGINARY ‘> 
FORMAT(1X,15,F10. 4, F10. 1, 2E12. 4) 
FORMAT (1HO) 
IF(I12.NE.1)¢0O TO 22 
IF (CN. E@. 19$(1)=S(1)/C(1, 1) 
IF(N.EQ@.1)GO TO 39 
IFCISYM. NE. ODO)GO TO 8 
DO 6 I=1,N 
DO & J=I,N 
& C(J, Td=CC(I, J) 
8 F=HC(1, 1) 
DO 10 L=2,N 
10 C(1,Ld=C(1,L)/F 
DO 20 L=2,N 
LLL=L~1 
DO 20 I=L,N 
F=C(I,L) 
DO 11 Kel,tLbLbL 
11 F=F-C¢(I,K)#C(K,L) 


Wi G) Me 


C(I,LI=F 

IF(L.EQ@. 1)GO TO 20 
P=C(L,L) 
IFCISYM. EQ. ODO)GO TO 15 
F=C(L, I) 


DO 12 K=1,LLL 
12 FHF-C(L,K)#C(K, I) 
C(t. I) =F /P 
GO TO 20 
15 F=CC(I,L) 
C(t, Il) =F/P 
20 CONTINUE 
22 DO 30 L#=1,N 
P=C(L,t) 
T=S(L) 
IF(L.EQ@. 1)GO TO 30 
LLL=t-1 
DO 25 A=1,tbb 
25 T=T-C(L, K)#SCK) 
30 S(L)=T/P 
DO 38 L=2,N 
I=N-L+1 
II=I+1 
T=S(1) 
DO 35 K=II.N 
35 T=T-C¢(I,K)#SCK) 
38 SC(I)=T 
39 IFCIWR. LE. ODO) GO TO 100 
WRITE(9, 5S) 
WRITEC(9, 1) 
WRITES, 2) 
CNOR=. ODO 
DO 40 I=i,N 
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SA=CDABS(S(T)) 

40 IF(SA. GT. CNOR)CNOR#SA 
IF (CNOR. LE. 0. )CNOR#=1. 
DO 44 I[=1,.N 
SS#S(1) 
SA=C DABS (SS) 
SNOR=SA/CNOR 
PH=. ODO 
IF (SA. GT. 0. )PH#®57. 29578#DATAN2 (DIMAG(SS),. DREAL(SS) ) 
WRITE (9, 3)1, SNOR, PH. SS 

44 WRITE(4, 111)SS 

111 FORMAT (2E12. 4) 
WRITEC(9, 5) 

100 RETURN 
END 
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OTE LESS S PETE PTT TT TT PP PTT TTT TPT TTT TT TTT So 


C 


SUBROUTINE EXPJ 


# 


C Ht HE EEE EE EE EEE EEE EEE EEE EEE EERE ERE REE EE EEE ERREREEEE EEE S 


10 


15 


20 


30 


et | 


32 


SUBROUTINE EXPJU(V1, V2. W12) 
IMPLICIT REAL#8 (A-H), (P-Z) 


COMPLEX#16 EC,.E15.S, T, UC, VC. V1. V2, W122, Z 
DIMENSION V(21),W(21),D(16),E(16) 


DATA V/ 0. 22284667E OO, 


20. 11889321E 01,0. 29927363E 01,0. 57751436E 01,0. 
20. 15982874E 02,0. 93307812E-01, 0. 49269174E OO, O. 
20. 22699495E 01,0. 36676227E 01,0. 54253366E 01,0. 
20. 10120228E 02,0. 13130282E 02.0. 16654408E 02,0. 
20. 25623894E 02,0. 31407519E 02,0. 38530683E O2, O. 


DATA W/ 0. 45896460E OO, 


20. 41700083E 00, 0. 11337338E 00,0. 10399197E-01, 0. 
20. 89854791E-06. 0. 21823487E 00.0. 34221017E 00.0. 
20. 12642582E 00. 0. 40206865E-01, 0. 85638778E-O2, O. 
20. 11167440E-03, 0. 64599267E-05, 0. 22263159E-064, O. 
20. 39218973E-10, 0. 14565152E-12, 0. 14830270E-15,. 0. 


DATA D/ 0. 22495842E 02, 
2 0. 74411568E O02, -0. 41431576E 
2 0.16021761E 03, -0. 23862195E 
2 0. 12254778E O02, -0. 10161976E 
2-0. 21069574E O02, 0. 22046490E 
DATA E/ 0.21103107E O2., 
2-0. 37959787E 03. -0. 97489220E 
2-0. 12910931E O03, -0. 55705574E 
2 0.17949528E O02. -0. 32981014E 
2 0. 22236961E O02. 0. 39124892E 
zZ=V1 

DO 100 JIM=1.,2 

X=DREAL(Z) 

Y=DIMAG(Z) 

EiS=(. ODO, . ODO) 

AB=CDABS(Z) 

IF (AB. EQ. 0. )GO TO 90 


03, ~-0 
03, -0 
O02, -O 
01, O 
O02, O 
O03, O 
00, O 
O02, O 


. 78754339E 
. 90094687E 
4721959 1E 
. 89728244E 


. 12900672E 
.135924801E 
. 31028836E 
. 81636799E 


IF (X.GE.O. . AND. AB. GT. 10. )GQ TQ 80 


YA=DABS(Y ) 


IF(X.LE.O. . AND. YA. GT. 10.)GO TO 80 
IF (YA-X. GE. 17. 5. OR. YA. GE. &. 5. OR. X+VA. GE. 5. 5. OR. X. GE. 3. 2GO TO 20 


IF (X.LE. -9. )GO TQ 40 

IF (YA-X. GE. 2.5)G0O TO SO 
IF (X+YA. GE. 1.5)9¢0 TO 30 
N=6. +3. #AB 

E1S=1./(N-1. )-Z/N#¥#2 
N=N-1 

E1581. /(N-1. )-Z#E15/N 
IF(N. GE.3)GO TO 15 


O02, O. 
03, -O. 
O1, O 
O1/ 

03, O 
O02, O 
O02, O 
O1/ 


E15=Z#E15-DCMPLX(. 577216+DLOG (AB), DATAN2(Y, X)) 


GO TO 90 

J1=1 

J2=6 

GO TO 31 

J1=7 

J2=21 
S=(. ODO, . ODO) 
YS=Y#Y 

DO 32 I=J1, J2 
XI=VCI) +X 
CFEWCI)/(XI*XI+YS) 
S=S+DCMPLX (XI#CF., -YA#CF ) 


150 


9B374674E O1, 
12155954E O1, 
7T56S59162E O1, 
20776479E Oe. 
48026086E O02/ 


26101720E-03, 
26302758E OO, 
12124361E-02, 
42274304E-O8. 
16005949E~-19/ 


11254744E 
68487854E 
79729681E 


. 17949226E 
. 14696721E 
.81657657E 


02, 
02, 
O1, 


02, 
03, 
O1, 
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40 


50 


32 


34 


36 


80 


90 


100 


GO TO 54 

TSEX#X—-YHY 

T4=2. ¥X#YA 

TS*X*T3-YA#T4S 

T&=X#T4+YAHTS 

UC=DCMPLX (D(11)4D(12)4#X+D(13) 4T3+TS-E(12) #YA-E(13)4T4, 


2 E€(11)4E(12)#X+E(13) #T3+T64+D(12) #YA+D (13) #T4) 


VC=DCMPLX (D(14)4+D(15)#X+D(16) #T3+T5-E(15) #YA-E (16) 4#T4, 


2 E€(14)4E(15) #X+E (16) #T3+T64+D(15) #YAt+D(16) #74) 


GO TO S52 

T3BX#X-YRY 

T4222. #X#YA 

TS=X#TS-YA#T4S 

T&=X#T4+YAHTS 

T7=X#TS-YA#TS 

TB=X*#TS+YA#TS 

T9I=X*#T7—-YA#TS 

T1O=X#*TE+YA#T7 

UC=DCMPLX (D(1)4+D(2)4#X4+D(3) #T34+D(4) #TS54+D( 5) *#T7+T9I-(E (2) #YAtE (3) 4T4 


2tE C4 #TSO+E (SI#TS), EC LIFE (CQ) #X+E (SD #TSAVE (4) # TSE (5) #T7+T10+ 
3(D(2)*#YA+D (3) #T44+D(4)#7T64+D(5)#T8) ) 


VC =DCMPLX (D(6)4+D(7) #X+D(8)#T34+D(9) #T54+D(10) #7 7+T9-(E (7) #YA+E (8) *T4 


24+E(9)*#T64E(10)#T8),E(6) +E (7) ¥X4E (8B) #TS+E (9) #TS+E (10) #T74+T10+ 
3(D(7) #YA+D(8)#T44+D(9)#T64+D(10) #T8) ) 


EC=s=UC/VC 

S=EC/DCMPLX(X, YA) 

EX=DEXP(-X) 

T=EX*DCMPLX (DCOS( YA), -DSINCYA)) 

E15#S+*T 

IF CY. LT. O. 1E15=#DCONJG (E15) 

GO TO 90 

E15=. 409319/(Z+. 193044)+. 421831/(Z+1. 02666)+. 147126/(Z+2. 56788) + 


2. 206335E-1/(2Z+4. 90035)+. 107401E-2/(Z+8. 18215)+. 15B8654E-4/ (Z+ 
312. 7342)+. 317031E-7/(24+19. 3957) 


E1S5=E15*CDEXP (-Z) 

IF (JIM. E@. 1)W12"E15 

Z=V2 

Z#V2/V1 

TH=DATAN2 (DIMAG(Z), DREAL(Z) )-DATAN2(DIMAG (V2). DREAL (V2) ) 


2+DATAN2(DIMAG(V1), DREAL (V1) ) 


AB=DABS( TH) 

IF(AB LT. 1. ) TH=. ODO 

IF (TH. GT. 1. )TH#®6. 2831653 

IF (TH. LT. -1. ) TH=-6. 28631853 
W12=W12-E15+DCMPLX(. ODO, TH) 
RETURN 

END 
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CHEER HH 


C 


SUBROUTINE GRILL * 


OLAS AAALALAL LSS L ETE ST SPT ET TPP TT TTT TT TT eT CPP eT TP ee eee ee SS 


20 


40 


SUBROUTINE GRILL(AK, BAR. DAD, DKW, NEG, NSD. NSW, TK, VU) 
IMPLICIT REAL#8 (A-H), (P-Z) 
COMPLEX#16 EGZ, GM, GP, GI(20),VU(1), GII, @ST. WST, V1 
DOUBLE PRECISION DZS,RS1.RS2 
DATA PI, TP/3. 14159265359, 6. 28318530718/ 
IDM=20 

DO 20 I=1, NEG 
VICI) =¢(. ODO, . ODO) 
VI(1)=C1.,. 0D0) 

IF (BAR. LE. 1. )RETURN 
VJ(1)=(. ODO, . ODO) 

DK=DKW 

SDK=DSIN(DK) 

CDK=DCOS( DK) 

BAL=DLOG( BAR) 
QST=DCMPLX(. ODO, 1. /(4. *BAL*SDK) > 
BA=AK#BAR 

AKS=AK#AK 

BAS=BK#BK 

LIM=NSW+1 

IF (LIM. GT. IDM)LIM=IDM 

NPH=6 

NPH=2+ (NPH/2) 

NPP=NPH+1 

PHA=. 0174533#20. 

DPH=PHA/NPH 

PH=. ODO 

DO 70 LPH=1,2 
WST=DPH#Q@ST/(3. #PI) 

SGN=-1. 

DO 80 IPH=1, NPP 

WF=3. +SGN 

IF CIPH. EG. 1)WF=1. 

IF (IPH. EQ. NPP )WF=1. 
CPH=DCOS(PH) 

IF (IPH. GT. 1)G0 TO 40 

IF (LPH. GT. 1)9GO TO 40 
CPH=DCOS(DPH/10. ) 

RS1=2. *AKS#(1. -CPH) 
RS2=AKS+BKS-2. *AK*#BK#CPH 
RH1=DSQGRT(RS1 ) 
RH2=DSGRT(RS2) 

CALL CISI(CA, CIN, SA, RH1 >) 
CALL CISI(CB,. CIN, SB, RH2) 
GI(1)=2. *DCMPLX (CB-CA, SA-SB) 
DO 50 I=2,LIM 

DZ=DK#(I-1) 

DZS=DZ*DZ 

RA=DSGRT(RS1+DZS) 
RB=DSQRT(RS2+DZS) 

CALL CISI(C1,CIN, S1, RA+DZ) 
CALL CISI(C2, CIN, S2, RB+DZ) 
GP=DCMPLX (C2-C1,.S1-S2) 
RAM=RS1/(RA+DZ) 

RBM=RS2/ (RB+DZ) 

CALL CISI(C1.,CIN, S1, RAM) 
CALL CISI(C2, CIN, S2, RBM) 
GM=DCMPLX (C2-C1,51-S2) 
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EGZ=DCMPLX (DCOS(DZ), DSIN(DZ)) 
30 GI(I)=GP#ECZ+OGM/EGZ 
VJ(1)=VJ(1)+WF#WST#(GI(2)-CDK#GI(1)) 
IF (NSW. LE.1)¢60 TO 78 
K1=0 
IA=NSD+1 
DO 60 I=IA, NEG 
K1=K1+1 
K2=K1+1 
K3=K2+1 
IF(K3. GT. IDM)GO TO 60 
GP=GI(K1)—-2. #CDK#GI (K2)+GI (KS) 
VICTI=VUCT) +WE#WST#GP 
60 CONTINUE 
78 SGN=-SGN 
80 PH=PH+DPH 
DPH=(PI-PHA) /NPH 
90 PH=PHA 
CALL CISI(CA, CIN, SA, AK) 
CALL CISI(CB, CIN, SB, BK) 
R2=AK+DKD 
SR2=DSIN(R2) 
CR2=DCOS(R2) 
SDKD=DSIN(DKD) 
V11=(SR2*(CB-CA)-CR2#(SB-SA) )/ (2. *BAL#SDKD) 
VJC1)2V114V5(1) 
IF (NSD. LE. 1)RETURN 
V22=(DSIN( AK) #(CB-CA) —-DCOS( AK) #(SB-SA) )/(2. *BAL*SDKD) 
VJ(2)=DCMPLX (V22, . ODO) 
RETURN 
END 
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SUBROUTINE QDD(CDKD. SDKD. $1, S3, T1, TS, TK, IWZ, NPH, Z22) 
IMPLICIT REAL#8 (A-H). (P~-Z) 

COMPLEX#16 Z12, Z22, ZD 

DATA PI,P3/3. 14159265359, 9. 42477796077/ 
FORMAT (1X, 8F10. 2) 

FORMAT (1HO) 

FORMAT(SX, ‘PRINTOUT FROM QDD‘) 
FORMAT(S5X, ‘DISK DIPOLE TO DISK DIPOLE’) 
IF(IWZ.LE.O0)GO TO 10 

WRITE(8&, 6) 

WRITE(&, 7) 

WRITE(8, 5) 

PHA=. 0174533*2. 

PHB=. 0174533#20. 

NPH=2% (NPH/2) 

DP H=PHA/NPH 

PH=. ODO 

NPHP=NPH+ 1 

Z22=(. ODA, . ODO) 

DO 80 IPH=1,3 

ZD=(. ODO, . ODO) 

SGI=-1. 

DO 70 I=1,NPHP 

WF=3. +SGI 

IF CT. EQ. 1)WF=1. 

IF CT. EQ. NPHP)WF=1. 

CPH=DCOS (PH) 

IF (I. EQ@.1 . AND. IPH. GT. 1)9GO TO 60 


NEXT: DISK DIPOLE TO DISK DIPOLE. 


60 


70 


80 


CALL SKEWS(S1,S3, T1, T3, TK, CDKD. SDKD,. CDKD, SDKD. CPH, Z12) 
PHD=57. 29578#PH 

IF CIWZ. GT. ODOIWRITE(8, 2)PHD, Z12 
SCI=-SGI 

PH=PH?+DPH 

ZD=ZD*+WF#Z12 

Z22=Z22+DPH*ZD/P3 

PH=PHA 

DP H= (PHB-PHA) /NPH 
IFCIPH EQ 1)6G0 TO 80 

PH=PHB 

DPH=(PI-PHB) /NPH 

CONT INUVE 
IFCIWZ. GT. ODOIWRITE(CSB. 5) 

RETURN 

END 
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C46 36 SE FE FE aE aE HE HE aE TE EE IE HE SE SE SE TE JE EE HE HE JE SE EE JE SE SE SE SE HS EE JE SE JS SEE JE SS SE J SE SE IE SH EE J HE HE EE gE EE HE 
C SUBROUTINE GDM S 
Ce 9 He HE J He aE Fe TE EH SE IE SE SE SE SE IE SE ESE HE SE SE EE ESE SE HS EE SE ES OE HES HE EJ ES EE SE SE EE dH Ed Ht 
SUBROUTINE QDM(AK. DKD., DKW, CDKD, SDKD, SDK, S11, S3, TK, IWZ, NPH, Z12) 
IMPLICIT REAL#8 (A-H), (P-Z) 
COMPLEX#16 2D, Z12, Z21, ZDM, ZM. P12, ZDD. PDM 
DATA PI,P3/3. 14159265359, 9. 42477796077/ 


2 FORMAT (1X, 8F 10. 2) 
-) FORMAT (1HO) 
& FORMAT( 5X, ‘PRINTOUT FROM Q@DM’) 
7 FORMAT( 5X, ‘DISK DIPOLE TO DISK MONOPOLE ’ ) 
8 FORMAT(SX, ‘DISK DIPOLE TO WIRE MONOPOLE ’) 
IF(IWZ.LE.0)GO TO 10 
WRITE(S, 6) 
WRITE(8, 7) 
WRITE (8, 5) 


10 PHA=. 0174533*2. 
PHB=. 0174533#20. 
NP H=2% (NPH/2 ) 
NPHP =NPH+ 1 
DPH=PHA/NPH 
IDM=1 
IF(S1. GT. 10. #AK) IDM=0O 
PH= ODO 
ZDD=(. ODO. . ODO) 
T2=AK+DKD 
DO 40 IPH#1,3 
ZD=(. ODO. . ODO) 
SGI=-1. 
DO 30 11. NPHP 
WF=3. +SGI 
IFCI. EG. 1)WF=1. 
IF (I. EG. NPHP)WF=1. 
CPH=DCOS (PH) 
IFCI.EG@.1 . AND. IPH. GT. 1)GO TO 20 
C NEXT: DISK-DIPOLE TO DISK-MONOPOLE. 
CALL ZSDM(S1,S3. AK. T2, TK, CDKD, SDKD, SDKD, CPH, P12) 
20 PHD#=57. 29578#PH 
IF CIWZ. GT. OOWRITE(8, 2) PHD, Pil2 
SGT=-SGI 
PH=PH+DPH 
30 ZD=ZD+WF#P12 
ZDD=ZDD+DPH*ZD/P3 
PH=PHA 
DPH=(PHB-PHA) /NPH 
IF (IPH. EG@.1)GO TO 40 
PH=PHB 
DPH=(PI-PHB)/NPH 
40 CONTINUE 
IFC IWZ. GT. O)WRITE(CS, 5) 
W2=TK+DKW 
IF (IDM. EQ@.0)GO TO 100 
IFCIWZ. GT. O)WRITE(S., 8) 
LPH=6 
LPH=2#(LPH/2) 
LPP=LPH+1 
PHA=. 0174533#20. 
DPH=PHA/LPH 
PH=. ODO 
ZDM=(. ODO, . ODO) 
RMN=AK/ 100. 
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DO 90 IPH#i,2 
PDM=(. ODO. . ODO) 
SGI=-1. 
DO 80 11, LPP 
WF=3. +SGI 
IF C(I. EQ. 1)WF=1. 
IF (I. EQ. LPP )WF#1. 
CPH=DCOS(PH) 
RH=AK#DSIN(PH) 
IF (CT. EQ. 1) RH=RMN 
AC =AK#CPH 
V1i=S1-AC 
V3=S3-AC 
IFCI.EQ@.1 . AND. IPH. GT.1)9GO TO 70 
C NEXT: DISK-DIPOLE TO WIRE-MONOPOLE. 
CALL ZSDM(V1, V3, TK, W2, RH, CDAD, SDKD, SDK. . ODO, Pi2) 
70 PHD=57. 293578*#PH 
IFCIWZ. GT. ODO)WRITE(S8, 2) PHD, Pi2 
SGI=-SGI 
PH=PH+DPH 
80 PDM=PDM+WF#Pi2 
ZDM=ZDM+DPH*#PDM/P3 
DPH=(PI-PHA) /LPH 
90 PH=PHA 
Z1i2=ZDM-ZDD 
IFCIWZ. GT. ODO)WRITE(S, 5) 
RETURN 
C NEXT: DISK-DIPOLE TO WIRE-MONOPOLE. 
100 CALL ZSDM(S1,S3, TK, W2,. AK, CDKD, SDAD,. SDK. . ODO. ZDM) 
Z12=ZDM-ZDD 
RETURN 
END 
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SUBROUTINE QMM(AK, DAD, DAW, CDAD, SDKD, CDK, SDK, TK, IWZ. NPH, Z11) 
IMPLICIT REAL#8& (A-H), (P-Z) 
COMPLEX*16 FDM, FMD, FMM, PDM, PMD, PMM 
COMPLEX#16 ZD, ZDD, ZDM, ZMD, ZMM, Z11,P11 
DOUBLE PRECISION DRG 

DATA PI,P3/3. 14159265359, 9. 42477796077/ 
FORMAT (1X, SF8. 0) 

FORMAT (1X, 8F10. 2) 

FORMAT (1HO) 

FORMAT (5X, ‘PRINTOUT FROM QMM’) 
FORMAT( SX, ‘DISK MONOPOLE TO DISK MONOPOLE ”) 
IF(IWZ.LE.0)GO TO 10 

WRITES, 6) 

WRITE(S, 35) 

AKS=AK#AK 

DAMP =TK+DKAW 

CDAMP=DCOS(DKMP ) 

SDAMP=DSIN (DAMP ) 

ZMM=(. ODO, . ODO?) 

ZMD=(. ODO, . ODO? 

ZDM=(. ODO, . ODO) 

LPH=6 

LPH=2%(LPH/2) 

LPP=LPH+1 

PHA=. 017453320. 

DPH=PHA/LPH 

RMN=AK#DPH/ 10. 

PH=. ODO 

DO 44 IPH=1,2 

FDM=(. ODO, . ODO?) 

FMD=(. ODO, . ODO? 

FMM=(. ODO, . ODO) 

SGI=-1. 

DO 40 I=1, LPP 

WF=3. +SGI 

IF (I. EQ. 1)WF=1. 

IF (I. EQ. LPP )WF=1. 

IF(I.€Q@.1 . AND. IPH. GT. 1)GO TO 38 
CPH=DCOS(PH) 

SPH=DSIN(PH) 

DRG=2. *(1. -CPH) 

R=AK*DSGRT (DRG) 

IF (I. EQ. 1)R=R+RMN 


C NEXT: WIRE MONOPOLE TO WIRE MONOPOLE. 


CALL ZSMM(-TK, DKW, . ODO, DAW, R, CDAMP, SDAMP, SDK, 1. .PMM) 
R=AK#SPH 

IF (I. EQ@. 1)R=R+RMN 

TI=AK#(1. -—CPH) 

T2=T1i+DKD 


C NEXT: WIRE MONOPOLE TO DISK MONOPOLE. 


CALL ZSMM(-TK, DAW, T1, T2,R, CDKMP, SDKMP, SDKD, . ODO, PMD) 
S1=AK*(1. -CPH) 
S2=S1+DKD 


C NEXT: DISK MONOPOLE TO WIRE MONOPOLE. 


38 


CALL ZSMM(S1, S2. TK, TK+DKW, R, CDAD, SDKD; SDK, . ODO, PDM) 
FMM=FMM+WF #PMM 
FMD=FMD+WF #PMD 
FDM=F DM+WF #P DM 
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PHD=357. 29578*#PH 
IF CIWZ. GT. OOWRITE(&. 1) PHD, PMM, PMD, POM 
SGI=~-SGI 
40 PH*PH+DPH 
ZMM=ZMM+DPH#FMM/PS 
ZMD#ZMD+DPH#FMD/PS 
ZDM*ZDM+DPH#FDM/P3 
RMN=. ODO 
DPH=(PI-PHA)/LPH 
44 PH=PHA 
IF CIWZ. GT. O)WRITE(S, 5) 
IF CIWZ. GT. OOWRITE(S, 7) 
PHA=. 017453932. 
PHB=. 017453320. 
NPH=2* (NPH/2) 
NPHP =NPH+1 
DPH*PHA/NPH 
PH=. ODO 
ZDD=(. ODO. . ODO) 
S22=AK+DKD 
DO 60 IPH=1,3 
ZD=(. ODO. . ODO) 
SGI=~1. 
DO 50 I=1,NPHP 
WF=3. +SGI 
IF CI. EQ. 1)WF=1. 
IF (CT. EQ. NPHP )WF #1. 
CPH=DCOS(PH) 
IF CI. EQ. 1 . AND. IPH. GT. 1)GO TO 48 
C NEXT: DISK MONOPOLE TO DISK MONOPOLE. 
CALL ZSMMCAK, S2, AK, S2,. TK, CDKD, SDKD. SDKD. CPH, P11) 
48 PHD=57. 29578#PH 
IF CIWZ. GT. OOWRITE(&, 2)PHD, Pli 
SGI=-SGI 
PH#PH+DPH 
3O ZD=ZD+WF#P1ii 
ZDD=ZDD+DPH#*ZD/P3 
PH=PHA 
DPH=(PHB~PHA) /NPH 
IFCIPH.EQ@.1)GO0 TO 60 
PH=PHB 
DPH=(P I-PHB ) /NPH 
60 CONTINUE 
Z11=ZDD-ZDM~ZMD+ZMM 
IF CIWZ. GT. O)WRITE(CS. 5) 
RETURN 
END 
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30 


SUBROUTINE SKEWS(S1,S3, Ti, T3, RHK, CDK, SDK, CDKT, SDKT, CPSI, Z12) 
IMPLICIT REAL#& (A-H), (P-Z) 
COMPLEX#14 Z12, EIN, EGDZ, CQX, EUXX, EM, EP 
DIMENSION S(3),T(3) 

DOUBLE PRECISION R.SI,SIS, TJ, TUS, XXD.S,T,.CI, CU, FK, FL, DZ 
DOUBLE PRECISION RHS, DZS, DDZ, RIV, DPSI, CPSS, SPSI, RH2 
DATA ETA. PI, TP/376. 730366239, 3. 14159265359, 6. 28318530718/ 
S§(1)=51 

S(2)=(61+53)/2. 

S$(3)=S3 

T(1)#T1 

T(2)2=(T14TS) /e2. 

T(Qd)=T3 

Z1i2=(. ODO, . ODO) 

DPSI=CPSI 

IF (DABS(CPSI).LT. . 999999)G0 TO 10 
RHS=RHKA#RHK 

CPSS=DPSI#*DPSI 

IF (CPSS. GT. 1. DO)CPSS=1. DO 
SPSI=DSGQGRT(1. DO-CPSS) 
RH22=SPSI#(T1+T3) /2. 

RHS#RHS+RH2#*RH2 

SGN=1. 

IF(CPSI. GT. .0D0)GO TO 60 

SGN=-1. 

T(1)=-T3 

T(2)8-T(2) 

T(3)=-Tl 

GO TO &0 

DO 350 11,3 

ST=S(1I) 

SIS=SI#SI 

Ci=1. 

IF(I. EQ. 2)CI=—-2. #CDK 

TJU=T(U) 

TUS#=TU#TJ 

R=DSQRT(SIS+TJS—-2. *SI*#TJ#DPSI ) 

CU=1. 

IF (J. EQ. 2)CJ=-2. *CDAKAT 

CQ@x=(. ODO, . ODO) 

DO 40 Ke1,2 

FA=(-1. )##K 

DO 40 L=1.2 

FL=(-1. )##L 

XXD=FA#SI+FLe#TU 

XX=XXD 

EUXX=DCMPLX (DCOS (XX), DSIN( XX) ) 
XXX#R+XXD 

X=DABS (XXX) 

CALL CISI(COSI,CIN, SINI, X) 

IF (XXX. LT. . ODOISINI=-SINI 
CQX=CQX+DCMPLX (COSI, -SINI ) #EUXX#FA#FL 
CONT INUE 

Z12=Z12+CQX#CI#CJU 

CONTINUE 

Z12=-ETA#Z12/(8. #PI¥#SDK*#SDKT ) 

RETURN 
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60 


70 
80 


DO 80 [1,3 

SI=S(I) 

CI=1. 

IF (CI. EQ. 2)CI=-2. *#CDA 
CQ@x=(. ODO, . ODO?) 

DO 70 J=1,3 

TJUH=T (JU) 

CJ=1. 

IF (J. EG. 2)CJ=-2. *#CDAT 
DZ=TJ-SI 

R=DSGRT (RHS+DZ#DZ ) 

X=R+DZ 

IF (DZ. LT. . ODO) X*®RHS/ (R-DZ) 
CALL CISI(COSI, CIN, SINI, X) 
EP=DCMPLX(COSI, -SINI) 

X=R-DZ 

IF (DZ. GT. . ODO) X=RHS/ (R+DZ) 
CALL CISI(COSI, CIN. SINI, X) 
EM=DCMPLX (COSI, -SINI) 
EGDZ=DCMPLX (DCOS( DZ), DSIN(DZ) ) 
CQX=CGX+CJ# (EP #EGDZ+EM/EGDZ) 
Z12=Z12+C I1*CQXx 

Zi 2=SGCN#ETA#Z12/(8. #P 1 #SDK#SDKT > 
RETURN 

END 
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SUBROUTINE SKEWT(AK, $1,553, T1. T3. CDK, SDK, CDAD, SDAD, IWZ, Z12) 
IMPLICIT REAL#G (A-H), (P-Z) 

COMPLEX#16 P12.,Z12,Q12 

DATA PI,P3/3. 14159265359. 9. 42477796077/ 
FORMAT (1X, 8F10. 2) 

FORMAT (1HO) 

FORMAT(5SX,. ‘PRINTOUT FROM SKEWT ’) 
FORMAT (SX. ‘WIRE DIPOLE TO DISK DIPOLE’) 
IF(IWZ.LE.0)GO TO 20 

WRITE(G, 6) 

WRITE(@, 7) 

WRITE(G, 5) 

RMN=AK/100. 

NPH=6 

NPH=2# (NPH/2) 

NPP =NPH+ 1 

PHA=. 0174533#20. 

DPH=PHA/NPH 

Z1i2=(. ODO, . ODO) 

PH=. ODO 


WIRE DIPOLE TO DISK DIPOLE. 


50 


60 


80 


DO 80 IPH=1,2 

Q12=(. ODO, . ODO) 

SGI=~-1. 

DO 60 I=1., NPP 

WF=3. +SGI 

IF(I. EQ. 1)WF=1. 

IF (I. EQ. NPP )WF=1 

RH=AK#DSIN(PH) 

IF (IT. EQ. 1) RH*RH+RMN 

AC=AK#DCOS (PH) 

V1=T1-AC 

V3=T3-AC 

IF (I. EQ@.1 . AND. IPH. EQ@.2)¢G0 TO 50 
CALL SKEW (S1,S3, V1. V3. RH, CDK, SDK, CDAD, SDKD, . ODO. P12) 
Q12=Q12+WF#P12 

PHD=57. 29578#PH 

IFC IWZ. GT. ODO)WRITE(&., 2)PHD. Pl2 
SGI=-SGI 

PH=PH+DPH 

Z12=Z12+DPH#*G12/P3 

RMN=. ODO 

DPH=(PI-PHA) /NPH 

PH=PHA 

IFCIWZ. GT. ODO)WRITE(S, 5) 

RETURN 

END 
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SUBROUTINE SPART (AK, DAD, DAW. MAX, IWZ, Z, Z1) 

IMPLICIT REAL#8 (A-H), (P-Z) 

COMPLEX#16 EID( 20), EM(20), EP(20),2Z(1),Z1(1) 
COMPLEX*14 CEM, CEP. EMD, EPD, EMD2,. EPD2, Z11, Z22, G11, @11 
DIMENSION CID(20), SID(20), CM( 20), CP (20), SM( 20), SP(20) 
DOUBLE PRECISION DRG 

DATA GAM, P2/. 577215664, 1. 5370794632/ 

DATA ETA, PI/376. 727, 3. 141597 

IDM=20 

FORMAT(3X. ‘MUST INCREASE DIMENSIONS IN SUBROUTINE SPART’) 
FORMAT (3X, ‘ACTUAL DIMENSION IDM = ’, 15, 6X, 


2’‘REGUIRED DIMENSION MAX2 = ‘’, 135) 


FORMAT (1X, 8F10. 2) 
FORMAT (1HO) 
FORMAT (SX, ‘PRINTOUT FROM SPART’ ) 
FORMAT (5X, ‘FIRST: WIRE DIPOLE TO MODE ONE’) 
FORMAT (SX. ‘ THEN: WIRE DIPOLE TO WIRE DIPOLE’) 
IF (MAX. LE. OVRETURN 
MAX2=MAX+2 

DO 14 I=1,MAX 
Z1(¢1)=¢(. ODO, . ODO) 
Z(I)=(. ODO, . ODO) 

IF (MAX2. LE. IDM)GO TO 16 
WRITE(9, 1) 
WRITE(9, 2) IDM, MAX2 
RETURN 

DK=DKW 

IF(IWZ. LE.0)GO TO 18 
WRITE(86, 6) 

WRITECSG, 7) 

WRITE (6, 8) 

WRITE(SG, 5) 

TDK=2. #DK 
SDKD=DSIN(DKD) 

S11=. ODO 

S13=TDK 

S21=DK 

S23=3. *DK 

DO 20 N=1, MAX2 

I=N-1 

DZ=I#DK 

CID(N)#DCOS(DZ) 
SID(N)=DSIN(DZ) 
EID(N)=DCMPLX(CID(N), SID(N) ) 
CDK=DCOS (DK) 
SDK=DSIN(DK) 
EPD=DCMPLX(CDK,. SDK) 
EMD=DCMPLX(CDK,. -SDK) 
EPD2=EPD#EPD 
EMD2=EMD*#EMD 

CEM=2. *CDK+EMD 

CEP=2. #+CDK+EPD 
AK2=AK#AK 
CSS=ETA/ (8. #PI1*#SDK*#SDK ) 
NPH=& 

NPH=2* (NPH/2) 

NPP=NPH+1 

PHA= 0174533+#20 
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DPH=PHA/NPH 
PH=. ODO 
DO 100 JPH=1,2 
CST=DPH#ETA/ (24. #PI#PI*#SDK*#SDK) 
C22=DPH/ (3. *PI) 
SGN=-1. 
DO 80 IPH=1, NPP 
CPH=DCOS(PH) 
SPH=DSIN(PH) 
IFCIPH. GT. 1)C0O TO 30 
IF (JPH. GT. 1)G0O TO 30 
PHO=DPH/ 10. 
CPH=DCOS (PHO) 
SPH=DSIN(PHO) 
30 T1=AK#(1. -CPH) 
T2=T1+DKD 
RH=AK*#SPH 
DRG=2. *AK2*# (1. -CPH) 
RA=DSGQGRT (DRG) 
RS=DRG 
WF=3. +SGN 
IF (IPH. EG. 1)WF=1. 
IF CIPH. EG. NPP )WF=1. 
WST=WF #CST 
W22=WF#C22 
DO 40 N=1, MAX2 
I=N-1 
DZ=I1#DK 
DZS=DZ*DZ 
R=DSGRT(RS+DZS) 
ARG=R+DZ 
IFCN. EG 1)ARG=RK 
CALL CISI(CP(N), CIN, SP(N), ARG) 
EP (N) =DCMPLX(CP(N),. -SP(N) ) 
IF(N. GT. 1)G0O TO 38 
CM(1)=CP (1) 
SM¢C1)=SP (1) 
EM(1)=EP (1) 
GO TO 40 
38 ARG=RS/ARG 
CALL CISI(CM(N), CIN, SM(N), ARG) 
EM(N)=DCMPLX (CMON), -SMCN) ) 
40 CONTINUE 
R=4, #(-CM(2)+2. #CP(1)-CP(2)) 
At]. #CID(3) #(4+CM(3)-2. #CM(2)4+2. #CP(1)-2. #CP(2)4+CP (3) ) 
B+2 *SID(3)*(-SM(3)42. #SM(2)-2. #SP(2)+SP (3) ) 
X=4. #(SM(2)—-2. #SP(1)4+SP (2) ) 
C+2. *#CID(3) #(-SM(3) 42. #SM(2)-2. #SP(1)+2. #SP(2)-SP(3)) 
D+2. #SID(3) #(-CM(3) +2. #CM(2)-2. #CP(2)4+CP(3)) 
@11=CSS*#DCMPLX(R, X) 
Z(1)=Z¢1)+WST*#DCMPLX(R, X) 
Z11=-6. *CDA#EP(1)+2. #(EPD+CDK) #EP (2) -EPD*EP (3) + 
22. *(EMD+CDK) #EM( 2) -EMD#EM( 3) 
CALL ZSDM(S11, S13, T1, T2, RH, CDK, SDK, SDKD. . ODO, Z22) 
G11=CSS#*Z11-Z22 
Z1(1)=Z1(1)4+WST#Z11-W22*Z22 
PHD=57. 29578#PH 
IFC IWZ. GT. O)WRITE(8, 3)PHD, G11,Q11 
IF (MAX. EQ. 12960 TO 70 
R=2. #CID(2)#(-CM(3)+3. #CM(2)-4. #CP(1)4+3. #CP(2)-CP (3) ) 
E+2. *SID(2)#(4+SM(3)-2. #SM(2)+2. *SP(2)-SPC(3) ) 
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F+CID( 4) #(4+CM(4)—-2. #CM(3)+CM(2)4+CP(2)-2. #CP(3)+CP(4)) 
G+S1D(4)#(-SM(4)+2. #SM(3)-SM(2)4+SP (2) -2. #SP(3)+SP (4) ) 
X=2. *CID(2)#(SM(9)-3. #SM(2)+4. #SP(1)-3. #SP(2)+SP (3) ) 
H +2. #SID(2)#(CM(9) <2. #CM(2)42. #CP(2)-CP(3)) 
I+CID(4)#(-SM(4)+2. #SM(3) -SM(2)-SP(2)+2. #SP(3)-SP(4)) 
J+#S1ID(4)4#(-CM(4)42. #CM(S) -CM(2)4+CP(2)~-2. #CP (3) +CP (4) ) 
Z(2)=2Z(2)4WST#DCMPLX(R, X) 
Z11=2. #EP(1)+CEM*#( EMD#EM(3) -EPD#EP (2) )+ 
2CEP#(EPD#EP (3) -EMD#EM(2) ) -EPD2#EP (4)-EMD2#EM( 4) 
CALL ZSDM(S21, S239. T1., T2,. RH, CDK, SDK, SDKD, . ODO, Z22) 
Z1(2)=2Z1(62)4WST#Z11-W224#Z22 
IF (MAX. EQ. 2)9G0 TO 70 
Si=DK 
DO 60 N#3, MAX 
Mi=N-1 
M2=N-2 
Ni=N+1 
N2=N+2 
Z11=EP(M1)#EID(M1)+EM(M1) /EID(M1 )+CEP#(EP(N1)#EID(N)-EM(N) /EID(N) ) 
2-CEM*#*(EP(N) #EIDCN)-EM(N1)/EID(N) )-EP(N2)#EID(N1 )-EM(N2) /EID(N1 ) 
S$1=S1+DK 
S3=S1+TDK 
CALL ZSDM(S1, S93. T1, T2, RH, CDK, SDK, SDKD, . ODO, Z22) 
Z1(N)=Z1(N)4WST#Z1 1-W22*Z22 
CPA=CP (M2)-2. #CP (M1 )+CP(N) 
CPB=2. #*CP(N)-CP(M1)-CPC(N1) 
CPC=CP (N2)-2. #CP(N1)+CP(N) 
CMA=CM(M2)-2. #CM(M1)+CM(N) 
CMB=2. #CM(N)-CM(N1)-CM(M1 > 
CMC=CM(N2)—-2. #CM(N1)+CM(N) 
SPA=SP (M2)—-2. #SP(M1)+SP(N) 
SPB=2. #SP(N)-SP(M1)-SP(N1) 
SPC=SP (N2)-2. *SP(N1)+SP(N) 
SMA=SM(M2)-2. #*SM(M1)+SM(N) 
SMB=2. #*SM(N)—-SM(N1)-SM(M1 > 
SMC =SM(N2)-2. #SM(N1)*SMON) 
R=CID(M2) #(CPA+CMA) +2. #CID(N) #(CPB+CMB) +2. #*SID(N) #(SPB-SMB) 
KA *#+CID(N2)#(CPC+CMC )+SID(N2) #(SPC-—SMC ) 
IF (CN. GT. 3)R=R*+SID( M2) *(SPA-SMA) 
X=-CID(M2)#(SPA+SMA)—-2. *CID(N) #(SPB+SMB)+2. #*SID(N)#(CPB-CMB) 
L -CID(N2)#(SPC+SMC )+SID(N2)#(CPC-CMC ) 
IF(N. GT. 3)X=X+SID(M2)*#( CPA-CMA) 
60 Z(N)=Z(N)4+WST#DCMPLX(R, X) 
70 PH=PH+DPH 
80 SGN=-SGN 
DPH=(PI-PHA) /NPH 
100 PH=PHA 
IFC IWZ. GT. OOWRITE(S, 3) 
RETURN 
END 
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SUBROUTINE ZSDM(S1, S3, T1,. T2., RHK,. CDK, SDK, SDAT, CPST, Z12) 
IMPLICIT REAL#8 (A-H), (P-Z) 

COMPLEX#16 E€(2,2),F (2, 2), ES2,ET2, EXPA, EXPB, EGZI,ES1 
COMPLEX*16 CQ@X, EUXX, 712, EP1, EM1,P11,P21,EGDZ., EM. EP 
DIMENSION S(3) 

DOUBLE PRECISION TD1. TD2. SDI.ST1.ST2, R1,R2, SK, ZD, TL1, TL2, SA, SB 
DOUBLE PRECISION R.SI.,SIS, TJ, TUS, XXD, S, T, CI, CU. FK, FL 
DOUBLE PRECISION RHS. DZS, DDZ.RIvJ. DPSI. CPSS. SPSI., RH2 

DATA ETA, PI, TP/376. 730366239, 3. 14159265359, 6. 28318530718/ 
$(1)=9S1 

S$2=(51+S3) /2. 


S(2)=S2 
$(39)=S3 

Z12=(. ODO. . ODO) 
DPSI=CPSI 
CPSS=DPSI#DPSI 


IF (CPSS. GT. 1. DO)CPSS=1. DO 
SPST=DSQGRT(1. DO-CPSS) 

IF (DABS(CPSI). LT. .999999)G0 TO 10 
RHS=RHK#RHK 
RH2=SPSI#(T1+T2) /2. 
RHS=RHS+RH2#RH2 

SGN=1. 

IF (CPSI.GT..OQDO)GO TO 80 
§$(1)=-S3 

S(2)=-S2 

$(3)=-S1 

SGN=-1. 

GO TO 80 

D=RHK 

DSQ=D*D 

CD=D/SPSI 

BD=CD*DPSI 

EB=DEXP(-BD) 

EC=DEXP(-CD) 

TD1=T1 

TD2=T2 

TS1=TD1*TD1 

TS2=TD2*TD2 
CST=-ETA/ (16. #PI1#SDK*#SDKT ) 
SA=S1 

SB=S2 

ET2=DCMPLX (DCOS(T2)., DSIN(T2) ) 
DO 60 ISS=1,2 

IFCISS. €@. 1)ES1=DCMPLX (DCOS(SA),. DSIN(SA) ) 
IFCISS. €@ 2)ES2=DCMPLX (DCOS(SB), DSIN(SB)) 
DO 20 K=1,2 

DO 20 L=1,2 
E(K,L)=(. ODO, . ODO) 

SI=SA 

DO SO I=1,2 

FI=(-1)##]I 

SDI=SI 

SIS=SDI#SDI 

ST1=2. *SDI#TD1*DPSI 

ST2=2. *SDI*TD2#DPSI 
R1=DSGQGRT(DSQ+S1IS+TS1-ST1) 
R2=DSQGRT (DSQ+S1IS+TS2-ST2) 
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EK=EB 

DO 40 K#1,2 

FA=(—-1 ) ##K 

SKA=FK#SDI 

EL=EC 

DO 30 L=1,2 

FL=(-1) ##L 

EKL=EKsEL 

XX=FA#BD+FL#*CD 

TL1=FL*TD1 

TL2=FL*#TD2 

RR1=R1+SK+TL1 

RR2=R2+SK+TL2 

CALL EXPJUCDCMPLX (XX, RR1), DOMPLX (XX, RR2), EXPA) 
CALL EXPJCDCMPLX(-XX, RR1), DCMPLX(-XX, RR2), EXPB) 
E(K, L#E(K,L)+FI* (CEXPA#EKL+EXPB/EKL ) 
EL=1. /EC 

EK=1. /EB 

IF (I. €@. 185)¢0 TO 350 

ZD=SDI#*DPSI 

ZC=ZD 

EGZI=DCMPLX(DCOS(ZC), DSIN(ZC)) 
RR1=R1+ZD-TD1 

RR2=R2+ZD-TD2 

CALL EXPJCDCMPLX(. ODO. RR1), DCMPLX(. ODO, RR2)., EXPB) 
RR1=R1-ZD+TD1 

RR2=R2-ZD+TDe2 

CALL EXPJU(DCMPLX(. ODO. RR1), DCMPLX(. ODO, RR2), EXPA) 
FCI, 1)=#¢. ODO, 2. )#®SDA*EXPA/EGZI 
F(I,2)=(. ODO, 2. )#SDA#EXPB#EGZI 

Sl=S3 

IFCISS. €Q@. 1) 
AP21=CST#((-FC2, 1)-E (2. 2) #ES1+E(1, 2)/ES1)*ET2 
B+( FC2, 2)+E (2, LI #ESI-E(1,1)/ES1)/ET2) 
IFCISS. €Q@ 2) 
CP11=CST#(( F(1, 1)+E C2, 2) *#ES2-E(1, 2) /ES2)#ET2 
D+(-F (1, 2)-E(2, 1 *#ES2+E (1, 1)/ES2)/ET2) 
SA=S2 

SB=S3 

Z12=P11+P21 

RETURN 

DO 100 1#1,3 

CI=1. 

IFC I. EQ 2)CI=-2. #CDK 

STI=S(1) 

TJ=TI1 

DO 90 J#1,2 

DZ=TJU-SI 

R=DSGQGRT(RHS+DZ#DZ) 

X=R+DZ 

IF (DZ. LT. . ODO) X#RHS/(R-DZ) 

CALL CISI(COSI.,CIN, SINI, X) 
EP=DCMPLX(COSI, -SINI) 

X=R-DZ 

IF (DZ. GT. . ODO) X=RHS/(R+DZ) 

CALL CISI(COSI,CIN, SINI, X) 

EM=DCMPLX (COSI, ~SINI) 

IF(J.EQ@.2)¢G0 TO 90 

EP1=EP 

EmM1=EM 

TJ=T2 
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100 


X=T2-SI 

EGDZ=DCMPLX (DCOS(X),. DSIN(X)) 
Z12=Z12+CI#( (EP-EP1 )*#EGDZ+(EM-EM1 )/EGDZ) 
CONTINUE 

Zi2= SGCN#ETA#Z12/(8. #PI*#SDK#SDKT ) 

RETURN 

END 
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SUBROUTINE ZSMM(S1, S2, Ti. T2. D. CDS, SDS, SDT, CPSI, P11) 
IMPLICIT REAL#@ (AH), (P-Z) 

DOUBLE PRECISION BD. CD. DD, DO. DZ. DPQ, DSG,. DPSI, CPSS 
DOUBLE PRECISION R,R1,R2, SDI. SIS, SK, SPSI,ST1, ST2 
DOUBLE PRECISION Tu, TD1, TD2. TL1., TL2, TS1, TS2, ZD 
COMPLEX#16 E(2,2),F(2,2),GAM,P11,P12, P21. P22 
COMPLEX#16 EGZI,ES1, ES2, ET1, ET2, EXPA, EXPB 
COMPLEX#16 EGDZ., EM, EP, EM1, EP1 

DATA ETA, GAM, PI/376. 730366239, (. ODO, 1. 1,3. 14159265359/ 
DD=D 

DP Q=DD*DD 

DPSI=CPSI 

CPSS=DPSI#DPSI 

IF (CPSS. GT. 1. DO) CPSS=1. DO 

SPSTI=DSQRT(1. DO-CPSS) 

SGDS=SDS 

IF (S2. LT. $1)SGDS*-SDS 

SGDT=SDT 

IF (T2. LT. T11)SGDT=-SDT 

IF (DABS(CPSI).LT. . 999999)G0O TO 6 
DO=SPS1I#*(T1+T2)/2. 

DSG=DP Q+D0#DO 

GO TO 110 

ES1=DCMPLX (DCOS(S1), DSIN(S1) ) 
ES2=DCMPLX(DCOS(S2), DSIN(S2) ) 

ET1=DCMPLX (DCOS(T1),DSINC(T1)) 

ET2=DCMPL.X (DCOS(T2), DSIN(T2) ) 

ID=1 

IF(D. E@. . ODO) ID=O 

TD1i=T1 

TD2=Te 

CD=DD/SPSI 

Cc=CD 

BD=CD*DPSI 

B=BD 

EB=. ODO 

EC=. ODO 

IF(CID E@ 0)GO TO 8 

EB=DEXP(-B) 

EC=DEXP(-C) 

DO 10 K=1,2 

DO 10 L=1,2 

E«K,L)=¢(. ODO, . ODO) 

TS1=TD1*TD1 

TS2=TD2*TD2 

SI=S1 

DO 100 I 1,2 

FI=(-1)##] 

SDI=SI 

SIS=SDI#SDI 

ST1i=e2. *SDI*#TD1L#DPSI 

ST2=2. *SDI#TD2*DPSI 
R1=DSQRT(DPQ@+SI1S+TS1-ST1 ) 

R2=DSQGRT (DPG+S1S+TS2-ST2) 

EK=EB 

DO SO K=1,2 

FA=(-1) ##K 

SKA=FK*SDI 
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EL=EC 
DO 40 L=1,2 
FL=(—-1) ##L 
EAL =EK#EL 
XX=2FK#BD+FL*CD 
TLI=FL#TD1 
TL2=FL#*#TD2 
RR1=R1+SK+TL1 
RR2=R2+SK+TL2 
AXX=DABS (XX) 
IF (AXX. GT. DABS(RR1)/100. )GO TO 28 
IF (AXX. GT. DABS(RR2)/100. )G60 TO 28 
IF (AXX. GT. .001)G0 TO 28 
IFCRR1/RR2.LT..0D0)GO TO 28 
CALL CISI(COS1, CIN, SIN1, RR1) 
CALL CISI(COS2, CIN, SIN2, RR2) 
EXPA=DCMPLX (COS2-COS1, SIN1-SIN2) 
ECA, L) SE CK, L) +FI#EXPA# (EKL4+1. /EKL) 
GO TO 40 
2868 CALL EXPJ(DCMPLX( XX, RR1), DCMPLX( XX, RR2), EXPA) 
CALL EXPJ( DCMPLX(-XX, RR1), DCMPLX(—-XX, RR2),. EXPB) 
E(K,L=E(K, L)+FI*# CEXPA#EKL+EXPB/EKL ) 
40 EL#1. /EC 
50 EK=1. /EB 
IFCI.E@. 2)6G0 TO 100 
ZD=SDI#DPSI 
ZC=ZD 
EGZI=DCMPLX(DCOS(ZC), DSIN(ZC) ) 
RR1=R1+ZD-TD1 
RR2=R2Z+ZD-TD2 
CALL CISI(COS1,CIN. SIN1, RR1) 
CALL CISI(COS2, CIN, SIN2, RR2) 
EXPB=DCMPLX (COS2-C0OS1. SIN1-SIN2) 
RR1=R1-ZD+TD1 
RR2=R2-ZD+TD2 
CALL CISI(COS1, CIN, SINI, RR1) 
CALL CISI(COS2, CIN. SIN2, RR2) 
EXPA=DCMPLX (COS2-COS1., SIN1I-SIN2) 
FCI, 1)=2. «SGDS#*(. ODO, 1. ) #EXPA/EGZI 
F(I,2)=2. #«SGDS#(. ODO, 1. ) #EXPB#EGZI 
100 SI=S2 
CST=-ETA/(16. #PI#SGDS#SGDT ) 
P11=CST#(( F(1,1)+E (2. 2)#ES2-E(1, 2) /ES2) #ET2 
A+(-F(1,2)-E(2, 1) #ES2+E (1, 1)/ES2)/ET2) 
RETURN 
110 IF(CPSI.LT. 0. »GO TO 120 
TA=T1 
TB=T2 
GO TO 130 
120 TA=-T1 
TB=-T2 
SGDT=-SGDT 
130 SI=S1 
CI=-cCDS 
P11=(. ODO. . ODO) 
DO 150 I=1,.2 
TJ=TA 
DO 140 J=1,2 
DZ=TU-SI 
R=DSGRT(DSG@+DZ#DZ) 
X=R+DZ 
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IF (DZ. LT. . ODO) X#DSG/ (R-DZ) 
CALL CISI(COSI. CIN. SINI, X) 
EP=DCMPLX (COSI, -SINI) 
X=R-DZ 
IF (DZ. GT. . ODO) X*DSG/ (R+DZ) 
CALL CISI(COSI., CIN. SINI,. X) 
EM=DCMPLX( COSI, -SINI) 
IF(J.EQ@.2)60 TO 140 
EP i1=EP 
EM1=EM 

140 TJ=TB 
X=TB-SI 
EGDZ=DCMPLX(DCOS(X),. DSIN(X)) 
Pi1=P11+CI#( (EP-EP 1) *EGDZ+(EM-EM1) /EGDZ) 
CI=1 

150 SI=S2 
PLL=ETA#P11/(8. #PI*SGDS*SCDT) 
RETURN 
END 


170 


RICHMD1 (28 of 30 pages) 


C 2 tt 2H Mt Hb a Ie HE EE HE OE ee HG EE Oe EE HE He aE HE EE OE a ee HE ae He HE HE EE SE EG EE HE aE aE HE EG aE ae aE HE a aE HE aE 
C SUBROUTINE SKEW bal 
C 2 ot tb at 3b Mt Hb aE Ie He TE IE He HE OE He aE HE HE HE aE HE aE HE ae SE 2G HE aE HE EE ae a aE Ea EE HE EE He Ea Ea a a a a OE a a a a 
SUBROUTINE SKEW (S1,S3, T1, T3, RHK, CDK, SDK, CDKT, SDKT, CPSI, Z12) 
IMPLICIT REAL#8 (A-H), (P-Z) 
COMPLEX#16 Z12, EIN, EGDZ,. CQX, EJXX, EM. EP 
COMPLEX#16 E(2, 2),F(2,2),ES1.ES2, ET1, ET2. EXPA,. EXPB, EGZI 
COMPLEX#16 P11,P12,P21,P22 
DIMENSION S(3),T(3) 
DOUBLE PRECISION R.SI,.SIS, TJ. TUS. XXD.S.T,.CI,CJU,FK, FL, DZ 
DOUBLE PRECISION RHS, DZS. DDZ. RIV. DPSI. CPSS, SPSI,. RH2 
DOUBLE PRECISION TD1,. TD2., SDI, ST1.ST2.R1.,R2, SK, ZD. TL1. TL2, SA. SB 
DATA ETA, PI. TP/376. 730366239, 3. 14159265359, 6. 28318530718/ 
$(1)=S1 
S2=(S1+S3) /a2. 
S(2)=S2 
$(3)=S3 
T(1)=T1 
T2=(71+T3) /2. 
T(2)=T2 
T(Q)=T3 
Z1i2=(. ODO. . ODO) 
DPSI=CPSI 
CPSS=DPSI#*DPSI 
IF (CPSS. GT. 1. DO) CPSS=1. DO 
SPSI#DSQRT(1. DO-CPSS) 
IF (DABS(CPSI). LT. . 999999)G0 TO 10 
RHS=RHK#RHA 
RH2=SPSI#(T1+T3) /2. 
RHS=RHS+RH2#RH2 
SGN=1 
IF(CPSI. GT. .O0D0O)GO TO 80 
SGN=-1. 
T(1)2-T3 
T(2)=-T2 
T(3)=-T1 
GO TO 80 
10 D=RHK 
DSG=D#D 
CD=D/SPSI 
BD=CD#DPSI 
EB=DEXP(-BD) 
EC=DEXP(-CD) 
CST=-ETA/ (16. #PI1#SDK#SDKT ) 
TA=T1 
TB=T2 
DO 70 ITT#1.2 
IFCITT. E@. 1)ET1=DCMPLX(DCOS(TA), DSIN(TA)) 
IF CITT. EQ. 2) ET2=DCMPLX(DCOS( TB), DSIN( TB) ) 
TD1=TA 
TD2=TB 
TS1=TD1#TD1 
TS2=TD2*TD2 
SA=S1 
SB=S2 
DO 60 I[SS=1,2 
IFC ISS. EQ. 1)ES1=DCMPLX (DCOS(SA). DSIN(SA) ) 
IF (ISS. EQ. 2)ES2=DCMPLX(DCOS(SB),. DSIN(SB)) 
DO 20 K=1,2 
DO 20 L=1,2 
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BO 


SI=SA 

DO 50 I=1,2 

FI=(-1) x] 

SDI=SI 

SIS=SDI#SDI 

STi=2. *SDI#TD1*#DPSI 

ST2=2. *SDI#TD2*DPSI 

Ri=DSQGRT (DSQ@+SIS+TS1-ST}2) 

R2=DSQGRT (DSQ+SIS+TS2-ST2) 

EK=EB 

DO 40 K=1, 2 

FA=(-1 )##K 

SK=FAK#SDI 

EL=EC 

DO 30 L=i.,2 

FL=(—-1)##L 

EAL=EK#EL 

XX=FA#BD+FL#CD 

TLI=FL#TD1 

TL2=FL#TD2 

RR1i=R1+SK+TL1 

RR2=R2+SK+TLS 

CALL EXPJ(DCMPLX (XX, RR1), DCMPLX( XX, RR2), EXPA) 
CALL EXPJ(DCMPLX(-XX, RR1), DCMPLX(-XX, RR2), EXPB) 
E(K,L)=E(K,L)+FI*#(EXPA#*EKL+EXPB/EKL ) 

EL=1. /EC 

EK=1. /EB 

IF (I. E@. ISS)GO TO $0 

ZD=SDI#DPSI 

ZC=ZD 

EGZI=DCMPLX(DCOS(ZC), DSIN( ZC) ) 

RRi=R1+ZD-TD1 

RR2=R2+ZD-TD2 

CALL EXPJ(DCMPLX(. ODO. RR1),. DCMPLX(. ODO, RR2), EXPB) 
RR1=R1-ZD+TD1 

RR2=R2-ZD+TD2 

CALL EXPJ(DCMPLX(. ODO. RR1), DCMPLX(. ODO, RR2), EXPA) 
FCI, 1)=(€. ODO. 2. )*SDA#EXPA/EGZI 
F(I,2)=(. ODO, 2. )*SDA#EXPB#EGZI 

SI=SB 

IFCITT. EQ. 2)6G0 TO 34 

IFC ISS. EQ. 1)P22=CST#( (F (2, 1)+E(2, 2)*#ES1-E(1,2)/ES1)#ETI1 


A+(-F (2, 2)-E(2, 1) #ES1+E(1,1)/ES1)/ET1) 


IF (ISS. EG. 2)P122CST#((-F (1, 1)-E(2. 2) #ES2+E(1,2)/ES2) #ET 1 


B+( F(1,2)+E(2, 1)*#ES2-E(1,1)/ES2)/ET1) 


GO TO 586 
IF CISS. EQ. 1)P21=CST#( (-F (2. 1)-E(2. 2) #ES1+E(1,2)/ES1 )#ET2 


C+( F(2, 2)+E(2, 1)*#ES1-E( 1. 1)/ES1)/ET2) 


IF (ISS. EQ. 2)P11=CST#( (F (1, 1)+E(2, 2)8ES2-E(1, 2)/ES2) #ET2 


D+(-F (1, 2)-E (2, 1) #ES2+E (1, 1)/ES2)/ET2) 


SA=S2 

SB=S3 

TA=T2 

TB=T3 
Z12=P11+P12+P21+P22 
RETURN 

DO 100 1=1,3 

ST=S(I) 

CIi=1. 

IF (I. EQ@. 2)CI=-2. #CDK 
CQ@x=(. ODO, . ODO) | 
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DO 90 J=1,3 

TUZ=T (J) 

CUri1. 

IF (J. EQ. 2)CUe-2. #CDKT 

DZ=TUJ-SI 

R=DSGRT (RHS+DZ+#DzZ ) 

X=R+DZ 

IF (DZ. LT. . ODO) X=RHS/(R-DZ) 
CALL CISI(COSI., CIN. SINI, xX) 
EP=DCMPLX (COSI, -SINI) 

X=R-DZ 

IF (DZ. GT. . ODO) X=RHS/ (R+DZ) 
CALL CISI(COSI,CIN. SINI, X) 
EM=DCMPLX(COSI, -SINI) 
EGDZ=DCMPLX (DCOS(DZ)., DSIN( DZ)? 
CQX=CQGX+C U# (EP#EGDZ+EM/EGDZ ) 
Zi2=Zi2+CI#CQx 

Zig= SGN#ETA#Z12/ (98. #PI*#SDK#SDKT ) 
RETURN 

END 
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PROGRAM "RICHMD2" 
At AE EE TE Ae HE a a a aE a a HE aE a ee a a a Ea a a a a a a 
* THIS COMPUTER PROGRAM. IN FORTRAN LANGUAGE, UTILIZES 
THE CURRENT DISTRIBUTIONS DETERMINED FROM THE “RICHMD1" 
PROGRAM TO OBTAIN THE RADIATION PATTERN OF A MONOPOLE 
ELEMENT OF LENGTH h AND RADIUS b AT THE CENTER OF A 
CIRCULAR GROUNDPLANE OF RADIUS a. AN INPUT FILE, 

INPUT1. DAT. IS AUTOMATICALLY GENERATED FROM THE 

“"RICHMD1" PROGRAM. 


THIS PROGRAM REGUIRES THE USE OF THE MMBSJO 
SUBROUTINE FUNCTION OF THE IMSL SOFTWARE LIBRARY. MMBSJ1 
IS THE MODIFIED BESSEL FUNCTION OF THE FIRST KIND AND 
ORDER ZERO. 
9G ET Te Se Tb aE a a HH a aE HE ae aE a aE aE a a a a a aE a a a a a aa a a a a a a a a a at a 
IMPLICIT REAL#8 (A-H), (P-Z) 

INTEGER#4 X, IER 

REAL#8 THETAD(89), THETA(89), PLOT(89),L 

REAL#8 BRKT(89), MBRAT(89), DIREC (89), DIRECD(89), DIRECR(89) 

REAL MMBSJO 

COMPLEX#16 C(50),E(50), ACUR1, ACUR2 

COMPLEX#1& ACUR( 300), UX, UY, U, USUM 

COMPLEX#16 SX, SY,S, SSUM 

COMPLEX#16 GP, GSUM, TSUM 

OPEN (UNIT=5, STATUS= ‘OLD’, FILE=’INPUT1. DAT’, READONLY ) 

OPEN(UNIT=6, STATUS= ‘NEW’, FILE=’ TEMP “ ) 

CH HE aE HE aE Se HE aE ETE EH UE a ab EE OE a aE aE HE aE a aE a a a a a a aH a Od a a a ae 

C READ AND WRITE STATEMENTS ial 

Ct HE Se a EE ae EE ETE ETE EE EH EI EE Oa IE HE IE EE IE aE aE a aE a Oa a aE a a a a a a a 
READ(S, 8) AL. H. EPSLN, NSD, NSW, N, X 


*eaeuue KR KKK EK eK 


ANAaAANANAAANRANAANH 
x“eeenenuagR ERK Ee HK & 


8 FORMAT(E12. 4, 2X, F8. & 2X, F7. 4, 2X, 4412, 2X) ) 
FLAG=1.0 
GO TO 70 

71 WRITE(&, 74) 

74 FORMAT(/,1X, “INPUT GROUNDPLANE CURRENTS. I ’) 
WRITE(&, 75) 

73 FORMAT(1X, “ RECI) %, 7X, ‘IMCI) %) 


DO 1 I=1,N-1 
READ(S,2) C(I) 

2 FORMAT (2E12. 4) 
IF (1. EQ. 1) THEN 
Ec I)=-Cc(I) 
ELSE 
EcI)=#C(I) 

1 CONTINUE 
DO 77 I#1,NSD 
WRITE(G6, 76)EC(1), 1 


7& FORMAT (2E12. 4, 2X, ‘AT RHO “, 12) 
77 CONTINUE 
WRITE (6, 78) 
78 FORMAT (/, 1X. ‘INPUT MONOPOLE CURRENTS. I°%) 


WRITE (6, 75) 
DO S55 I=N, NSD+1,-1 
C(Id)=CCI-1) 
oe Bbw Bo CONTINUE 
C (NSD+1)=C(1) 
DO 41 I=NSD+1,N 


MM=1-NSD 
WRITE (&, 42)C(1), MM 


42 FORMAT(2E12. 4. 2X, ‘AT Z ‘%, 12) 
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41 CONT INUE 
PI=3. 1415926536D0 
D=H/NSW 
DD=(EPSLN/(2#PI)—-AL) /NSD 
E(NSD+1)=DCMPLX(0. 0, 0. 0) 
ee eee. eee. ee eee eee eee eee eee ree ee ee ee Ce eee rer ee ee er ee ee 
C COMPUTE THE AVERAGE CURRENT IN EACH ANNULAR ZONE * 
ee ee ee ee oe ee ee eee cere eee eee ee er ee ee ce eo 
DELTA=DD/Xx 
K=1 
J=mtl 
P1=AL 
P2=AL+DELTA 
10 ACUR1=E (K) #(DCOS(2#P1#(K*#DD+AL-P2) )-DCOS(2*#PI#(K#*DD+AL-P 1) )) 
ACUR2=E (K+1) #(DCOS(2*PI#(P2-(K-1)#DD-AL ) )-DCOS (2#P1#(P1—(K—-1)#DD-AL) )) 
ACUR (J) =(ACUR1—-ACUR2) / (2#P I#DELTA#DSIN(2#PI#DD) ) 
JmJ+1 
P1=P1+DELTA 
P2=P2+DELTA 
G=(U-1)/ (K#X) 
IF(G.EG@.1) K=K+1 
IF(K.LE.NSD) GO TO 10 
Ck. ee eee ee ee ee eee eee ee eT PPP PC Pee eee eee ee 2 
C CONVERT DEGREE INTO RADIAN # 
C 9 He He He He ee He He He He He He He Hee He HE He He He a a a a HH HH a He SH a eH 2 HH 2 
DO 12 II=1,89 
K=II#2 
THETAD( 11) =DFLOAT(K) 
THETA(II)=PI#THETAD(II)/180. ODO 


12 CONTINUE 
C9 46 tb HE ae aE ae aE aE aE aE SE a ea a EE EE a aE a aE Ee EE aE a a a a aE a a Ee 
C COMPUTE WITH RESPECT TO EVERY TWO DEGREE * 


C HH Hb He He ee He He He He EH eH ae a ae He a ae He a He a a a ee 
DO 30 INDEX=1, 89 

C He 9 ee ee He a He He He HE HE He a a a a HH HH a 2 2 a a a a 

C FIELD DUE TO THE MONOPOLE % 

C He He Me ee EERE 
USUM=DCMPLX (0. 0, 0. O) 
NN=1 

40 UXR=DSIN(NN#2#P I #D#DCOS ( THETA( INDEX) ) ) 

UX IM=DCOS ( NN#2#P I #D#DCOS ( THETA( INDEX) ) ) 
UX=DCMPLX (-UXR, UXIM) 
UYR1=1-DCOS(2#P1#D) #DCOS (2#P I#D#DCOS( THETA( INDEX) ) ) 
UYR2=DCOS( THETA( INDEX) ) #DSIN(2#P1#D) #DSIN(2#P I #D#DCOS( THETA ( INDEX) )) 
UYR=UYR1-UYR2 
UYIM1=DCOS(2#P1#D) *DSIN(2#P 1#D#DCOS ( THETA( INDEX) )) 
UY IM2=DCOS( THETA ( INDEX) ) #DSIN(2#*P1#D) #DCOS ( 2#P I #D#DCOS( THETA( INDEX) ) ) 
UY IM=UYIM1-UYIM2 
UY=DCMPLX(UYR, UYIM) 
UsUX#UY#C (NN+NSD) 
USUM=USUM+U 
NN=NN+1 
IF(NN. LE. NSW) GO TO 40 
USUM=USUM/DSIN(2#PI#D) 

C He He Hee He He He HE ee He HE He ae ee ae He a He a HH a a 0 a a a a a 
SSUM=DCMPLX (0. 0, 0. 0) 
MM=2 

50 SXR=DSIN(2#P 1#D#(MM—2) #DCOS( THETA( INDEX) ) ) 

SX IM=DCOS( 2#P 1#D#(MM—2) #DCOS( THETA( INDEX) )) 
SX=DCMPLX (SXR, -SXIM) 
SYR1=DCOS(2#P1#D) #DCOS (2*#P I1#D#DCOS( THETA( INDEX ) ) ) 
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SYR2=DCOS( THETAC INDEX) ) *#DSIN(C2*P1#D) *¥DSIN(2#P1*#D#DCOS( THETACINDEX)) )-1 
SYR=SYR1+SYR2 
SYIM1=DCOS(2#P1I#D) #*DSIN(2*P1I#D#DCOS( THETAC INDEX) )) 
SY IM2=DCOS ( THETA( INDEX) )#DSIN(2#PI#D) #DCOS (2*P 1#D#DCOS ( THETAC INDEX) ) ) 
SY IM=SYIM1-SYIM2 
SY=DCMPLX (SYR, SYIM) 
S=C (NSD+MM ) #SY#SX 
SSUM2SSUNMN+S 
MM=MM+1 
IF (MM. LE. NSW) GO TO 50 
SSUM=SSUM/DSIN(2*#PI#D) 
CO aE HE AE AE AE AE TE HE AE TE HE HE AE AE AE HE HE HE HE IE HE HE HE HE IE HE IH Ae HE A aE HE EH I aE HEH a a a aH HH a Ot 
G FIELD DUE TO THE GROUNDPLANE #* 
C4 4 OE aE HE He aE HEE IEE OE AE IE a OE aE EOE AE EE HE aE EE A HE EH EE A AE HE I HE a HH aH a aH a a a at 
GSUM=DCMPLX (0. 0,0. 0) 
M=1 
P1=AL 
P2=AL+DELTA 
60 X1i=2e#PIl#P1#DSIN( THETACINDEX) ) 
Bi=MMBSJO(X1, IER) 
X2=24Pl#P2e#DSIN( THETAC INDEX) ) 
B2=MMBSJO(X2, IER) 
GP=DCOS( THETA( INDEX) )#ACUR(M) #(B1i-B2) 
GSUM=GSUM+GP 
Pi=P1+DELTA 
P2=P2+DELTA 
M=M+1 
IF(M. LE. (NSD*¥X)) GO TO 60 
C96 OE SE 3 aE HE OE aE He IE A OE AEE aE OE SE EAE HE IE HE HE aE OE OE EHH EH IEE HE OE a a aE Ea aE aaa aa a a aH a 


C COMPUTE THE RADIATION RESISTANCE AND DIRECTIVE GAIN #* 
C 4 3 St a St HE aE aE aE aE EE HE UE A HE aE TE HE HE HE HE HE HE UE IE EH HE I HE IE HE a EH EE aaa a a aaa a a aaa at 
TSUM=GSUM+USUM+SSUM 
RP=DREAL (TSUM) 
P=DIMAG(TSUM) 


BRKT C INDEX) =RP#RP+P#P 
MBRKT CINDEX ) =BRAT( INDEX) /DSIN( THETA( INDEX) ) 

30 CONT INVE 
DENOM=0. ODO 
DO 3 K=1,89 
DENOM=DENOM+MBRKT (HK) 

re CONT INUE 
DE NOM=DENOM*# THE TA(C1 ) 
DIRMAX=0. ODO 
DO 4 K=1,89 
DIREC (K)=2. ODO#BRKT (K) /C(CDSIN(CTHETACK) ) )##2. ODO#DENOM) 
DIRECD(K)=10. ODO#DLOGIO(CDIREC(K) ) 
IF (DIRMAX. LT. DIRECD(K) )DIRMAX=DIRECD(K) 

+ CONTINUE 
DO 15 J=1,89 
DIRECR(J)=DIRECD(J)-DIRMAX 

15 CONT INVE 
RADRES=15. 0#DENOM/ (CDABS(C (NSD+1) )) ##2 

C0 EOE Se He AE HE EE aE a aE aE ae AE AEE aE aE aE aE a Ha HE a HE a HH a a HH a a a aa a aa a Ha a Ha a aaa aaa aaa HE 


C WRITE STATMENTS ” 
C0696 Ob SESE SE Sb St He SE AE aE OE a OE EOE AE a Ea EI OE a I aE aE aE OE Aa A aE A aE HEH a aE aE aE a aa aa a aa aaa aa aH aE at 
7O WRITE(6. 80) EPSLN 
B80 FORMAT(1H1, “EPSLN= ‘’,F7. 4 ‘(NORMALIZED GROUNDPLANE RADIUS) ’) 
WRITE(46, 81) H 
81 FORMAT (1X, ‘MONOPOLE LENGTH= ’,F8. 6 ’ (WAVELENGTH) ‘) 
WRITE (6, 82) AL 
B2 FORMAT(1X, ‘MONOPOLE RADIUS= ‘’,E12. 4, ’ (WAVELENGTH) “) 
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WRITE (46, 83) NSW 
83 FORMAT( 1X, ‘NUMBER OF MONOPOLE SEGMENTS= ‘, I2) 
WRITE(6,84) NSD 
84 FORMAT( 1X, ‘NUMBER OF DISK SEGMENTS= ‘, I2) 
FLAG#FLAG+1. 0 
IF (FLAG. LE. 2)G0 TO 71 
WRITE (6, 73)X 
73 FORMAT(1X, ‘NUMBER OF ANNULAR ZONES PER DISK SEGEMENT= ‘, I2) 
WRITE (6,96) RADRES 
96 FORMAT( 1X, “RADIATION RESISTANCE DETERMINED BY PATTERN#= ‘, 
1 E12. 6, ° OHMS’) 
WRITE (6, 47) 


47 FORMAT(//,1X, ‘ELEVATION’, 53X, ‘ELEVATION ‘) 
WRITE(4, 85) 
85 FORMAT (1X, 2( “ANGLE ’, 6X, ‘DIRECTIVITY’, SX, ‘DIRECTIVITY’, 5X 
1 » “RELATIVE POWER’, 5X)) 
WRITE( 6. 86) 
86 FORMAT (1X, 2(°(DEG) °, 6X. ‘C(NUMERIC) ‘, 7X. ‘CdBid’. 11X 
1 » ‘(dB) °, 15X)? 
WRITE (6, 87) 
87 FORMAT(//,2X, ‘ O. 4 9X, ‘0. O0O000%, SX. ‘INFINITY’, 6X, “~INFINITY ’) 


DO 16 L=1, 44 
WRITE(6,. 88) THETAD(L), DIREC(L). DIRECD(L), DIRECR(L) 


1 » THETAD(L+45)., DIREC (L+45), DIRECD(L+45)., DIRECR(L+45) 
88 FORMAT (1X, 2(01X,F5. 1, 6X, F9. 5, SX, F9. 5S, OX, FF. 5, 11X)?) 
16 CONTINUE 


WRITE(&, 89) THETAD( 45), DIREC (45), DIRECD( 45), DIRECR(45) 
89 FORMAT (1X, 1X, FS. 1, 6X, F9. 5, SX, F9. 5, OX, FF. 5, 12X, 6180. ° 
1 .9X, ‘0.00000’, 5X. ‘~INFINITY’, 6X, ‘~ INFINITY ‘) 
23 CONTINUE 
CLOSE (UNIT#=5, STATUS# ‘KEEP “ ) 
CLOSE (UNIT=6, STATUS ‘KEEP *” ) 
END 
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CHRKEKKKKEKKEKKKKKEKEKEKKEKEKKEKKEKEKKEEKEKEEKKKEKKKKKKKKKKKEEKEEKEKEEKEEKEKKKEEEKETEE 


C * 
C PROGRAM '"MONOPL"” * 
C * 
CHRKKKAKKKKEKKKEEKKEKEKEKEEKKEEKKEEKKKEKEKEKEKEKEKEKEKEKEEKREEEKEEEKEKEKKEKKKEEKKE 
C * 
THIS PROGRAM IS WRITTEN IN FORTRAN FOR A DEC PDP 11/70. * 

IT CALCULATES, BY THE METHOD OF OBLATE SPHEROIDAL WAVE FUNCTIONS, * 

THE DIRECTIVITY (D,DDBI) AND THE RADIATION RESISTANCE * 
(RADRES) IN OHMS FOR A QUARTERWAVE MONOPOLE ANTENNA MOUNTED * 

AT THE CENTER OF CIRCULAR GROUNDPLANE OF RADIUS a IN FREE * 
SPACE, WHERE a IS OF THE ORDER OF A WAVELENGTH. * 

* 


KEKKKEKKEKKEKEKEEKEKEKEEKEKEKKEKEKEKEKKEKEKKEKKKEKKEKKEKKEEEKKEEKKKEKEKEEKKE 


THE INPUT VARIABLE IS: € (EPSLN) GIVEN BY 

€ = 277a/) (NORMALIZED GROUNDPLANE RADIUS) 
WHERE A IS THE EXCITATION WAVELENGTH. THIS IS ENTERED IN THE 
MAIN PROGRAM ON LINE 12. 


+ + + + H+ 


KEKEKEKKEKEKEKEKEKEEEKEKEKEKEKEKKKEKKEKEKEKKKEKKKKKEEKREEKKKEKKEKKKKKKEKKKKK 


THE OUTPUT VARIABLES ARE: 


RADRES -- THE ANTENNA’S RADIATION RESISTANCE (OHMS) 
D -- THE ANTENNA’S DIRECTIVITY (NUMERIC) 

DDBI -- THE ANTENNA’S DIRECTIVITY (DB) 

THETA -- ANGLE FROM THE ZENITH (RADIANS) 

THETAD -- ANGLE FROM THE ZENITH (DEGREES) 


++ ¢+ + + & H 4 F 


REFERENCES ; 


1) A. LEITNER AND R.D. SPENCE,”THE OBLATE SPHEROIDAL WAVE 
FUNCTIONS”, JOURNAL OF THE FRANKLIN INSTITUTE, 
VOL. 249, NUMBER 4, APRIL 1950, PAGES 299-321. 


2) A. LEITNER AND R.D. SPENCE, "EFFECT OF A CIRCULAR 
GROUNDPLANE ON ANTENNA RADIATION”, JOURNAL OF APPLIED 
PHYSICS, VOL. 21, OCTOBER 1950, PAGES 1001-1006. 


3) M.M. WEINER, S.P. CRUZE, C.C. LI, W.J. WILSON, *MONOPOLE 
ELEMENT AT THE CENTER OF A CIRCULAR GROUNDPLANE OF 
ARBITRARY RADIUS", MITRE TECHNICAL REPORT MTR-9622, 
SECTION 3.6, MAY 1985. 


+ + + + & & F HF & HF HH HH FH FH HF FH 


KEKKAKEKKEKEKEKEKEKEEKEKEKKKEKEKEKEEKEKKEKEKEKKEKEKEKEKEKEKREKKEKKKKKKEKEKEK 


THIS PROGRAM CONSISTS OF THE MAIN ROUTINE AND THE FOLLOWING 
SUBROUTINES: GAMMA, AALPHA, CEE, ANGLAR, FACTOR, ODDEVN, 
NORMAL, QUE, FACTR2, BAETA, QU, RAD30, DBLYU2, BEE, 
DIRECT, AEY, BBEE, GGEE, RADI, RAD2, RADRES, IMPROV, 

ODD, EVEN. 


+ + + & eH F 


C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
CHARI EKER TATE KEKE 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 


KEKKKEEKKKKEEKKEKEKKEKKEEKKEKEKKKKKKEKKKEKEEKEEEEKKKEKEKEKKEKKKEKEKKKKEKEKKKKK 
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** MAIN ROUTINE ** 


AOaAaAANH 


COMPLEX*8 Q1(20,1),B(20),V30( 20,1) 

REAL*8 U(20,1),NUE,C(40,20,1),V1(20,1), PSI, BETA( 20, 1) 
REAL*8 GAM(20,1),F(8,20,1),ALPHA(20,1),EPSLN,N( 20, 1) 
REAL*8 Q2(20,1),V2(20,1),W(20,1),D,A(40, 20,1) 

REAL*8 BB(40, 20,1),B0(20,1),G6(20,1),R,V22( 20,1), GAMIMP( 20,1) 
REAL*8 INCR, THETA, DDBI , THETAD 

INDEX=10 

JNDEX=24 

KNDEX=1 2 

LNDEX=18 

MNDEX=18 


INPUT: 


aa Qe 


EPSLN=3.0D0 


€2a' <3 


WRITE(4,4) EPSLN 
WRITE(4, 5) 
5 FORMAT(1X,”° ANGLE FROM ZENITH(DEGREES)* ,3X,°ANGLE FROM ZENITH 
* (RADIANS)* ,3X, “DIRECTIVITY (DBI) * ,3X, “DIRECTIVITY (NUMERIC) “/) 
INCR=1.0D0 
2 CONTINUE 
THETA=INCR*3.1415926536D0/90.0D0 
THETAD=THETA* 180.0D0/3.1415926536D0 
NUE=DCOS (THETA) 
C 
C CALL SUBROUTINES 
Cc 
CALL GAMMA(F,GAM,EPSLN, INDEX) 
CALL IMPROV(GAM,GAMIMP, EPSLN, INDEX) 
CALL AALPHA(ALPHA, GAMIMP, INDEX) 
CALL CEE(C,GAMIMP,EPSLN, INDEX, JNDEX) 
CALL ANGLAR(C,U,NUE, INDEX, KNDEX) 
CALL NORMAL(C,N, INDEX ,KNDEX) 
CALL QUE(Q1,C,EPSLN, INDEX ,KNDEX) 
CALL BAETA(C, BETA, INDEX ,KNDEX) 
CALL QU(Q2, BETA,Q1,EPSLN, INDEX) 
CALL AEY(A,ALPHA,EPSLN, INDEX, LNDEX) 
CALL RADI(EPSLN,A, V1, INDEX, LNDEX) 
CALL BBEE( BB,EPSLN,ALPHA,A, BETA, BO, INDEX, LNDEX, MNDEX) 
CALL GGEE(EPSLN, BB,G, BO, INDEX, MNDEX) 
CALL RAD2(V2,EPSLN,G,V1,Q2, INDEX) 
CALL RAD30( BETA,Q2, V30, INDEX) 
CALL DBLYU2(W,EPSLN,V1,V2,Q2, BETA, INDEX) 
CALL BEE(B,W,ALPHA,N,Q1, V30, INDEX) 
CALL DIRECT(B,U,EPSLN,D,N, INDEX, DDBI ) 
CALL RADRES(B,N,EPSLN,R, INDEX) 


a 


OUTPUT: 
WRITE(4,1) THETAN, THETA, DDBI ,D 


IF(INCR.EQ.89.0D0) GOTO 3 

INCR=INCR+1.0D0 

GOTO 2 

FORMAT(///1X,° RADIATION RESISTANCE = °,F10.5,1X,° OHMS”) 
FORMAT(1X,°EPSILON = *,F5.2///) 
FORMAT(10X,F5.1,24X,F7.5,16X,F9.5,11X,F9. 5) 

CONTINUE 

WRITE(4,6) R 

STOP 

END 
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é 

C ** SUBROUTINE GAMMA ** 

C 

C THIS SUBROUTINE CALCULATES THE RAW ELCENVALUES FROM EQUATION 21, 
C REFERENCE 1. THE F’S CORRESPOND TO THE F’S FROM EQUATIONS 22,23,24,AND 
C 25. THE H CORRESPONDS TO THE H FROM EQUATION 26. 

C THIS SUBROUTINE IS CALLED FROM MAIN AND DOES NOT CALL ANY OTHER 
C SUBROUTINES. 

C INPUTS: EPSLN, INDEX 

C OUTPUT: GAM 

c 


SUBROUTINE GAMMA(F, GAM, EPSLN, INDEX) 
REAL*8 F(8,20,1),GAM(20,1),H,A,B,EPSLN,AJ 
REAL*8 Z1,22,23,24,25,26,27,28 
H(A)=(A*A-1.0D0) /(4.0D0*A*A-1.0D0) 
DO 4 N=1, INDEX 
CAM(N,1)=0.0D0 
4 CONTINUE 
B=1.0D0 
DO 1 L=1, INDEX 
F(2,L,1)=2.0D0*(B*(B+1.0D0) )/((2.0D0*B-1.0D0)*(2.0D0*B+3.0D0) ) 
F(4,L,1)=H(B-1.0D0)*H(B)/(2.0D0*(2.0D0*B-1.0D0) )- 
*H(B+1.0D0)*H( B+2.0D0) /(2.0D0*(2.0D0*B+3.0D0)) 
F(6,L,1)=(3.0D0/((2.0D0*B-1.0D0)*(2.0D0*B+3.0D0) ) )*(H(B-1.0D0) 
**H(B)/(((2.0D0*B-1.0D0)**2.0D0)*(2.0D0*B-5.0D0) )-H( B+1.0D0) 
**H(B+2.0D0)/(((2.0D0*B+3.0D0)**2.0D0)*(2.0D0*B+7 .0D0) ) ) 
Z1=H(B-1.0D0)*H(B) /(2.0D0*(2.0D0*B-1.0D0)) 
Z2=36.0D0/((2.0D0*B-5.0D0)**2.0D0*(2.0D0*B-1.0D0)**4.0D0 
** (2. 0D0*B+3.0D0)**2.0D0) 
Z23=F(4,L,1)/(2.0D0*(2.0D0*B-1.0D0) ) 
Z4=H(B-3.0D0)*H(B-2.0D0)/(8.0D0*(2.0D0*B-1.0D0)* 
*(2.0D0*B-3.0D0) ) 
Z5=H(B+1.0D0)*H(B+2.0D0) /(2.0D0*(2.0D0*B+3.0D0) ) 
26=36 .0D0/((2.0D0*B-1.0D0)**2.0D0*(2.0D0*B+3.0D0)**4.0D0 
**(2.0D0*B+7.0D0)**2.0D0) 
Z7=F(4,L,1)/(2.0D0*(2.0D0*B+3.0D0) ) 
Z8=H(B+3.0D0)*H(B+4.0D0) /(8.0D0*( 2.0D0*B+3.0D0)* 
*(2.0D0*B+5.0D0) ) 
F(8,L,1)=Z1*(2Z2-23+Z4) -Z5*(Z64+Z7+Z8) 
B=B+1.0D0 
1 CONTINUE 
DO 2 L=1, INDEX 
AJ=2.0D0 
DO 3 J=2,8,2 
CAM(L,!)=GAM(L,1)+F(J,L,1)*EPSLN*®*AJ 
AJ=AJ+2.0D0 
3 CONTINUE 
CONTINUE 
RETURN 
END 
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** SUBROUTINE AALPHA ** 


THIS SUBROUTINE CALCULATES ALPHA FROM EQUATION 15, REFERENCE 1. 
GAM IS THE REFINED EIGENVALUE WHICH IS INPUT TO THE SUBROUTINE. 
THIS SUBROUTINE IS CALLED FROM MAIN. 

INPUTS: GAM, INDEX 

OUTPUT: ALPHA 


SUBROUTINE AALPHA(ALPHA, GAM, INDEX) 
REAL*8 ALPHA( 20,1),GAM(20,1),AL 
AL=1.0D0 
DO 1 L=1, INDEX 
ALPHA(L,1)=AL*(AL+1 .0D0)+GAM(L, 1) 
AL=AL+1.0D0 

1 CONTINUE 
RETURN 
END 


** SUBROUTINE CEE ** 


THIS SUBROUTINE CALCULATES THE C’S IN EQUATIONS 16 AND 17 IN 
REFERENCE 1. 

THIS SUBROUTINE IS CALLED FROM MAIN AND DOES NOT CALL ANY OTHER 
SUBROUTINES. 

INPUTS: GAM,EPSLN, INDEX, JNDEX 

OUTPUT: C 


SUBROUTINE CEE(C,GAM, EPSLN, INDEX, JNDEX) 
REAL*8 C(40,20,1),GAM(20,1),AL,AJ,EPSLN 
AL=1 .0D0 
DO 1 L#1,(INDEX-1),2 
C(2,L,1)=((1.0D0-AL)*(2.0D0+AL)-GAM(L,1))/8.0D0 
C(4,L,1)=(((3.0D0-AL)*(4.0D0+AL)-GAM(L,1))*C(2,L,1)+ 
*EPSLN**2.0D0)/24.0D0 
AJ=6.0D0 
DO 2 J=6, INDEX, 2 
C(J,L,1)=(((AJ-1.0DO-AL)* (AJ+AL) -GAM(L,1))*C(J-2,L,1)+ 
*C(J-4,L,1)*EPSLN**2.0D0) /(AJ*(AJ+2.0D0) ) 
AJ=AJ+2.0D0 
2 CONTINUE 
AL«AL+2.0D0 
1 CONTINUE 
AL=2.0D0 
DO 3 L=2, INDEX, 2 
C(2,L,1)=((2.0D0-AL)*(3.0D0+AL)-GAM(L,1))/8.0D0 
C(4,L,1)=(€((4.0D0-AL)*(5.0D0+AL)-GAM(L,1))*C(2,L,1)+ 
*EPSLN**2.0D0)/24.0D0 
AJ=6.0D0 
DO 4 J=6,JNDEX,2 
C(J,L,1)=(((AJ-AL)*(AJ+1.0D0+AL) -GAM(L,1))*C(J-2,L,1) 
*+C(J-4,L,1)*EPSLN**2.0D0) /(AJ*(AJ+2.0D0) ) 
AJ*AJ+2.0D0 
4 CONTINUE 
AL=AL+2.0D0 
3 CONTINUE 
RETURN 
END 
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€ 
c ** SUBROUTINE ANGLAR ** 
Cc 
C THIS SUBROUTINE CALCULATES THE ANGULAR SPHEROIDAL WAVE FUNCTIONS 
C AS PERK EQUATION 13 AND 14, REFERENCE 1. THE ANGULAR SPHEROIDAL 
C WAVE FUNCTIONS ARE REPRESENTED BY U AND THE COEFFICIENTS BY C. 
C THIS SUBROUTINE IS CALLED FROM MAIN AND DOES NOT CALL ANY SUBROUTINES. 
C INPUTS: C,NUE, INDEX,KNDEX 
C OUTPUT: U 
Cc 
SUBROUTINE ANGLAR(C,U,NUE, INDEX,KNDEX) 
REAL*8 C(40,20,1),U(20,1),NUE,AJ 
DO 3 N=], INDEX 
U(N,1)=0.0D0 
3 CONTINUE 
DO 1 L#=1,(INDEX-1),2 
AJ=1.0D0 
DO 2 J=],KNDEX 
UCL, 1)=UC(L,1)+C(2*J,L,1)*(1.0DO-NUE**2.0D0)**AJ 
AJ=AJ+1.0D0 
2 CONTINUE 
UCL, 1) =DSQRT(1.O0D0-NUE**2.0D0)*(1.0D0+U(L,1)) 
1 CONTINUE 
DO 4 L=2,INDEX,2 
AJ=1.0D0 
DO 5 J=1,KNDEX 
UCL, 1)=UC(L, 1)+C(2*J,L,1)*(1.0DO0-NUE**2.0D0)**AJ 
AJ=AJ+1.0D0 
5 CONTINUE 
U(L,1)=DSQRT(1.0D0-NUE** 2.0D0)*NUE*(1.0D0+U(L,1)) 
4 CONTINUE 
RETURN 
END 
Cc 
Cc ** SUBROUTINE FACTOR ** 
Cc 
© THIS SUBROUTINE CALCULATES THE FRACTIONAL PART OF EQUATIONS 28 AND 
C 29 IN REFERENCE 1. 
C THIS SUBROUTINE IS CALLED FROM NORMAL AND DOES NOT CALL ANY 
C SUBROUTINES. 
C INPUTS: 1,J,1TERM,AI,AJ 
C OUTPUTS: AF, DENOM,B 
Ec 


SUBROUTINE FACTOR(ITERM,AI ,AJ,DENOM,AF,1I,J,8) 
REAL*8 AI,AJ,DENOM,AF,AN,B,AM 
AF=].0D0 
DENOM=1 .0D0 
AN=1.0D0 
DO 1 N=1,(1l+I+J) 
AF=AF*AN 
AN=AN+1] .0D0 
1 CONTINUE 
AM=1.0D0 
DO 2 M=1,(2*(1+J+1)+ITERM) ,2 
DENOM=DENOM*AM 
AM=AM+2.0D0 
2 CONTINUE 
B=2.0D0**(2.0D0+AI+AJ) 
RETURN 
END 
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** SUBROUTINE ODDEVN ** 


THIS SUBROUTINE FIGURES OUT WHETHER A ONE OR A THREE BELONGS IN THE 
DENOMINATOR OF EQUATIONS 28 AND 29 OF REFERENCE 1. THE PROGRAM INPUTS 
L, DETERMINES WHETHER L~-] IS ODD OR EVEN AND OUTPUTS THE CORRESPONDING 
3 OR 1. 

THIS SUBROUTINE IS CALLED FROM NORMAL AND CALLS NO OTHER SUBROUTINES. 
INPUTS: L 

OUTPUTS: ITERM 


SUBROUTINE ODDEVN(L, ITERM) 
ITERM=1 

B=FLOATI(L) 

I=L/2 

C=B/2.0+0.6 

J=LINT(C) 

IF(J.EQ.1) ITERM=3 

RETURN 

END 


** SUBROUTINE NORMAL ** 


THIS SUBROUTINE CALCULATES THE NORMALIZATION FACTOR,N, FOR THE 
ANGULAR FUNCTIONS,U. THIS SUBROUTINE CALCULATES EQUATIONS 28 

AND 29 OF REFERENCE 1. 

THIS SUBROUTINE CALLS SUBROUTINES FACTOR AND ODDEVN, AND IS CALLED 
FROM MAIN. 

INPUTS: C, INDEX,KNDEX 

OUTPUTS: N 


SUBROUTINE NORMAL(C,N, INDEX,KNDEX) 
REAL*8 N1,AJ,AI,AF, DENOM,B,N2,N(20,1),C(40, 20,1) 
DO 4 L=], INDEX 
CALL ODDEVN(L,ITERM) 
[20 
AJ=1.0D0 
AI=0.0D0 
N1=0.0D0 
DO 1 J=1,KNDEX 
CALL FACTOR(ITERM,AI,AJ, DENOM,AF,I,J,B) 
Nl=N1+C(2*J,L,1)*B*AF/DENOM 
AJ=AJ+1.0D0 
1 CONTINUE 
Nl=2.0D0*N1 
AI=1.0D0 
N2=0.0D0 
DO 2 I=1,KNDEX 
AJ=1.0D0 
DO 3 J=]1,KNDEX 
CALL FACTOR(ITERM,AI,AJ,DENOM,AF,I1,J,B) 
N2=N2+C(2*J,L,1)*C(2*1,L,1)*B*AF/DENOM 
AJ=AJ+1.0D0 
3 CONTINUE 
AI=AI+1.0b0 
2 CONTINUE 
CALL FACTOR(ITERM,0.0D0,0.0D0, DENOM,AF,0,0,B) 
N(L,1)=AF*B/DENOM+N1+N2 
4 CONTINUE 
RETURN 
END 
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C ** SUBROUTINE QUE ** 

c 

C THIS SUBROUTINE CALCULATES THE JOINING FACTOR Q, USING EQUATIONS 
C 69 AND 70 IN REFERENCE 1. 

C THIS SUBROUTINE CALLS SUBROUTINE FACTR2 AND IS CALLED FROM MAIN. 
C INPUTS: C,EPSLN, INDEX,KNDEX 

€ OUTPUTS: Q] 

C 


SUBROUTINE QUE(Q1,C,EPSLN, INDEX,KNDEX) 
REAL IMPRT 
COMPLEX Q1( 20,1) 
REAL*8 AF,DENOM,B,AK,TERM,C(40, 20,1) 
DO 1 L=#1,(INDEX-1),2 
AK=] .0D0 
TERM=0 .0DO 
DO 2 K=]1,KNDEX 
CALL FACTR2(AK,K,DENOM,AF,B,!) 
TERM=TERM+C( 2*K,L,1)*AF*B/DENOM 
AK=AK+1.0D0 

2 CONTINUE 
CALL FACTR2(0.0D0,0, DENOM,AF,B,1) 
TERM=TERM+ B*AF/DENOM 
RLPRT#=1 .0/TERM 
Q1(L,1)=CMPLX(RLPRT,O.0) 

1 CONTINUE 
DO 3 L=2, INDEX, 2 
AK=] ..0D0 
TERM=0 .0D0 
DO 4 K=1,KNDEX 
CALL FACTR2(AK,K, DENOM,AF,B,3) 
TERM=TERM+C( 2*K,L, 1) *AF*B/DENOM 
AK=AK+1 .0D0 

4 CONTINUE 
CALL FACTR2(0.0D0,0,DENOM,AF,B, 3) 
TERM=TERM+ B*AF /DENOM 
IMPRT=-1 .0/ (EPSLN*TERM) 
Q1(L,1)=CMPLX(0.0, IMPRT) 

3 CONTINUE 
RETURN 
END 


** SUBROUTINE FACTR2 ** 


THIS SUBROUTINE CALCULATES THE FRACTIONAL PART OF EQUATIONS 69 
AND 70 IN REFERENCE 1. 

THIS SUBROUTINE IS CALLED FROM SUBROUTINE QUE, AND CALLS NO 
OTHER SUBROUTINES. 

INPUTS: AK,K,ITERM 

OUTPUTS: DENOM,AF,B 
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SUBROUTINE FACTR2(AK,K, DENOM,AF, B, ITERM) 
REAL*8 AK, DENOM,AF,B,AN 
AF=1 .0D0 
DENOM=1.0D0 
AN=1 .ODO 
DO 1 N=1,(1+4K) 
AF=AF*AN 
AN=AN+1 .0D0 

1 CONTINUE 
AM=1 .0D0 
DO 2 M=1,(2*(K+1)+ITERM) , 2 
DENOM=DENOM*AM 
AM=AM+2.0D0 

2 CONTINUE 
B=2.0D0**(1.0D0+AK) 
RETURN 
END 
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** SUBROUTINE BAETA ** 


THIS SUBROUTINE CALCULATES AN INTERMEDIATE RESULT BASED ON EQUATIONS 
44 AND 45 FROM REFERENCE 1. BETA CAN ALSO BE DERIVED FROM EQUATIONS 
84,87, AND 88 IN REFERENCE 1. 

THIS SUBROUTINE IS CALLED FROM MAIN, AND CALLS NO OTHER SUBROUTINES. 
INPUTS: C, INDEX, KNDEX 

OUTPUT: BETA 
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SUBROUTINE BAETA(C, BETA, INDEX, KNDEX) 

REAL*8 C(40,20,1),BETA(20,1),TERM 

DO 2 L=1,(INDEX-1),2 

TERM#1 .0D0 

pO 1 I#1,KNDEX 

TERM=TERM+C( 2*1,L,1) 
1 CONTINUE 

TERM=TERM*TERM 

BETA(L,1)=(0.5D0-2.0D0*C(2,L,1))*TERM 
2 CONTINUE 

DO 4 L=2, INDEX, 2 

TERM=1.0D0 

DO 3 1#1,KNDEX 

TERM=TERM+C( 2*1,L,1) 
3 CONTINUE 

TERM=TERM*TERM 

BETA(L,1)=(1.5D0-2.0D0*C(2,L,1)})}*TERM 
4 CONTINUE 

RETURN 

END 


** SUBROUTINE QU ** 


THIS SUBROUTINE EVALUATES EQUATIONS 87 AND 88 IN REFERENCE 1. 
THIS SUBROUTINE IS CALLED FROM MAIN AND CALLS NO SUBROUTINES. 
INPUTS: BETA,Q1,EPSLN, INDEX 

OUTPUT: Q2 
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SUBROUTINE QU(Q2, BETA, Ql, EPSLN, INDEX) 
REAL*8 Q2( 20,1), BETA(20,1),EPSLN 
COMPLEX Q1( 20,1) 


DO 1 L=1, INDEX 
Q2(L,1)#4.0D0* BETA(L,1)*(CABS(QI(L,1))**2.0)/(EPSLN**3.0D0) 


1 CONTINUE 
RETURN 
END 


** SUBROUTINE RAD30 ** 


THIS SUBROUTINE CALCULATES THE DERIVATIVE OF THE RADIAL FUNCTION 
OF THE THIRD KIND, EVALUATED AT ZERO. SEE EQUATION 57 REFERENCE 1. 
THIS SUBROUTINE IS CALLED FROM MAIN AND CALLS NO OTHER SUBROUTINES. 
INPUTS: Q2, BETA, INDEX 

OUTPUTS: ¥30 
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SUBROUTINE RAD30( BETA, Q2, V30, INDEX) 
REAL*8 BETA(20,1),Q2(20,1) 
COMPLEX V30( 20,1) 
pO 1 L#=1,(INDEX-1),2 
A=Q2(L,1)/BETA(L,1) 
V30(L,1)=CMPLX(0.0,A) 

1 CONTINUE 
DO 2 L=2, INDEX, 2 
B=-Q2(L,1)*1.570796327D0 
V30(L,1)=CMPLX(1.0,B) 

2 CONTINUE 
RETURN 
END 
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** SUBROUTINE DBLYU2 ** 


THIS SUBROUTINE EVALUATES EQUATION 26 OF REFERENCE 2, USING THE 
PORTION OF EQUATION 17, REFERENCE 2, WHICH IS IN BRACKETS. 

THIS SUBROUTINE IS CALLED FROM MAIN AND CALLS NO OTHER SUBROUTINE. 
INPUTS: EPSLN, V1,V2,Q2, BETA, INDEX 

OUTPUT: W 
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SUBROUTINE DBLYU2(W, EPSLN, V1] ,V2,Q2, BETA, INDEX) 
REAL*8 W(20,1),PS1,V1(20,1),V2(20,1) 
REAL*8 Q2(20,1),BETA(20,1),EPSLN,PI2 
PI2=1.570796327D0 
PS1=P12/EPSLN 
DO 1 L#],(INDEX-]),2 
W(L,1)=DSQRT(1.0D0+PSI*PSI)*(Q2(L,1)*VI(L,1)/ 
*BETA(L,1)) 
1 CONTINUE 
DO 2 L=2, INDEX, 2 
W(L,1)=DSQRT(1.0D0+PSI*PSI)*(-Q2(L,1)*1.570796327D0 
**V1(L,1)-V2(L,1)) 
2 CONTINUE 
RETURN 
END 


** SUBROUTINE BEE ** 


THIS SUBROUTINE CALCULATES THE B’S IN EQUATION 35 OF REFERENCE 2. 
THIS SUBROUTINE IS CALLED FROM MAIN AND CALLS NO SUBROUTINES. 
INPUTS: W,ALPHA,N,Q1,V30, INDEX 

OUTPUT: B 
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SUBROUTINE BEE(B,W,ALPHA,N, QI] , V30, INDEX) 

COMPLEX Q1(20,1),V30(20,1),B(20) 

REAL*8 W( 20,1) ,ALPHA(20,1),N(20,1) 

DO 1 L=],(INDEX-1) ,2 

B(L)=W(L,1)/(ALPHA(L,1)*N(L,1)*Q1(L,1)*V30(L,1)) 
1 CONTINUE 

DO 2 |.=2, INDEX, 2 

B(L)=-wW(L,1)/(ALPHA(L, 1)*N(L,1)*Q1(L,1)*V30(L,1)) 
2 CONTINUE 

RETURN 

END 


** SUBROUTINE DIRECT ** 


THIS PROGRAM CALCULATES THE DIRECTIVITY,D, AND THE DIRECTIVITY, 
DDBI, IN DBI, USING EQUATIONS FOR DIRECTIVITY DEVELOPED IN THE 
WORKING PAPER. THESE EQUATIONS, IN TURN, COME FROM EQUATIONS 43 AND 
45 OF REFERENCE 2. 

THIS SUBROUTINE IS CALLED FROM MAIN AND CALLS NO SUBROUTINES. 
INPUTS: U,EPSLN,N, INDEX,B 

OUTPUTS: D,DDBI 
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SUBROUTINE DIRECT(B,U,EPSLN,D,N, INDEX, DDBI) 
COMPLEX B(20),SI 
REAL*8 U(20,1),EPSLN,D,S,P,N(20,1) ,DDBI 
SI=CMPLX(0.0,0.0) 
DO 1 L=)], INDEX 
SI=S1+B(L)*U(L,1) 
1 CONTINUE 
S=(CABS(SI))**2.0D0 
P=0.0D0 
DO 2 J=1, INDEX 
P=P+N(J,1)*(CABS(B(J))**2.0) 
2 CONTINUE 
D=S/P*2.0D0 
DDBI=10.0D0*DLOG10(D) 
RETURN 
END 
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** SUBROUTINE AEY ** 


Oa 


C THIS SUBROUTINE EVALUATES THE A’S IN EQUATION 35 REFERENCE 1. 
C THIS SUBROUTINE IS CALLED FROM MAIN AND CALLS NO OTHER SUBROUTINES. 
, C INPUTS: ALPHA, EPSLN, INDEX, LNDEX 
C OUTPUT: A 
Cc 
SUBROUTINE AEY(A,ALPHA, EPSLN, INDEX, LNDEX) 
REAL*8 A(40,20,1),ALPHA(20,1),EPSLN,AN 
DO 1 L=#=1,(INDEX-1),2 
A(2,L,1)=(ALPHA(L, 1)-EPSLN*EPSLN-2.0D0)/2.0D0 
A(4,L,1)=( (ALPHA(L, 1) -EPSLN*EPSLN-12.0D0)*A(2,L,1) 
*-EPSLN*EPSLN) /12.0D0 
AN=6.0D0 
DO 2 N=6,LNDEX, 2 
A(N,L,1)=( (ALPHA(L, 1)-EPSLN*EPSLN-2.0D0-(AN-2.0D0) 
**(AN+1.0D0))*A(N-2,L,1)-EPSLN*EPSLN*A(N-4,L,1))/ 
*(AN*(AN-1.0D0)) 
AN=AN+2.0D0 
CONTINUE 
CONTINUE 
DO 3 L=2, INDEX, 2 
A(1,L,1)#1.0D0 
A(3,L,1)=(ALPHA(L, 1)-EPSLN*EPSLN-6 .0D0)/6.0D0 
AN=5.0D0 
DO 4 N=5,(LNDEX-1),2 
A(N,L,1)=( (ALPHA(L, 1) -EPSLN*EPSLN-2.0D0-(AN-2.0D0) 
**(AN+1.0D0))*A(N-2,L,1)-EPSLN*EPSLN*A(N-4,L,1))/ 
* (AN*(AN-1,0D0)) 
AN=AN+2.0D0 
4 CONTINUE 
3 CONTINUE 
RETURN 
END 
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** SUBROUTINE BBEE ** 


THIS SUBROUTINE EVALUATES THE B’S, HERE CALLED BB, IN EQUATION 50, 
REFERENCE 1. 

THIS SUBROUTINE IS CALLED FROM MAIN AND DOES NOT CALL ANY SUBROUTINES 
ITSELF. 

INPUTS: EPSLN,ALPHA,A, BETA, INDEX, LNDEX, MNDEX 

OUTPUTS: BB,BO 
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SUBROUTINE BBEE( BB, EPSLN,ALPHA,A, BETA, BO, INDEX, LNDEX, MNDEX) 
REAL*8 BB(40,20,1),EPSLN,AN,ALPHA( 20,1) ,A(40, 20,1) 
REAL*8 BETA(20,1),B0(20,1),DENOM 
DENOM=6 .0DO 
DO 1 L=1,(INDEX-1),2 
BB(1,L,1)=1.0D0/BETA(L,1)-1.0D0 
BB(3,L,1)=(-4.0D0*(A(2,L,1)+1.0D0)+BB(1,L,1)* 
* (ALPHA(L, 1)-EPSLN*EPSLN-4 .0D0) ) /DENOM 
AN=5.0D0 
DO 2 N=5,(MNDEX-1),2 
BB(N,L,1)=(BB(N-2,L,1)*(ALPHA(L, 1 )-EPSLN*EPSLN- 
*(AN-2.0D0)*(AN+1.0D0) )-2.0D0*(AN-1.0D0)* 
*(A(N-1,L,1)+A(N-3,L,1))-EPSLN*EPSLN*BB(N-4,L,1)) 
* / (AN*(AN-1.0D0)) 
AN=AN+2.0D0 
2 CONTINUE 
1 CONTINUE 
DO 3 L=2,INDEX,2 
BO(L,1)=-1.0D0/BETA(L,1) 
BB(2,L,1)=(-2.0D0+(ALPHA(L, 1) -EPSLN*EPSLN)*BO(L,1))/2.0D0 
BB(4,L,1)=(-6.0D0*(A(3,L,1)+1.0D0) -EPSLN*EPSLN*BO(L,1)+ 
*BB(2,L,1)*(ALPHA(L,1)-10.0D0-EPSLN*EPSLN))/12.0D0 
AN=6 .0D0 
DO 4 N=6,MNDEX, 2 
BB(N,L,1)=(BB(N-2,L,1)*(ALPHA(L, 1)-EPSLN*EPSLN- 
* (AN-2.0D0)*(AN+1 .0D0) )-2.0D0*(AN-1 .OD0)* 
*(A(N-1,L,1)+A(N-3,L,1))-EPSLN*EPSLN*BB(N-4,L,1)) 
* /(AN*(AN-1.0D0)) 
AN=AN+2.0D0 
4 CONTINUE 
3 CONTINUE 
RETURN 
END 
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** SUBROUTINE GGEE ** 


THIS SUBROUTINE CALCULATES THE G°S IN EQUATION 48 REFERENCE 1. 

THIS SUBROUTINE IS CALLED BY MAIN AND DOES NOT CALL ANY SUBROUTINES 
ITSELF. 

INPUTS: EPSLN, BB, BO, INDEX, MNDEX 

OUTPUT: G 


SUBROUTINE GGEE(EPSLN, BB,G, BO, INDEX, MNDEX) 
REAL*8 EPSLN, BB(40,20,1),G(20,1),AN, TERM, BO(20,1),PSI 
PSI#=3.1415926536D0/(2.0D0*EPSLN) 
DO 1 L=2,INDEX,2 
TERM=BO(L,1) 
AN=2.0D0 
DO 2 N=2,MNDEX,2 
TERM=TERM+BB(N,L,1)*(PSI**AN) 
AN=AN+2.0D0 
2 CONTINUE 
GCL, 1) =TERM/DSQRT(1.0D0+PS1**2.0D0) 
1 CONTINUE 
DO 3 L=1,(INDEX-1),2 
TERM=0 .ODO 
AN=1.0D0 
DO 4 N=1,(MNDEX-1),2 
TERM=TERM+ BB(N,L,1)*(PSI**AN) 
AN=AN+2.0D0 
4 CONTINUE 
G(L, 1)=TERM/DSQRT(1.0D0+PSI**2.0D0) 
3 CONTINUE 
RETURN 
END 


** SUBROUTINE RADI] ** 


THIS SUBROUTINE CALCULATES THE RADIAL FUNCTION OF THE FIRST KIND, 

Vl, USING EQUATION 34 IN REFERENCE 1. 

THIS SUBROUTINE IS CALLED FROM MAIN BUT DOES NOT CALL ANY SUBROUTINES 
ITSELF. 

INPUTS: EPSLN,A, INDEX, LNDEX 

OUTPUT: Vl 


SUBROUTINE RADI1(EPSLN,A,V1, INDEX, LNDEX) 
REAL*8 EPSLN,A(40,20,1),V1(20,1),AN, TERM, PSI 
PS1I=3.1415926536D0/ (2.0D0*EPSLN) 
DO 1 L=2,INDEX,2 
TERM=0 .0D0 
AN=1.0D0 
DO 2 N=1,(LNDEX-1),2 
TERM=TERM+A(N,L,1)*(PSI**AN) 
AN=AN+2.0D0 
2 CONTINUE 
V1(L, 1) =TERM*DSQRT(1.0D0+PSI**2.0D0) 
1 CONTINUE 
DO 3 L#=1,(INDEX-1),2 
TERM=1 .ODO 
AN=2.0D0 
DO 4 N=2,LNDEX,2 
TERM=TERM+A(N,L,1)*(PS1**AN) 
AN=AN+2.0D0 
4 CONTINUE 
V1(L,1)=TERM*DSQRT(1.0D0+PS1**2.0D0) 
3 CONTINUE 
RETURN 
END 
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** SUBROUTINE RAD2 ** 


THIS SUBROUTINE CALCULATES THE RADIAL FUNCTION OF THE SECOND KIND, 
V2, USING EQUATION 46, REFERENCE 1. 

THIS SUBROUTINE IS CALLED BY MAIN BUT DOES NOT CALL ANYTHING ITSELF. 
INPUTS: EPSLN,G,V1,Q2, INDEX 

OUTPUT: V2 


SUBROUTINE RAD2(V2,EPSLN,G, V1,Q2, INDEX) 
REAL*8 V2(20,1),6(20,1),V1(20,1),Q2(20,1),EPSLN,C, PSI 
PSI=3.1415926536D0/( 2.0D0*EPSLN) 
C=-3.1415926536D0/2.0D0 
DO 1 L=1, INDEX 
V2(L,1)=Q2(L,1)*(VI1CL,1)*(DATANC( PSI )4+C)+G(L,1)) 

1 CONTINUE 
RETURN 
END 


** SUBROUTINE RADRES ** 


THIS SUBROUTINE CALCULATES THE RADIATION RESISTANCE, R,USING 
EQUATION 45, REFERENCE 2. 

THIS SUBROUTINE IS CALLED BY MAIN AND DOES NOT CALL ANY 
SUBROUTINES ITSELF. 

INPUTS: N,EPSLN,B, INDEX 

OUTPUT: R 


SUBROUTINE RADRES(B,N,EPSLN,R, INDEX) 

REAL*8 R,N(20,1),EPSLN, TERM 

COMPLEX*8 B(20) 

TERM=0 .0D0 

DO 1 L=1, INDEX 

TERM=TERM+N(L, 1)*(CABS(B(L))**2.0) 
1 CONTINUE 

R=(EPSLN**4.0D0)*TERM*60.0D0 

R=R*0. 99930833D0 

RETURN 

END 
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** SUBROUTINE IMPROV ** 


THIS SUBROUTINE CALCULATES THE EXACT EIGENVALUES USING THE 
RECURSIVE HALF INTERVAL METHOD. FIRST IT COMPUTES AN ERROR 

BETWEEN A PREVIOUS EIGENVALUE(GAM) AND AN EIGENVALUE RETURNED 

FROM THE CONTINUED FRACTION(ZE OR ZO). ZE OR ZO IS OBTAINED BY 
PLUGGING GAM INTO THE CONTINUED FRACTION. POSITIVE OR NEGATIVE 
ERROR TELLS THE PROGRAM TO SEARCH DOWN OR UP. ONCE A POSITIVE 

ERROR AND NEGATIVE ERROR HAVE BEEN FOUND, CORRESPONDING TO TWO INPUTS 
TO THE CONT. FRAC., THE PROGRAM PECLARES THAT IT HAS FOUND THE 
INTERVAL WHERE THE ROOT(PRESUMABLY THE EIGENVALUE) OF THE CONTINUED 
FRACTION LIES. 

FROM HERE THE HALF INTERVAL METHOD TAKES OVER. THE HALF INTERVAL 
METHOD OPERATES BY CUTTING THE INTERVAL IN HALF, FIGURING OUT IN 
WHICH INTERVAL THE ROOT LIES AND CUTTING IT IN HALF AGAIN. THIS 
REPEATS UNTIL THE ERROR IS LESS THAN SOME TOLERANCE(TOLER). 

THIS SUBROUTINE IS CALLED BY MAIN, AND CALLS SUBROUTINE EVEN AND 
SUBROUTINE ODD. 

INPUTS: EPSLN, INDEX, GAM 

OUTPUT: GAMIMP 


SUBROUTINE IMPROV(GAM,GAMIMP,EPSLN, INDEX) 

REAL*8 GAM(20,1),Z0(20,1),ZE(20,1),AL,EPSLN, XL( 20,1) 
REAL*8 GAMIMP(20,1),ERROR( 20,1), TOLER, XM( 20,1) ,xU( 20,1) 
TOLER=0.00000001D0 


L ODD,(L-1) EVEN 


AL=1..0D0 

DO 1 L=1,(INDEX-1),2 

I=] 

J=] 

CALL EVEN(ZE(L,1),AL,L,GAM(L,1),EPSLN) 
ERROR(L,1)=GAM(L,1)-ZE(L,1) 
IF(ERROR(L,1).LE. TOLER. AND. ERROR(L,1).GE.-TOLER) GOTO 2 
IF(ERROR(L,1).LT.-TOLER) GOTO 3 


ERROR IS POSITIVE 


5 GAM(L,1)=GAM(L,1)-0.5D0 
I=I+¢l 
IF(1.GT.1000) GOTO 4 
CALL EVEN(ZE(L,1),AL,L,GAM(L,1),EPSLN) 
ERROR(L,1)=GAM(L,1)-ZE(L,1) 
IF(ERKOR(L,1)}.LE.TOLER.AND.ERROR(L,1).GE.-TOLER) GOTO 2 
IF(ERROR(L,1).GT. TOLER) GOTO 5 


FOUND THE INTERVAL 
WRITE(4,7) I,L 


7 FORMAT(1X,°FOUND THE INTERVAL” ,2X,215) 
GOTO 6 
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an 


ERROR IS NEGATIVE 


3 GAM(L,1)=GAM(L,1)+0.5D0 
I=I¢] 
IF(1.GT.1000) GOTO 4 
CALL EVEN(ZE(L,1),AL,L,GAM(L,1),EPSLN) 
ERROR(L,1)=GAM(L,1)-ZE(L,1) 
IF (ERROR(L,1).LE.TOLER.AND. ERROR(L,1).GE.-TOLER) GOTO 2 
IF(ERROR(L,1).LT.-TOLER) GOTO 3 


FOUND THE INTERVAL 
WRITE(4,7) 1,L 


HALF INTERVAL METHOD 


AOC eo 6: Ca«G 


XU(L,1)=GAM(L,1) 
XL(L,1)=GAM(L,1)-0.5D0 
GOTO 8 

6 XU(L,1)=GAM(L,1)+0.5D0 
XL(L,1)=GAM(L,1) 

8 XM(L,1)=(XU(L,1)+XL(L,1))/2.0D0 
J=J+1 
IF(J.GT.100) GOTO 14 
CALL EVEN(ZE(L,1),AL,L,XM(L,1),EPSLN) 
ERROR(L,1)=XM(L,1)-ZE(L,1) 
IF(ERROR(L,1).LE.TOLER.AND. ERROR(L,1).GE.-TOLER) GOTO 2 
IF(ERROR(L,1)}.LT.-TOLER) GOTO 9 
XU(L,1)=XM(L,1) 
GOTO 8 

9 XL(L,1)=XM(L,1) 
GOTO 8 


oO 


CONVERGENCE 


2 CONTINUE 
11 FORMAT(1X, “CONVERGENCE” ,2X,£12.5,1X,F20.15,1X,14) 
GAMIMP(L,1)=XM(L,1) 
AL=AL+2.0D0 
GOTO 1 


aa 


DID NOT FIND THE INTERVAL 


4 WRITE(4,12) 

12 FORMAT(1X,°DID NOT FIND THE INTERVAL’ ) 
AL*=AJ.+2.0D0 
GOTO 1 


aan 


NONCONVERGENCE 


14 WRITE(4,13) 

13 FORMAT(1X,°HALF INTERVAL METHOD DID NOT CONVERGE’ ) 
AL=AL+2.0D0 

1 CONTINUE 


L EVEN,(L-1) ODD 


OaqQura n 


AL=2.0D0 

DO 2] L=2,INDEX,2 

I=] 

J=] 

CALL ODD(ZO(L,1),AL,L,GAM(L,1),EPSLN) 
ERROR(L,1)=GAM(L,1)-ZO(L,1) 
IF(ERROR(L,1).LE.TOLER.AND.ERROR(L,1).GE.-TOLER) GOTO 22 
IF(ERROR(L,1)}.LT.-TOLER)} GOTO 23 
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C ERROR IS POSITIVE 


Cc 
25 GAM(L,1)=GAM(L,1)-0.5D0 
I=I+] 
IF(I.GT.1000) GOTO 24 
CALL ODD(ZO(L,1),AL,L,GAM(L,1) , EPSLN) 
ERROR(L,1) =GAM(L,1)-ZO(L, 1) 
IF(ERROR(L,1).LE.TOLER.AND.ERROR(L,1).GE.-TOLER) GOTO 22 
IF(ERROR(L,1).GT.TOLER) GOTO 25 
Cc 
C FOUND THE INTERVAL 
Cc 
Cc WRITE(4,27) I,L 
27 FORMAT(1X,°FOUND THE INTERVAL’ , 2X, 215) 
GOTO 26 
C 
C ERROR IS NEGATIVE 
Cc 
23 GAM(L,1)=GAM(L,1)+0.5D0 
I=I+] 
IF(I.GT.1000) GOTO 24 
CALL ODD(ZO(L,1),AL,L,GAM(L,1),EPSLN) 
ERROR(L,1)=GAM(L,1)-Z0(L,1) 
IF(ERROR(L,1).LE.TOLER.AND.ERROR(L,1).GE.-TOLER) GOTO 22 
IF(ERROR(L,1).LT.-TOLER) GOTO 23 
Cc 
C FOUND THE INTERVAL 
Cc 
Cc WRITE(4,27) I,L 
ec 
C HALF INTERVAL METHOD 
Cc 
XU(L,1) =GAM(L, 1) 
XL(L, 1) =GAM(L,1)-0.5D0 
GOTO 28 
26 XU(L,1)=GAM(L,1)+0.5D0 
XL(L,1)=GAM(L,1) 
28 XM(L,1)=(XU(L,1)+XL(L,1))/2.0D0 
J=J+1 
IF(J.GT.100) GOTO 214 
CALL ODD(ZO(L,1),AL,L,XM(L,1),EPSLN) 
ERROR(L,1) =XM(L,1)-ZO(L,1) 
IF(ERROR(L,1).LE.TOLER.AND. ERROR(L,1).GE.-TOLER) GOTO 22 
IF(ERROR(L,1).LT.-TOLER) GOTO 29 
XU(L,1)=XM(L,1) 
GOTO 28 
29 XL(L,1)=XM(L,1) 
GOTO 28 
Cc 
C CONVERGENCE 
C 
22 CONTINUE 
211 FORMAT(1X, “CONVERGENCE” , 2X, E12.5,1X,F20.15,1X, 14) 
GAMIMP(L,1)=XM(L,1) 
AL#=AL+2.0D0 
GOTO 21 
c 
C DID NOT FIND THE INTERVAL 
Cc 
24 WRITE(4, 212) 
212 FORMAT(1X,“DID NOT FIND THE INTERVAL’ ) 
AL=AL+2.0D0 
GOTO 21 
Cc 
C NONCONVERGENCE 
Cc 


214 WRITE(4, 213) 

213 FORMAT(1X,“HALF INTERVAL METHOD DID NOT CONVERGE’) 
AL#AL+2.0D0 

21 CONTINUF 
RETURN 
END 
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** SUBROUTINE ODD ** 


THIS SUBROUTINE CALCULATES THE CONTINUED FRACTION AS PER 
EQUATION 20,REFERENCE 1. TERM} AND TERM2 ARE THE FRACTION TERMS 
IN THE CONTINUED FRACTION. EVALUATION IS STARTED WITH THE 

LAST FRACTION TERM AND IS RECURSIVELY WORKED FORWARD UNTIL 

THE ENTIRE CONT. FRAC. IS EVALUATED. THE TWO CONTINUED FRACTIONS 
ARE THEN ADDED TOGETHER TO GIVE THE FINAL RESULT. 

THIS SUBROUTINE IS CALLED BY IMPROV AND DOES NOT CALL ANY 
SUBROUTINES. 

INPUTS: AL,L,GAM,EPSLN 

OUTPUT: 20 


SUBROUTINE ODD(Z0,AL,L, GAM, EPSLN) 
REAL*8 Al ,A2,A3,A4,Z0, TERM] , TERM2,AL, GAM, EP, EPSLN, ALM 
REAL*8 M 
M=1.0D0 
ALM=AL*2.0D0 
Al=201.0D0 
A2=200.0D0 
A3=199.0D0 
A4=198.0D0 
EP=EPSLN*EPSLN 
TERM1=0.0D0 
TERM2=0.0D0 
pO 1 I=1,100 
TERM! =(AL+A34M)* (AL+A4 )*EP/(A2*(ALM+A1 ) -GAM+TERM1 ) 
TERM2=(AL-A3+«M)* (AL-A2) *EP/(A2*(ALM-A3) +GAM+TERM2) 
Al=A1-2.0D0 
A2=A2-2.0D0 
A3=A3-2.0D0 
A4=A4-2.0D0 
1 CONTINUE 
Z0=TERM1-TERM2 
RETURN 
END 
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** SUBROUTINE EVEN ** 


THIS SUBROUTINE CALCULATES THE CONTINUED FRACTION AS PER 
EQUATION 19, REFERENCE 1. TERMI AND TERM2 ARE THE FRACTION TERMS 
IN THE CONTINUED FRACTION. EVALUATION IS STARTED WITH THE LAST 
FRACTION TERM AND IS RECURSIVELY WORKED FORWARD UNTIL THE ENTIRE 
CONTINUED FRACTION IS EVALUATED. THE TWO CONTINUED FRACTIONS 

ARE THEN ADDED TOGETHER TO GIVE THE FINAL RESULT. 

THIS SUBROUTINE IS CALLED BY IMPROV AND DOES NOT CALL ANY 

OTHER SUBROUTINES. 

INPUTS: AL,L,GAM, EPSLN 

OUTPUT: ZE 


SUBROUTINE EVEN(ZE,AL,L,GAM,EPSLN) 
REAL*8 Al ,A2,A3,A4,ZE, TERM3, TERM4,AL,GAM,EP,EPSLN,ALM 
REAL*8 M 
M=1 .ODO 
ALM=AL*2.0D0 
EP=EPSLN*EPSLN 
Al=201.0D0 
A2=200.0D0 
A3=199.0D0 
A4=198.0D0 
TERM3=0.0D0 
TERM4=0.0D0 
DO I I#1,100 
TERM3=(AL+A2+M)* (AL+A3) *EP/(A2*(ALM+ALl ) -GAM+TERM3 ) 
TERM4=(AL-A4+M)*(AL-A3 ) *EP/(A2*(ALM-A3 ) +GAM+TERMG ) 
Al=Al-2.0D0 
A2=A2-2.0D0 
A3=A3-2.0D0 
A4=-A4-2.0D0 
1 CONTINUE 
ZE=TERM3-TERM4 
RETURN 
END 
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APPENDIX B4. PROGRAM MONOSTOR 


THIS PROGRAM, IN FORTRAN LANGUAGE, CALCULATES THE INPUT 
IMPEDANCE 7 OF A 1/4 WAVELENGTH MONOPOLE ANTENNA MOUNTED 
ON A CIRCULAR GROUNOPLANF OF RADIUS A FOR E=2*PI+A/LAMBOA 
GREATER THAN APPROXIMATELY 30. WHERE LAMBOA IS THE FREE 
SPACE WAVELENGTH. 


REFERENCE: J.A.STORER, "THE IMPEOANCE OF AN ANTENNA OVER A 
LARGE CIRCULAR SCREEN” ,J.APPL.PHYS. ,VOL 22.N0O.8,PP. 1058-1066, 
AUGUST 1951,FQUATION PRECEDING EQUATION 20. 


LET FS-FREE SPACE WAVE IMPEOANCE 
R=RE(Z)=RAOIATION RESISTANCE 
X= IM(Z)=-REACTANCE 
ZINF-IMPEOANCE FOR F-INFINITY 
RINF-RAOIATION RESISTANCE FOR E-INFINITY 
XINF=REACTANCE FOR E=INFINITY 
K=DEPI/2)NUMERIC=ODIRECTIVE GAIN (NUMERIC) FOR 
FLEVATION ANGLE OF PI/2 RADIANS 
Kt=D(PI/2)DRI-OIRECTIVE GAIN (ORI) FOR 
FLEVATION ANGLE OF PI/2 RADIANS 
RO-=RAOLATION RESISTANCE OF MONOPOLE FOR E=0 
ON=DIRECTIVE GAIN AT AN ELEVATION ANGLE PI/2 
RADIANS OF MONOPOLF FOR E-O 
CHEE HOHE HEE EEE EHH HEHEHE HEHEHE HHH HEV HEHEHE HEHEHE HEHEHE HEHEHE THEE HH HEY 
REAL F(200) 
REAL K(200),K1(200),G(200) 
COMPLEX C1,C2,C3,C4,C5, ZINF ,7F(200),2(200) .C6 
FS$=376.73037 
PI-2. 141592654 
DO- 1.54259 
RO-19.43387 
WRITE (6,15) 
15 FORMAT (°1°./////7/) 
nO 10 M-1,150 
READ (5,20)E(M) 
20 FORMAI(F 15.7) 
10 CONTINUE 
M= 1 
AMY C1=CMPLX (0.,-FS/(4.+*PI+E(M))) 
C3-CMPLX (COS(2.'E(M)4+3.+*PI/4.),SIN(2.+*E(M)4+3.4PI/4. )) 
C2-CMPLX (COS(2.*E(M)),SIN(2. *E(M))) 
C4-C3/SART(4*PT+E(M)) 
C5=14+C4 
C6-=1/C5 
ZF(M)-C1*C2*CE 
ZiINF-CMPLX (36.539509,0. ) 
Z(M)-ZF(M)+ZINE 
G(M)-=REAL(ZD(M)) 
K(M)=DO+RO/G(M) 
K1i(M)=10. rALOGIOUK(M)) 
WRITE (6,25)E(M),ZF(M),7(M) KIM) K1(M) 
25 FORMAT (1X,7F 15.7) 
IF (M.G1.149) GO TO 50 
M=M+ 1 
GO 19 40 
50 STOP 
FND 
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I 


AANAANDAANAANAAANAAAANDANAANAANAAAAAAAAAAAAAANDNIYO ACY 


1 
A MONOPOLE ANTENNA AT THE CENTER OF A FINITE GROUNDPLANE", 
EEE TRANS. ON ANTENNAS AND PROPAGATION, VOL. AP-26 NO. 2. 
MARCH 1978. PP. 244-248 


PROGRAM AWADALLA 


THIS COMPUTER PROGRAM. IN FORTRAN LANGUAGE, 


WAS WRITTEN BY PROF. AWADALLA OF MANOUFTIA UNIVERSITY, EGYPT 

IT USES THE METHOD OF MOMENTS AND THE GEOMETRIC THEORY OF 
DIFFRACTION TO COMPUTE THE INPUT IMPEDANCE AND THE DIRECTIVITY 
PATTERN OF A MONOPOLE ELEMENT OF LENGTH h AND RADIUS b AT THE 
CENTER OF A CIRCULAR GROUNDPLANE OF RADIUS a. DETAILED 
DERIVATIONS ARE PUBLISHED BY IN REFERENCES 1 AND 2. 


THIS PROGRAM REQUIRES THREE INPUTS WHICH ARE ENTERED FROM 


AN INPUT FILE NAMED “AIN. DAT” WITH FORMAT 172. 


b = AL - MONOPOLE RADIUS IN WAVELENGTHS 
h = AH - MONOPOLE LENGTH IN WAVELENGTHS 


27a/hK= € = EPSLN —- NORMALIZED GROUNDPLANE RADIUS IN WAVENUMBERS 
THIS PROGRAM ALSO REQUIRES THE USE OF THE MMBSJ1 SUBROUTINE 
FUNCTION OF THE IMSL SOFTWARE LIBRARY. MMBSJ1 IS THE MODIFIED 
BESSEL FUNCTION OF THE FIRST KIND AND ORDER ONE. 


REFERENCE: 


K.H. AWADALLA AND T.S.M. MACLEAN, “INPUT IMPEDANCE OF 


2. K.H. AWADALLA AND T.S.M. MACLEAN, “MONOPOLE ANTENNA 
AT THE CENTER OF CIRCULAR GROUNDPLANE: INPUT IMPEDANCE 
AND RADIATION PATTERN”, IEEE TRANS. ON ANTENNAS AND 
PROPAGATION. VOL. AP-27, NO. 2, MARCH 1979, PP. 151-153 


NOTE: THE NUMBER OF MONOPOLE SEGMENTS USED IN THE PROGRAM 


IS 30. IF THE NUMBER OF MONOPOLE SEGMENTS IS CHANGED 
TO N. THEN THE DIMENSION OF ALL THE ARRAYS OF 30, AND 


30x30 SHOULD BE CHANGED INTO N. AND NXN. 


INTEGER IER 

REAL MMBSV1 

COMPLEX 2(30., 30), ZAZ(30, 30), B(30), ZC (30), FF, JHF (100) 
,EP,C1,C2,EPPM, EPMP, ZT, PSI, ZMG, YIN, ZIN, CDH. DH. EDP, SMP, Z2T2 
» SM2, CMEV, ZT1, JLD( 100), YAY 

DIMENSION IH(30), AUTDB(182), Y2(182), AUT(182), Y1(182) 
DIMENSION AUT2(182),. DIRECT(182)., DIRDB(182), Y(182) 
COMMON DL. RHO1.,PI,AK 

OPEN (UNIT=6, FILE=’AOUT ’, STATUS= ‘NEW’ ) 

OPEN(UNIT=5, FILE=‘AIN. DAT’. STATUS= ‘OLD’, READONLY ) 

DO 355 LIi=1,35 

N=30 

FMHZ=1000. O 

FREQ*FMHZ*1. OES 

ALAMDA=2. OE8/FREG 


C 3 Ot Mt St He Ht Mt Ht He HE ME He He ME Ot HE HE He He HE HE HE HE HG HE EE ee OE He EE EE OE ee a 


C 


INPUT READ STATEMENT at 


C 3 He He He He He HH He He He Ee He EE He EE He EE He He HE OE HE Ee He He He He 0 a a a a 


172 


990 


READ(S, 172) AL, AH, EPSLN 
FORMAT (3F 12. 8) 

WRITE(&, 990) 
FORMAT (1H1, “AWADALLA RESULT “> 
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*eeegnReREHRHE eH HER HEH HAH REHEKR ERA ERKR EHR EHR EH EH E 


2 Se SG ee AE aE aE HE aE Ee Ee aE aE aE ee ee Ee Ee ee ee he EE ee Ee ee Ee a eH 
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WRITE(46, 171) EPSLN, AH, AL 
171 FORMAT( 1X, ‘“EPSLN= ‘,1F 12.8. ° (NORMALIZED GROUNDPLANE RADIUS) ’, 
1 /1X, ‘“MONOPLE LENGTH= ’,1F12.8, ° (WAVELENGTH) ’ 
2 = 5s/1X, “MONOPLE RADIUS= °°, 1F12.8, ° (WAVELENGTH) ‘ ) 
R=AL#ALAMDA 
PI=4,. O#ATAN(1. 0) 
OM=2. O# PI#HFREQG 
AK=2. O#PI/ALAMDA 
AMU=4. O#PI#1, OE-7 
EPSI=1. 0/(36. O#PTI#1. OE9) 
DL=AH#ALAMDA/ (FLOAT (N) +40. 5) 
A=DL/2. O 
C1=CMPLX(0. 0, OM*#AMU#DL+*DL ) 
C2=CMPLX(0. O, -1. O/ (OM#EPSI ) ) 
RHO 1=EPSLN#ALAMDA/ (2. #PTI) 
DO 333 I=1.N 
B(I)=CMPLX(0.,.0. ) 
333 IHC(I)=0 
EP=PSI(0.0.R.,A, AK, PI) 
EPPM=PSI(DL.R, A, AK. PI) 
ZT=C1#EP+2. OX CQH% (EP—-EPPM) 
Z(1,1)2ZT 
RMNPM=DL 
NN=N+N 
DO 12 IN=2, NN 
EPMP=EP 
EP=EPPM 
RMNP M=RMNPM+ DL 
EPPM=PSI(RMNPM,.R, A, AK, PI) 
ZT=C2*( EP+EP-EPMP-EPPM ) 
IFCIN. GT.N)GO TO 13 
Z(1, IN) =C1#EP+ZT 


GO TO 12 
13 IC=IN-N 
ZCCIC)SCixEP+ZT 
12 CONTINUE 


DO 16 J=1,N 
DO 11 I=sJ,N 


TA=I-J+1 

Z(J, T)2Z01, TA) 
11 ZCI, J)@Z(U, 71) 
16 CONTINUE 


DO 22 K=2.,N 
DO 20 I=1,N 
IB=I+K 
IF (CIB. GT.N)GO TO 21 
Z(K, IT) SZCK, 1)94+Z(1, IB) 
GO TO 20 
ai IC=IB-N 
Z(K, I) =Z(K, 1) *#ZC CIC) 
20 CONT INUVE 
a2 CONT INVE 
NF =N~-1 
DO 200 I=1,NF 
IZ=I+1 
200 Z(¢1,1)2Z¢1, 194Z(1,1Z) 
Z(1,N)*Z¢1,.N)4ZC(1) 
DO 121 I=1,.N 
DO 122 J=#1,N 
ZAZ(I, J) @Z(1I, J) 
122 CONTINUE 


203 


AWADALLA (3 of 8 pages) 


121 CONT INUVE 
CH a He IE HE IE IE HE IE HE HE IE HE EI EE HE HE IE TE HE IE IE EE IE EE EH SEE HE IE IE HE IE IE EI IE HE IE I HEE SE aE HE Ha 
Cc DIFFRACTION EFFECT OF THE EDGE OF GROUNDPLANE * 


CE HE ae He I EE IE HEI HE IE IE HE aE I EE HE EE IE EEE I EE I EE a a EE IE HEE a aE Ea 
DO 30 I=1,N 
ZS=DL*(FLOAT(I)-O. 5) 
DO 32 J=I,N 
ZO=DL* (FLOAT (U)-O. 5) 
Z(I, J)=Z¢01, J)+ZMCCZS, ZO) 


32 Z(J,1T)*2Z¢1, J) 
30 CONT INUVE 
CHa He a He aE I HE aE IE I aE ea I aE a EI IE Ee I aE a EE a EE eI I HE Ie aE 
C CALCULATION OF THE ADMITTANCE * 
C99 9 A ae aE aE aE aE EE He ae aE a aE a EEE EE SE EE I EE EE JE EE I EEE HE EE a Hd 
150 CALL MTXINV(Z.N. B. IH) 
151 ZIN=1.0/Z(1,1) 
WRITE (6, 23)ZIN 
23 FORMAT(1X, ‘INPUT IMPEDANCE= ’,E15.9, ’ +yjy ‘,E15.9, 2X, ’COHMS) ‘//) 
WRITE (4, 999) 
99 FORMAT (/1X,. ‘ELEVATION’, 53X, ‘ELEVATION’, /1X, 2¢ 


1 ‘ANGLE’, 6X, ‘DIRECTIVITY’, 5X, ‘DIRECTIVITY’, 5X, 
2 ‘RELATIVE POWER’, 5X),/1X,2(¢’ (DEG) ’, 6X, ’(NUMERIC) ‘, 
3S 7X, CdB4) 4.11%, CB) %. 15%) 
CHM HEHE EEE 
C COMPUTING THE E FIELD RADIATION PATTERN + 
CPE Leeee eee etree PES SPT SE PEE ET EPC TC PPC PC CPT PPP PPT 
CDH=-CMPLX(1.,1. )/SGRT(2. #PI) 
AKA=AK#RHO1 
PID2=PI/2. 
PI4=4. PI 
PID4=PI/4. 
SAKP=SGRT(AK/ (PI#8. )) 
DTH=P1/180. 
TH=O. 
RHO2=RHO1#RHO1 
ZT=CMPLX(0.,0. ) 
NG=RHO1/(0. OS#ALAMDA) 
DR=RHO1/NG 
RG=DR/2. O 
CJUJHF=-2. #DL/PI4 
DO 40 IK=1,NG 
ZT1=CMPLX(O.,0. ) 
ZT=CMPLX(O.,0. ) 
H=A 
DO 41 J=1.N 
RG2=RG#RG 
H2=H#*H 
RA2=RG2+H2 
RA=SGRT(RA2) 
SNA=RG/RA 
AKRA=-AK#RA 
CSAKRA=COS (AKRA) 
SNAKRA=SIN(AKRA) 
ZT=ZT+Z(1, J) *(SNA/RA) #CMPLX (CSAKRA, SNAKRA) #CMPLX (1. O/RA, AK) 
RI2=RHO2+H2 
RI=SQGRT(RI2) 
AKRI=-AK#RI 
AK1=-AK# (RHO1-RG) 
ZT1=ZT1+Z(1, J) #RHO1*DL#CMPLX (COS(AKRI), SINCAKRI) ) *#CMPLX(1. /RI, AK) 
1 *CMPLX(COS(AK1), SINCAK1) )/(PI4#RI) 
H=H+DL 
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41 CONTINUE 
JHF (IK) =ZT#CUHF+ZT1#2. O 
YLD( IK) =ZT1#2. 0 
RG=RG+DR 
40 CONTINUE 
FCR=RHO1#SQGRT (2. #*AK) /PI*DL 
DO 33 11, 181 
SNTH=SIN(TH) 
CSTH=COS( TH) 
H=A 
EDP=CMPLX(0O.,0. > 
DO 34 J=1,N 
S=SGRT (H#H+RHO2 ) 
PHI=ATAN(RHO1/H) 
SB=PI-PHI 
RB=P ID2-PHI 
AIl=1.0 
100 ZAC =AK#H#CSTH 
EDP=EDP+Z (1, J) *SNTH#DL#AI #CMPLX (COS(ZKC),SIN(ZKC> > 
34 H=H+DL 
THA=TH#*180. O/PI 
ZT=CMPLX(O.,0. ) 
RG=DR/2. 
CONST=~2. #PI1#DR#*CSTH 
AKSN=AK#SNTH 
DO 42 K=1,NG 
X=4KSN#RG 
IF(X.LT.O. 1E~6)CG0 TO 44 
SABF=MMBSVJ1(X, IER) 


GO TO 45 
44 SABF=O. 0 
45 ZT=ZT+JHF (K) #SABF#CMPLX (0. , CONST) #RG 
47 RG=RG+DR 
42 CONTINUE 
43 EDP=EDP+ZT 
AEDP=CABS (EDP ) 


AUT (I }=AEDP 

AUT2( 1) =AEDP #AEDP 

IF (AEDP.EQ@.0.0) GO TO 33 

PHASE=180. O#ATAN2(AIMAG (EDP), REAL (EDP) )/PI 


33 TH=TH+DTH 
X=0. O 
DO 70 K=1, 181 
CHO=AUT (K) 

70 IF (CHO. GT. X) X=CHO 


DO 71 I#1,181 
YCI)SAUT(C I> /X 

YiICT =YCT)#YCT) 
IFCY1(1).EQ@. 0. > THEN 
Y2(1)=999 

AUTDB (1 )=999 

ELSE 

AUTDB( I) =1LO#ALOGCIOCAUT(I)) 
Y2CI)=1LO#ALOGIOCY1I¢1)) 


END IF 
71 CONT INUVE 
C96 46 0 aE a aE a ba a a a a a a a a a a a a a a OO OO 
Cc COMPUTE THE DIRECTIVITY PATTERN * 
C69 4 a a Ab a A a a HE a a HH a a a a aH a a aH aH HH 
PR=0.0 


DO 995 I=2, 181 
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995 


994 


997 


996 


998 


DEL=(I-1)*#DTH 

PR=PR+AUT2(1)#SIN(DEL)*#DTH 

CONTINUE 

DO 994 I=1, 181 

DIRECT(I )=2#AUT2(1)/PR 

IF (DIRECT(I). EG. O. ) THEN 

DIRDB(I)=999 

ELSE 

DIRDB( I )=10#ALOGIO(DIRECT(I)) 

END IF 

CONT INVE 

DO 998 I=1, 45 

Il=I-1 

Il 2=I+43 

II3#II2+1 

IF (I. EQ. 1) THEN 

WRITE (6, 997) II, AUT(T) 

FORMAT (2X, I3, 8X, F9. 5, SX, ‘- INFINITE‘, 6X, ‘-~INFINITE “ ) 
ELSE 

WRITE(6, 996) II, DIRECT(I), DIRDB( I), Y2(1I), II2, DIRECTC(II3) 
»DIRDBC(IIS), Y2C(II3) 

FORMAT( 1X, 2(¢1X, 13, 8X, FI. 5, SX, FF. 5, OX, FO. 5, 1X) ) 
END IF 

CONTINUE 

WRITE (6, 990) 

WRITE(6,. 171) EPSLN, AH, AL 

WRITE(6, 23)ZIN 

WRITE (6, 999) 

DO 992 1290, 136 

IIl=I-1 

II2=I1+44 

II3=IIl2+i 

WRITE(6, 996) 11, DIRECT(I), DIRDB(T), Y2(I), Ila, DIRECT(II3) 
»DIRDBC(II3), Y2C(II3) 

CONTINUE 

CONT INUE 

CLOSE (UNI T=6, STATUS= ‘KEEP “ ) 

CLOSE (UNIT=5, STATUS= ‘KEEP ’” ) 

STOP 

END 
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C Ht at Ht Ht Ht Ht TE A aE aE Ht HE HE aE OEE a OE HEE aE SEE a EH EE Se aE SE aE a SE SE ESE EE SE SE SE a a OE 2 aE aE EE a a 
Cc SUBROUTINE ZMG & 
C tH Ht HEH aE HE EE HE EE IE HE SE HE HE HE HE OE HE HE OEE HE SHE HEE HE HE HE HE EE HEE EH HE HE SHE HE EE HH a OE HE a a 

COMPLEX FUNCTION ZMG(ZS, ZO) 

COMMON DL. RHO1, PI. AK 

ZO2=Z0*Z0 

ZS2=ZS*7ZS 

RH2=RHO1#RHO1 

RS2=RH2+ZS2 

RO2Z2=RH2+7Z02 

RS=SQRT(RS2) 

RO=SQRT(RO2) 

AKA=AK# (RS+RO) 

PS=ATAN(ZS/RHO1 ) 

PO=ATAN(ZO/RHO1 ) 

ZMG=15. *DL#DL#RHO1 #RH2/ (AK#RS2#RO2Q) #CMPLX(O.,-1. )# 

1 CEXP(CMPLX(O. , -AKA))#(1. /COS((PO-PS)/2. >+1. /COS( (PO+PS) 7/2. )) 
2 *CMPLX(1. /RS, AK) *CMPLX (1. /RO, AK) 

35 RETURN 

END 


C2646 96 FE aE HE aE SE HE EE SE IEE Ee ab ae a ae a aE ae ae IE ae IE IE aE HE EE ae a He a ae a ae aE a HE HE a OE a HE HE HO OE OE OE a EE 


Cc SUBROUTINE PSI + 
C ra 
Cc THIS FUNCTION IS TO CALCULATE THE SELF AND MUTUAL + 
Cc IMPEDANCES FOR A STRAIGHT WIRE ANTENNA + 


C3 3b HH HE HE IE I aE HE a HE ae SE aE EH He HE SE HE HE SE EE SE EE aE HE a aE HH aE aE SE EH I HE EH SE OE HE HE a HH Ha Ht 
COMPLEX FUNCTION PSI(CH,R., A, AK. PI) 
IMPLICIT REAL#¥8 (A-H), (P-Z) 
COMPLEX#146 P1,P2,P3,P4,PII 
R2=R#R 
D=DSGRT (H#H+R2 ) 

S1=DSGRT (R2+(H+A) ##2) 

S2=DSGRT (R2+(H-A) ##2) 

AI1=DLOG( (H+A+S1) /(H-A+S2) ) 
AI2=2. #A 

AI3=(A+H) /2. #S14+(A-H) /2. #S24+R2/2. #AI1 
AI4=2. *A#R2Q4+(2. HAHAHATG. HAHHHH) /3. 
AKA=-AK#D 

PIA=6. #*PI#A 

P4=DCMPLX (0. ODO, AKA) 
P1=CDEXP(P4)/PIA 

PI=AI1-(AI2Z—-D#AI1 )*#DCMPLX(O. ODO, AK) 
P=(AK#AK/2. )#(AI3Z—-2. #D#AIQ+AI1#D#D) 
P3=DCMPLX(0. ODO, AK#AK#AK/4. ) 
PP=AI4-3. #D#AI343. #*D#D#AI2 
PP1=D#D#D+#AI 1 

PS5=PP-PP 1 

PII=P1# (P2-P+P3#P5) 

PSI=CMPLX(PII) 

RETURN 

END 
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OPEL EE SEES LESS eT ETP PL EP TP PTT TT PPT PT TT PP 


C 
C 
C 
C 


SUBROUTINE MTXINV 


* 
% 


THIS SUBROUTINE IS FOR THE INVERSION OF A MATRIX AC(N,N) # 


THE MATRIX SHOULD BE A SQUARE ONE OF ORDER 


Ag 


+ 


C 3 Ft HE FE HE IE EE ETE EE IE IE HE HE IE HE He FEE IE HE HE IE HE HEE IE He EEE IE IE IE EE IE HE HE EE TEE EE EE 


24 


12 


11 


22 


14 


SUBROUTINE MTXINV(A, N. B, IH) 
COMPLEX A(N, N),.BC(N),@ 
DIMENSION IH(N) 

DO 24 I#1,N 

BCI )=CMPLX(O. 0,0. 0) 
IH(1)=0 

DO 1 K=1,N 

IK=0 

Z=CABS(ACK, K) ) 
IF(Z.E@.0.0) GO TO & 
Q@=1. O/A(K, K) 


GO TO 7 

IF(K. E@.N)GO TO 13 
IH(K) =K 

DO 8 I=1,N 
BCI)=ACK, I) 

NA=N-1 

DO 9 I=K,NA 

II=I+1 


DO 10 J=1,N 

ACT, JI=ACTI, J) 

CONT INVE 

DO 11 J=1,N 
ACN, J) =BCU) 

IK=IK+1 

IFCIK.EQ@.1) GO TO 12 
WRITE (6, 22)K,K 


FORMAT(1X, ‘MTXINV FAIL BECAUSE ELEMENT Z(’, T2, ’%. ’, Ta, 
‘) AND THE NEXT IN THE SAME COLUMN ARE ZEROS’/) 


GO TO 23 

DO 2 I™1,N 
IF(I.E@ K) GO TO 2 
DO 3 J=1,N 

IF (J. EQ@.K)GO TO 3 
ACT, JI=ACT, JI-ACT, KD #ACK, J) #Q 
CONTINUE 

CONT INVE 

DO 4 I=1,N 
IF(IT.€Q@.K) GO TO 4 
ACT, K)#ACT, KA) #Q@ 
CONT INVE 

DO 5S I=1,N 

IF CI. E@. K) GO TO 5 
ACK, I) #-ACK, 1) *Q 
CONT INVE 

A(K. K)=@Q 

CONT INVE 

N1=N+1 

DO 15 I=1,N 
NN=N1-I 
IFCIHCNN).EQ@.0) GO TO 15 
DO 16 J#1,N 

BCU) =ACN, J) 
K=TH(NN) 

KK=N-K 
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DO 17 J=1, KK 


N2=N1-J 
NS=Ne2-1 
DO 18 L=1,N 
16 A(N2,L)=A(N3,L) 
17 CONT INVE 
DO 19 L=1,N 
19 ACK, LISBC(L) 
15 CONT INVE 
GO TO 20 
13 WRITE (6, 14) 
14 FORMAT( 1X. ‘AC(N,N)#0O: PROGRAM FAILIOUR. ‘) 
20 RETURN 
23 STOP 
END 
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